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SECTION  I 


INTRODUCTION 

The  study  of  ice  nucleation  has  been  a  topic  or  ’nterest  for  many 
years  [1-4].  Most  experimental  work  to  date  has  involved  measurements 
of  bulk  macroscopic  properties  (such  as  temperature  and  pressure).  Con¬ 
troversy  arose  among  proponents  of  competing  theories  designed  to  ex¬ 
plain  the  results  of  these  measurements  [5].  At  the  root  of  the  contro¬ 
versy  is  the  disagreement  over  the  microscopic  details  on  which  the 
macroscopic  models  are  based.  The  wide  discrepancy  among  the  predic¬ 
tions  of  these  models  suggests  the  need  for  experiments  from  a  micro¬ 
scopic  (molecular)  point  of  view.  The  original  goal  of  the  present 
study  was  to  design  and  perform  such  experiments  to  study  heterogene¬ 
ous  ice  nucleation  (nucleation  on  a  substrate  other  than  ice  itself). 

A  molecular  beam  scattering  system  seemed  the  ideal  sort  of  tool 
to  use  in  such  a  study.  The  fate  of  molecules  with  particular  incident 
velocities  and  directions  can  be  studied,  and  one  can  approach  the 
conditions  for  the  creation  of  stable  nuclei  arbitrarily  slowly  by  con¬ 
trolling  beam  intensity  and  substrate  temperature.  Some  progress  has 
been  made,  but  thus  far,  we  are  still  unable  to  correlate  microscopic 
observations  with  macroscopic  observations  in  any  detail.  Along  the 
way,  there  have  been  some  interesting  additions  to  the  knowledge  of  homo¬ 
geneous  nucleation  in  nozzle  beams,  evaporation  from  surfaces,  and, 
perhaps  most  importantly,  the  interaction  of  ice  clusters  with  surfaces. 

It  is  not  practical  to  give  a  complete  survey  of  work  on  ice  nucle¬ 
ation;  to  do  so  could  fill  a  sizeable  book.  Substantial  portions  of 
books  by  Mason  [1],  Hobbs  [2],  and  Abraham  [3],  for  example,  are  devoted 
to  ice  nucleation.  These  books  are  already  aging  by  modern  scientific 
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standards,  but  no  major  breakthroughs  have  occurred;  Hobbs  [2]  can  still 
be  used  for  a  reasonable  picture  of  the  current  state  of  knowledge  about 
all  aspects  of  the  physics  of  ice,  and  Abraham  [3]  gives  a  good  survey 
of  homogeneous  nucleation  in  general.  A  more  recent  survey  of  homoge¬ 
neous  nucleation  has  been  written  by  Springer  [4j,  whicn  includes  HO 
references  to  work  in  the  field.  We  restrict  discussion  here  to  work 
which  takes  a  molecular  point  of  view,  and  review,  in  addition,  some 
work  on  the  experimental  methods  involved  and  some  closely  related  work 
involving  species  other  than  water. 

1.1.  ICE  NUCLEATION  FROM  A  MOLECULAR  POINT  OF  VIEW 

To  date  relatively  little  work  on  ice  nucleation  has  been  done  from 
a  molecular  point  of  view,  and  a  large  share  of  what  has  been  done  is 
theoretical.  Most  of  this  work  has  been  concentrated  in  a  single  loca¬ 
tion:  the  Center  for  Cloud  Physics  at  the  University  of  Missouri  at 
Rolla.  Numerous  theoretical  papers  by  Kassner,  Plummer ,  Hale  and  others 
and  some  experimental  papers  by  Levenson  and  others  have  appeared  and 
are  discussed  below. 

1.1.1.  THEORETICAL  WORK 

Several  sorts  of  calculations  have  been  performed.  Hale  and  Plum¬ 
mer  first  constructed  molecular  models  of  small  clusters  of  molecules 
of  water  [6]  for  comparison  with  the  liquid  drop  theory  of  nucleation, 
which  is  built  up  from  consideration  of  surface  and  bulk  energies  of 
individual  clusters.  Such  a  theory  is  suspect  for  small  clusters  since 
the  extension  of  the  concept  of  surface  energy  down  to  clusters  of  a 
few  molecules  is  unrealistic.  (See  the  controversy  over  the  "replacement 
free  energy"  term  necessary  (?)  to  correct  the  model  to  agree  with 
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experimentally  observed  nucleation  rates  [5].  In  the  theory  of  Hale 
and  Plummer,  free  energy  of  formation  of  the  clusters  was  calculated, 
and  from  that,  macroscopic  quantities  (critical  supersaturation  ratio, 
critical  nucleus  size,  nucleation  rates)  could  be  calculated.  Most 
notable  in  the  results,  besides  the  small  but  significant  deviations 
from  the  results  of  the  liquid  drop  theory,  was  the  increased  stabil¬ 
ity  of  clathrate  structures  (see  Section  4.1)  over  structures  like  bulk 
ice. 

Statistical  mechanics  has  been  applied  to  various  model  systems 
which  include  models  of  growing  clusters  [7-10].  There  are  problems 
with  such  an  approach  at  both  ends  of  the  size  scale.  First,  a  satis¬ 
factory  statistical  model  requires  accurate  knowledge  of  the  possible 
microscopic  configurations  and  the  formation  energies  involved.  Such 
knowledge  is  still  actively  being  sought.  Boundary  conditions  at  the 
edges  of  a  system  must  be  established  though  the  choices  made  with  con¬ 
venience  of  calculation  in  mind  may  be  unrealistic.  A  kinetic  approach 
is  necessary  to  handle  such  situations  as  the  loss  of  mass  as  particles 
grow  big  enough  to  fall  out  of  a  system  (snow).  Statistical  mechanics 
is  most  suited  to  the  calculation  of  equilibrium  distributions,  which 
for  the  case  of  a  system  allowing  cluster  growth  generally  ends  up 
with  all  the  mass  in  a  single  large  cluster.  The  fact  that  such  a  state 
may  not  be  achieved  in  finite  time  or  that  other  vastly  different  states 
may  be  nearly  as  favorable  is  ignored.  Quantum  statistical  mechanics 
shows  some  promise  of  getting  closer  to  first  principles,  since  it  does 
not  require  any  ad  hoc  input  of  microscopic  details  [11]. 

Simulation  of  a  nucleating  system  by  "Molecular  Dynamics" 
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calculations  has  recently  been  performed  for  a  two-dimensional  100  mole¬ 
cule  system  by  Zurek  [12,  13].  McGinty  has  performed  a  molecular  dynam¬ 
ics  study  on  argon  clusters  [14],  and  Abraham,  Binder  and  others  [15-18] 
have  performed  some  Monte  Carlo  simulations  of  water  clusters.  The 
approach  provides  useful  insight  into  the  processes  which  may  occur  in 
a  nucleating  system  but  with  current  computer  technology,  computing  costs 
are  prohibitive  for  all  but  the  simplest  systems.  This  situation  may  be 
expected  to  improve  rapidly  for  at  least  a  couple  of  decades,  and  simu¬ 
lation  of  more  and  more  realistic  systems  should  become  possible. 

All  the  calculations  mentioned  thus  far  have  been  for  homogeneous 
nucleation  systems  typically  containing  only  a  single  molecular  species. 
Work  has  begun  on  theoretical  modeling  of  heterogeneous  nucleation.  Sta¬ 
tistical  mechanics  has  been  applied  to  a  system  including  a  surface, 
clusters  on  the  surface,  and  clusters  and  monomers  in  a  vapor  above  the 
surface  [19].  Results  which  compare  satisfactorily  with  experiments 
have  not  yet  been  achieved,  but  the  method  shows  promise.  Again  it  is 
the  microscopic  details  which  cause  problems.  Quantum  mechanical  model¬ 
ing  of  the  interaction  of  a  monomer  with  a  simulated  surface  (restricted 
to  10  atoms  to  make  the  calculation  manageable)  [20]  has  begun  to  shed 
some  light  on  exactly  where  bonding  is  favorable.  Such  an  approach  may 
shed  light  on  microscopic  details,  though  as  yet  calculations  for  anything 
but  the  simplest  systems  are  formidable. 

1.1.2.  EXPERIMENTAL  WORK 

It  is  not  possible  by  current  experimental  techniques  to  measure  the 
state  of  individual  water  clusters,  but  one  can  infer  the  behavior  of 
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single  clusters  more  readily  from  molecular  beam  measurements  than 
from  any  others  now  available.  Measurements  on  beams  themselves  are 
described  in  Section  1.2.  Very  little  work  on  heterogeneous  nuclea- 
tion  has  been  attempted  to  date.  The  few  exceptions  include  the  work 
of  Bentley  and  Hands,  who  studied  condensation  of  argon,  nitrogen,  and 
carbon-dioxide  on  gold  by  monitoring  reflected  flux  with  a  quadrupole 
mass  spectrometer  [21].  They  found  no  nucleation  barrier  for  the  con¬ 
densation  of  argon  and  nitrogen,  and  a  barrier  corresponding  to  a  crit¬ 
ical  nucleus  size  of  about  9  molecules  for  carbon  dioxide.  Nucleation 
apparently  took  place  on  an  adsorbed  layer  rather  than  on  the  bare  gold 
surface.  Similar  studies  were  performed  by  Levenson  et  al  [22,  23]  using 
a  water  beam,  gold  and  silver-iodide  surfaces  and  a  quartz  crystal  micro¬ 
balance  as  a  flux  detector.  Activation  energies  for  adsorption  and 
desorption  were  obtained,  and  a  critical  cluster  size  of  only  one  or  two 
molecules  was  found.  They  also  studied  carbon  dioxide  nucleation  on 
ice  [24]  finding  a  critical  nucleus  size  of  about  4  molecules.  Such 
small  critical  nucleus  sizes  (compared  to  that  for  homogeneous  nuclea¬ 
tion)  is  surprising,  but  agrees  with  that  calculated  for  the  epitaxial 
growth  of  metals  from  the  vapor  [25]. 

1.2.  CONDENSATION  IN  BEAMS 

Condensation  (clustering)  in  beams  has  become  a  widely  studied 
topic.  In  many  situations,  such  as  the  initially  envisioned  plan  for 
the  present  study,  it  is  desirable  to  avoid  condensation  in  the  beam 
and  to  have  only  a  high  intensity  monomer  beam.  In  high  pressure  nozzle 
sources  of  all  types  (whether  "sonic”  ("free  jet";  converging  nozzle 
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or  small  hole)  or  "supersonic"  (converging-diverging  nozzle)',  conden¬ 
sation  occurs  readily  for  most  species  because  of  the  extreme  cooling 
which  takes  place  during  the  expansion.  Applications  of  condensed  beams 
include  fusion  devices  [26],  film  deposition  [27],  and  the  study  of 
homogeneous  nucleation  (see  below).  A  good  general  picture  of  conden¬ 
sation  beams  can  be  had  by  referring  to  recent  proceedings  from 
symposia  on  Rarefied  Gas  Dynamics  [23]  which  included  sessions  devoted 
to  “Condensation  in  Expansions”. 

Condensation  has  been  observed  in  beams  of  Ar,  ,  C02[28,  29], 

Hg,  02,  02,  CO,  SF6,  UFg,  [28],  N02,  Kr[29] ,  CgHg,  CHCI3,  CC1 3F ,  CgHgOH 
[30],  CH^OH  [31],  H20  [28,  31-37].  Some  of  these  studies  used  free  jet 
expansions  and  some  used  supersonic  nozzles.  Obert  and  Hagena  [38,  39] 
have  studied  extensively  the  effects  of  various  parameters  of  super¬ 
sonic  nozzles  on  clustering.  Dankert  and  Koppenwallner  [40]  have  studied 
the  properties  of  clustering  in  free  jet  expansions.  These  results  are 
of  interest  from  the  point  of  view  of  designing  new  sources  (see  Section 
2.2),  but  unless  a  source  can  be  exactly  duplicated,  it  is  necessary 
to  determine  the  properties  of  each  new  source  to  the  extent  that  those 
properties  are  relevant  to  any  particular  experiment.  The  theory  of 
nozzle  expansion  is  not  sufficiently  developed  to  predict  accurately 
the  behavior  of  an  arbitrary  nozzle,  especially  if  there  is  the  possi¬ 
bility  of  condensation. 

A  variety  of  techniques  have  been  developed  for  the  study  of 
clusters  in  beams.  The  simplest  measurements  to  make  are  measurements 
of  pressure  (density)  or  flux  [39].  Discontinuities  in  these  measurements 


6 


as  a  function  of,  say,  source  pressure  have  been  interpreted  [39]  as  an 
indication  of  massive  condensation  in  the  beam.  Similar  but  more  sen¬ 
sitive  results  can  be  obtained  by  scattering  laser  light  [41].  The 
laser  axis  is  along  the  beam  axis,  and  the  scattered  light  is  measured 
perpendicular  to  the  beam.  It  is  also  possible  to  scatter  an  atomic 
beam  from  a  cluster  beam  [42],  but  the  results  are  more  difficult  to 
interpret.  Mass  spectrometer  techniques  (quadrupole  [32-34],  retarding 
field  [39]  or  ion  time-of-fl ight  [43])  show  the  most  promise  for  actu¬ 
ally  measuring  cluster  size  distributions,  but  instruments  with  suffi¬ 
cient  mass  range  to  study  clusters  of  thousands  of  molecules  generally 
have  very  poor  mass  resolution  and  severe  amplitude  calibration  prob¬ 
lems  (see  Section  3.1.1.).  High  energy  (50  keV)  electron  diffraction 
measurements  have  been  used  [35-37]  to  determine  crystalline  structure 
and  cluster  size  (from  diffraction  line  widths).  Electron  diffraction 
is  at  present  limited  to  the  study  of  clusters  of  more  than  ~300  molecules. 

1.3.  CLUSTER  SCATTERING 

As  yet  there  have  been  only  a  few  experiments  on  the  scattering  of 
clusters  from  surfaces.  Gspann  and  Krieg  [44]  have  studied  the  scatter¬ 
ing  of  helium,  hydrogen,  and  nitrogen  clusters  from  a  stainless  steel 
surface.  The  surface  was  carefully  polished  (the  measured  roughness 
depth  was  less  than  500  &  and  cleaned  by  baking  to  450°  K.  Their 

measurements  included  mass  spectrometry  and  crude  velocity  analysis. 

4  5 

The  beams  were  highly  clustered  consisting  of  clusters  of  10  to  10 
molecules.  A  grazing  incident  angle  (84.3°)  was  used.  Two  sorts  of 
scattering  were  observed.  For  hydrogen  on  a  target  above  215°  K  and 
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and  for  helium  at  all  target  temperatures  tested  (SO0  K  <  T  <  550°K) 
clusters  gained  velocity  by  a  few  per  cent  and  lost  mass.  The  reverse 
was  true  for  hydrogen  at  a  lower  target  temperature  and  for  nitrogen. 
Superspecular  scattering  occurred  for  and  subspecular  occurred 
for  He.  Nitrogen  clusters  apparently  partially  fragmented  on  collision, 
breaking  into  two  pieces.  Reflection  coefficients  (fraction  of  clusters 
surviving  collision)  were  around  0-0.6  for  He,  0. 4-1.0  for  Hj.  and 
0-0.1  for  N^.  (The  monomers  resulting  from  destruction  of  clusters 
were  not  studied.)  The  exceptionally  high  reflection  coefficient  for 
H 2  occurred  at  the  so-called  “optimal  reflector  temperature"  which 
divided  the  two  regions  of  scattering  behavior. 

No  previous  experiments  (to  our  knowledge)  have  studied  the  scatter¬ 
ing  of  water  clusters  from  surfaces.  We  describe  here  the  design  and 
construction  of  a  system  for  such  a  study,  and  give  a  report  on  experi¬ 
ments  which  we  have  performed  thus  far.  Detailed  theoretical  modeling 
is  incomplete,  and  will  require  some  additional  experimental  work. 
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SECTION  II 


EXPERIMENTAL  APPARATUS  AND  DATA  HANDLING  TECHNIQUES 

2.1.  PRIMARY  VACUUM  FACILITIES 

The  experimental  chamber  was  a  1.2  m  diameter  viton-sealed  stain¬ 
less  steel  vacuum  chamber  pumped  by  a  61  cm  oil  diffusion  pump  with  a 
freon-cooled  baffle.  Additional  pumping  for  water  and  other  condensa¬ 
bles  was  achieved  with  strategically  located  1 iquid-nitrogen-ccoled 

surfaces.  Pumping  speed  for  water  was  more  than  10®  1/s  versus  between 
3  4 

10  and  10  1/s  for  most  gases.  This  system  has  been  described  in  some 

detail  elsewhere  [45].  The  base  pressure  achievable  was  about  5  x  10~7  Pa 
with  a  typical  operating  base  pressure  of  around  4  x  10"®  Pa  (there  were 
no  provisions  for  bake-out,  so  pumpdown  to  5  x  10’7  Pa  took  weeks). 

Irv  a  typical  scattering  experiment,  a  rotatable  and  moveable  target 
was  mounted  in  the  middle  of  the  chamber  with  other  equipment  mounted 
in  various  positions  around  it  (see  Figure  1).  Other  equipment  included 
a  source  system,  a  time-of-fl ight  detection  system,  a  stagnation  gauge, 

»nd  a  quadrupole  mass  spectrometer.  These  are  described  separately  be¬ 
low. 

2.2.  SOURCE  SYSTEMS 

Early  measurements  were  made  using  the  source  assembly  shown  in 
Figure  2.  There  was  a  brass  nozzle  with  a  cylindrical  hole,  .08  cm  in 
diameter  by  .66  cm  long.  The  gas  was  supplied  either  from  an  external 
manifold  system  or  an  internal  water  reservoir.  A  motor  driven  needle 
valve  between  the  reservoir  and  the  nozzle  functioned  primarily  as  a 
shutoff  valve.  The  nozzle  was  located  3  cm  from  a  collimating  plate 
in  the  wall  of  a  small  liquid  nitrogen  cryopumped  collimating  chamber. 
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This  source  provided  the  beam  for  most  of  the  calibration  work 
on  the  time-of-fl ight  detection  system.  It  proved  unsatisfactory  for 
scattering  experiments  for  several  reasons.  The  liquid  nitrogen  pumped 
collimating  chamber  functioned  satisfactorily  only  for  a  water  beam 
(other  gases  used  were  not  condensable).  Even  for  a  water  beam,  the 
geometry  was  such  that  insufficient  pumping  speeds  were  achieved,  espec¬ 
ially  at  higher  beam  intensities.  As  a  result,  achievable  beam  intensi¬ 
ties  were  only  marginally  satisfactory  for  scattering  measurements. 

Source  pressure  and  temperature  were  often  poorly  controlled. 

A  new  source  system  was  designed  and  built  to  overcome  these  diffi¬ 
culties.  The  major  design  criteria  were  to  be  able  to  achieve  a  stream 
of  high  intensity  consistent  with  the  available  pumping  to  direct  the 
beam  so  that  most  of  it  strikes  the  target,  to  have  independent  control 
over  source  temperature  and  pressure,  and  to  make  the  source  readily 
adaptable  to  use  with  non-condensable  gases  in  addition  to  water. 

A  nozzle/single  collimator  system  was  selected  as  simplest  to  build 
and  adequate  for  the  task.  There  appeared  to  be  no  substantial  advan¬ 
tage  to  carefully  designed  skimmers  as  long  as  the  collimator  was  located 
well  downstream  of  the  Mach  disk  [45].  A  single  stage  collimation 
system  was  adequate,  provided  the  pumping  speed  in  the  collimating 
chamber  was  high  enough  to  render  scattering  in  the  collimating  chamber 
unimportant.  In  the  present  system  this  condition  was  achieved  for 
water,  but  probably  not  for  argon. 

The  properties  of  flows  through  nozzles  are  not  yet  well  understood 
in  any  detail,  especially  if  condensation  occurs.  The  only  "nozzle"  for 
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which  detailed  calculations  can  be  made  is  the  effusive  source.  De¬ 
fine  d  as  a  characteri stic  linear  nozzle  dimension  perpendicular  to  the 
beam  direction  and  £  as  a  characteristic  dimension  along  the  beam.  An 
effusive  beam  occurs  for  £<<d<<x,  where  x  is  tne  mean  free  path  for 
molecules  in  the  source.  Rough  calculations  show  that  for  practical 
scattering  systems,  it  is  not  possible  to  achieve  a  sufficient  intensity 
with  such  a  source.  With  a  nozzle  source,  on  the  other  hand,  one  can 
achieve  arbitrarily  high  intensity  by  increasing  the  source  pressure. 

In  addition  the  flow  tends  to  be  more  directional,  and,  if  high  pressure 
ratios  can  be  maintained  across  the  nozzle,  the  width  of  the  velocity 
distribution  of  the  beam  becomes  very  narrow. 

Several  shapes  are  possible  for  nozzles:  cylindrical  holes,  con¬ 
verging  ("sonic"),  and  converging-diverging  ("supersonic")  (either 
straight-edge  or  curved).  The  effect  of  some  of  these  configurations 
has  been  studied  in  some  detail  by  Obert  and  Hagena  [38,  39].  Cluster¬ 
ing  and  beam  intensity  were  enhanced  in  "supersonic"  nozzles.  Among 
supersonic  nozzles,  larger  throat  diameters  promoted  clustering  at  the 
expense  of  centerline  beam  intensity.  Large  £/d  promoted  directional¬ 
ity  (hence  centerline  intensity)  and  clustering  (by  extending  the  high 
pressure  supersaturated  region  in  the  flow).  Narrow  divergence  angles 
also  promoted  directionality,  and  a  growing  boundary  layer  made  the 
effective  divergence  angle  smaller  than  the  actual  physical  angle. 

At  the  time  we  designed  the  source  nozzle,  we  were  not  interested 
in  or  much  concerned  about  clusters.  We  chose  a  cylindrical  nozzle  with 
the  smallest  easily  available  hole  size  and  H d  a  20.  This  nozzle 
was  expected  to  give  high  intensity  and  good  directionality.  Boundary 
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layer  effects  probably  give  a  flow  pattern  of  the  type  shewn  -'n  rigure  3. 

This  is  not  an  ideal  supersonic  flow  profile,  but  experiments  indicated 
that  clustering  of  water  molecules  was  very  probable  and  actually  diffi- 
cult  to  avoid.  It  is  not  clear  where  in  the  nozzle  sufficient  super¬ 
saturation  occurred  for  condensation  to  begin,  and  we  nave  not  studied 
the  source  in  detail.  Qur  primary  interest  was  in  scattering,  and  for 
that  purpose,  the  detailed  behavior  of  the  flow  through  the  nozzle  was 
of  secondary  interest.  As  long  as  pumping  speed  limitations  were  not 
reached,  the  beam  intensity  and  clustering  could  be  controlled  over  a 
wide  range  by  adjusting  source  pressure  and  nozzle  temperature.  The 
only  unfortunate  consequence  of  our  choice  of  nozzle  shape  was  that  there 
was  not  a  body  of  data  on  performance  to  draw  on. 

The  nozzle  diameter  was  chosen  as  .1  mm,  which  was  the  smallest 
readily  available  hypodermic  tubing.  We  considered  drilling  methods, 
but  the  technology  for  drilling  holes  at  reasonable  cost  less  than  about 
.25  mm  was  mostly  limited  to  very  shallow  holes.  Electroplating  on 
mandrels  and  etching  out  of  the  mandrel  is  possible  [38]  but  more  diffi¬ 
cult. 

The  source  as  constructed  is  shown  in  Figure  4.  The  nozzle  open¬ 
ing  was  .1  mm  in  diameter  by  2  rtm  long.  It  was  stainless  steel  hypo¬ 
dermic  tubing  brazed  into  solid  copper.  A  thermocouple  (Cu/CuNi)  was 
mounted  in  the  copper  to  monitor  source  temperature.  A  50  watt  oven 
surrounded  the  copper.  This  oven  consisted  of  a  coil  of  nichrome  wire 
mounted  in  grooves  cut  on  the  inside  of  a  cylinder  of  lava.  The  noz¬ 
zle  was  connected  to  a  boiler  outside  the  vacuum  system  via  a  1.27  cm 
0.0.  stainless  steel  tube.  This  tube  was  installed  so  that  the  path 

1 


12 


from  the  nozzle  to  the  boiler  was  unobstructed  and  downnill.  There  was 
a  radiation  shield  around  the  tube  and  the  source  oven.  The  assembly 
functioned  as  a  heat  pipe.  With  no  added  heat  at  the  source  oven,  the 
source  temperature  rose  and  stabilized  a  few  degrees  below  the  boiler 
temperature. 

The  boiler  was  a  7.5  cm  stainless  steel  T  holding  approximately 
1.5  1  of  water.  Heat  was  supplied  by  a  1000  watt  cartridge  inmersion 
heater  and  temperature  was  monitored  by  a  thermocouple  (Cu/CuNi).  Tem¬ 
perature  controllers  on  the  boiler  and  source  oven  controlled  the  source 
pressure  and  temperature.  These  were  specially  built  to  hold  temperatures 
to  within  ±1/2°  C.  An  overpressure  -elief  valve  was  set  at  about  1.3  MPa. 
A  fill  tube  was  provided.  To  add  water  to  the  boiler,  the  temperature 
was  reduced  to  below  100°  C,  a  valve  was  opened,  and  water  was  driven 
in  by  atmospheric  pressure. 

The  inside  of  the  source  system  was  almost  entirely  stainless  steel. 
Exceptions  were  the  copper  at  the  nozzle,  some  silver  solder,  and  viton 
0-rings  at  the  flanges  of  the  boiler.  The  viton,  in  fact,  caused  some 
problems;  some  fragments  which  came  loose  were  transported  through  the 
source  system  and  ended  up  plugging  the  nozzle.  It  was  possible  to  clear 
the  nozzle,  but  a  future  design  should  include  copper  or  teflon  seals, 
or,  perhaps,  a  screen  between  the  boiler  and  nozzle. 

The  collimating  chamber  was  a  rectangular  region  of  welded  stainless 
steel  about  20  cm  square  extending  about  50  cm  into  the  main  chamber. 

The  region  was  pumped  by  a  45  cm  oil  diffusion  pump  connected  to  the 
collimating  region  by  a  large  T.  In  addition  there  was  a  complete 
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inner  wall  of  liquid  nitrogen  cooled  copoer  wnich  accomplisbec  most  of 
the  pumping  for  water.  A  large  access  port  was  provided  in  :ne  side 
of  the  collimating  chamber  and  a  mounting  platform  was  Put  in  the  cottcm. 
The  front  of  the  chamber  had  the  top  and  side  corners  cut  Pack  at  an 
angle  of  about  30°  to  provide  additional  manoeuver  •'g  room  for  tne 
chopper.  Several  different  collimating  openings  could  be  installed. 

These  were  circular  openings  in  76  micron  stainless  steel  *oil\  cor 
all  the  scattering  experiments,  a  .76  cm  diameter  opening  was  selected 
for  a  nozzle  located  12.4  cm  away. 

A  solenoid-controlled  beam  flag  was  installed  in  front  of  the 
collimator  to  make  it  possible  to  turn  the  beam  on  'd  off  quickly. 

2.3.  TARGETS 

2.3.1.  TARGET  MOUNT 

The  target  mount  is  shown  in  Figure  5.  Originally  it  was  designed 
to  fit  into  a  glass  enclosure  so  that  ice  crystals  could  be  grown  with 
the  assembly  in  place  in  the  main  vacuum  system.  This  feature  was  not 
used  in  the  present  work.  Motor  drives  allowed  the  target  to  be  rotated 
and  moved  up  and  down  (out  of  the  path  of  the  beam).  Calibrated  potenti¬ 
ometers  were  used  to  keep  track  of  target  position.  Temperature  control 
was  provided  by  resistive  heating  and  liquid  nitrogen  cooling.  An 
early  attempt  to  use  a  thermoelectric  device  proved  unsatisfactory  for 
the  temperature  range  of  interest.  The  final  design  described  here 
is  the  third  or  fourth  generation.  It  worked  well  over  a  temperature 
range  of  about  200°  K  <  T < 700°K.  Temperature  control  to  about  ±  1/2°  K. 
was  provided  by  feedback  control  to  the  heater.  A  thermocouple  (Cu/CuNi) 
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clamped  to  the  target  surface  provided  the  temperature  measurement,  ana 
it  was  this  temperature  which  was  controlled. 

The  assembly  was  very  compact.  Starting  from  the  copDer  mounting 
surface,  there  was  a  sheet  of  mica,  a  coil  of  nichrome  wire  (the  heater), 
a  sheet  of  mica,  and  another  sheet  of  copper.  Four  stainless  steel  screws 
with  alumina  spacers  held  this  assembly  to  a  liquid  nitrogen  cooled  cop- 
oer  block.  This  liquid  was  supplied  via  flexible  metal  tubing  from  a 
reservoir.  The  two  lines  both  connected  to  the  bottom  of  the  reservoir, 
and  the  liquid  self-circulated.  (No  pump  is  necessary;  there  are  two 
stable  flow  patterns:  the  1 iquid  can  flow  in  either  direction  but  cannot 

stay  still;  warm  fluid  rising  on  one  side  pulls  more  cold  fluid  down  on 

• 

the  other. ) 

The  target  itself  was  clamped  to  the  mount,  and  a  thermocouple 

clamped  to  the  face  of  the  target  measured  the  temperature  which  was  read 

and  controlled  as  the  target  temperature.  A  thin  film  of  a  heat  sink 

compound  (zinc  oxide  loaded  silicone  grease)  provided  thermal  contact 

for  all  but  the  ice  surfaces,  which  were  simply  frozen  on  (in  practice, 

% 

"glued"  with  water). 

2.3.2.  TARGETS 

Several  targets  were  used.  Their  preparation  and  characteristics 
are  described  separately. 

Aluminum:  a  standard  alloy  (2Q24-T4)  whose  surface  was  rough-polished 
by  a  dry  polishing  technique  developed  for  covellite  [46,  47]  was  blown 
clean  with  compressed  air  and  installed.  This  is  a  rough  dirty  ("engi¬ 
neering")  surface  from  the  point  of  view  of  a  molecular  beam. 
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Silicon:  a  (111)  cleaved  surface  was  used.  The  sample  had  been 
stored  in  air  for  about  ten  years.  It  had  been  used  by  Lewis  [48],  who 
had  obtained  diffraction  of  a  He  beam  from  the  surface.  The  surface 
should  be  quite  smooth,  but  may  well  be  oxidized. 

Platinum:  a  polycrystall ine  block  (whose  detailed  history  is  un¬ 
known)  was  cleaned  in  trichloroethylene,  acetone,  and  distilled  water. 

Glass:  a  pyrex  cover  slip  (for  microscope  slides)  was  placed  over 
an  aluminum  block  to  provide  the  ..tandard  target  thickness.  The  surface 
was  cleaned  by  the  process  used  for  platimun.  Thermal  conducting  grease 
was  put  between  glass  and  aluminum  and  between  aluminum  and  target  mount. 

Covellite:  no  practical  laboratory  methods  for  preparing  large 
single  crystals  of  CuS  are  available:  high  pressure  and  long  time  seem 
to  be  necessary.  Hence  it  was  found  necessary  to  use  the  natural  mineral 
as  it  comes  from  a  copper  mine  in  Butte,  Montana.  This  material  is 
known  to  provide  a  favorable  nucleating  surface  for  ice  [2,  46,  47], 
because  it  has  a  very  similar  lattice  structure  and  size.  A  piece  was 
selected  which  was  reasonably  strain  free  and  free  of  iron  pyrite  in¬ 
clusions.  A  surface  was  dry-polished  parallel  to  the  cleavage  planes 
(readily  visible).  The  dry  polishing  technique  [46,  47]  was  developed 
because  it  appeared  that  water  changes  the  surface  in  a  way  which  af¬ 
fected  the  growth  of  macroscopic  ice  crystals  from  the  vapor.  Auger 
spectroscopy  indicated  that  the  surface  was  indeed  copper  and  sulphur 
with  no  significant  impurities. 

Ice:  single  crystal  ice  was  obtained  from  the  Athabasca  glacier 
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in  Alaska.  Thermal  shock  induced  by  plunging  in  liquid  nitrogen  caused 
cracks  which  indicated  crystal  orientation.  A  sample  could  then  be 
cleaved  or  cut.  It  is  not  easy  to  characterize  the  surface  used,  and 
samples  were  not  saved.  Cleaving  ice  is  difficult.  3est  results  were 
achieved  for  temperature  near  the  melting  point,  but  even  then  results 
were  inconsistent.  No  attempts  were  made  with  any  of  the  samples  to  ob¬ 
tain  atomically  clean  surfaces. 

At  10'®  Pa,  a  monolayer  forms  on  a  surface  in  a  few  seconds.  There 
was  therefore  no  point  in  using  conventional  high  vacuum  surface  clean¬ 
ing  techniques  (such  as  argon  ion  bombardment).  One  has  to  be  aware, 
though,  in  interpreting  the  scattering  data,  that  the  surfaces  were  not 
initially  clean.  Adsorbed  water,  at  least,  was  initially  driven  off  by 
heating  targets  to  400°  K  for  at  least  15  minutes.  (330°  K  for  covellite 
which  undergoes  a  transition  at  376°  K  and  decomposes  at  around  490°  K  ) 
For  scattering  which  seemed  to  depend  on  an  accumulated  surface  layer 
of  some  sort,  an  additional  preparation  step  was  taken:  the  target  was 
cooled  to  240°  K  and  exposed  to  a  normally  incident  beam  for  at  least 
15  minutes. 


*  Here  again  high  pressure  and  long  times  are  normally  necessary  for  the 
formation  of  large  single  crystals.  The  best  laboratory  results  are  a- 
chieved  by  epitaxial  growth  from  the  vapor  [47],  but  the  practical  limit 
seems  to  be  crystals  a  few  millimeters  in  diameter;  which  is  too  small 
for  scattering  experiments.  A  zone  refining  technique  [49]  has  also  a- 
chieved  large  single  crystals  and  may  be  the  best  approach, 
t  Studies  in  our  lab  [50]  following  the  work  of  Knight  and  Knight  [51] 
suggest  that  an  alternative  to  thermal  cracking  and  cleaving  or  cutting  is 
the  use  of  ''negative"  ice  crystals.  If  a  hypodermic  tube  is  inserted 
into  a  single  crystal  sample  of  ice,  and  a  vacuum  is  maintained  in  the  tube, 
crystal  shaped  cavities  form.  This,  at  least,  gives  a  good  indication  of 
crystal  orientation.  Sizes  achieved  thus  far  are  still  in  the  millimeter 
range,  but  larger  smooth  crystal  faces  may  be  possible  if  one  can  avoid 
the  tendency  toward  a  stepped  structure.  For  orientation  determination, 
one  could  probably  also  make  do  with  the  crystal-shaped  etch  pits  formed 
when  a  single  crystal  ice  sample  is  put  in  a  vacuum  for  a  few  hours  at 
around  -250°  K  [50]. 
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2.4. 


DETECTION  EQUIPMENT 


2.4.1.  VELOCITY  SPECTROMETER 

The  primary  measuring  system  for  the  present  work  was  a  time-of-fl ight 
velocity  spectrometer.  A  brief  description  of  the  basic  features  of  the 
system  is  given  in  Appendix  A  which  is  a  copy  of  Reference  [53].  The 
components  of  the  system  were  essentially  as  described  in  Reference  [52]. 

A  photoetched  pseudorandom  chopper  disk  (see  Figure  2  of  Appendix  A)  re¬ 
placed  the  machined  one  described  in  Reference  [52].  The  chopper  and 
ionizer  were  mounted  on  independent  arms  which  could  be  revolved  around 
the  target  in  the  principal  plane  defined  by  the  target  normal  and 
incident  beam.  The  position  of  these  arms  could  be  controlled  and  meas¬ 
ured  from  outside  the  vacuum  system.  The  independent  motion  allowed  the 
chopper  to  be  positioned  to  chop  either  the  incident  beam  or  reflected 
stream.  To  obtain  precise  alignment  of  the  two  arms,  a  bare  wire-wound 
trim  potentiometer  was  mounted  on  one  arm  with  a  wiper  on  the  other. 
Alignment  to  better  than  .1°  was  easily  achieved. 

Distances  and  dimensions  are  indicated  in  Figure  6.  The  history 
of  the  various  baffles  is  given  in  detail  in  Section  3.1.2. 

The  ionizer  (the  same  one  as  was  used  in  Reference  [52])  was  mounted 
in  an  open  5.4  cm.  copper  T,  which  was  cooled  with  liquid  nitrogen 
to  keep  the  background  gas  density  at  a  minimum  in  the  ionizing  region. 

The  electron  multiplier  was  mounted  in  the  arm  of  the  T  (vertically), 
so  that  unionized  beam  particles  passed  unobstructed  out  of  the  T.  Early 
experiments  were  carried  out  using  an  18  dynode  beryl  1 ium-copper  Venetian 
blind  type  multiplier  (made  by  EMI).  This  was  replaced  with  a  continuous- 
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dynode  semiconductor  multi  pi ier  ("channeltron')(Gali1eo  Electron  Optics 
model  4716),  which  thus  far  apoears  to  have  a  much  longer  usable  lifetime. 

Signal  processing  after  the  electron  multiplier  is  described  in 
Section  2.5. 

2.4.2.  STAGNATION  GAUGE 

For  some  gas  flow  measurements,  a  glass-enclosed  ionization  gauge 
with  a  slit  opening  was  mounted  on  one  of  the  rotating  arms.  The  open- 
ing  was  .12  cm  ,  and  the  volume  about  151  cmJ5  which  gives  a  gauge 
time  constant  of  about  .08  sec.  This  does  not  agree  with  experimental - 
ly  measured  time  constants: for  pressures  of  around  10  Pa  the  measured 
time  constant  was  about  2  sec,  and  for  lower  pressures  (of  water,  in 
particular),  the  response  time  is  much  slower:  about  20  sec  at  10”4  Pa, 
probably  due  to  adsorption  on  the  glass  wall. 

2.4.3.  MASS  SPECTROMETER 

A  late  but  important  addition  to  the  detection  apparatus  was  a 
UTI  model  100C  quadrupole  mass  spectrometer  equipped  with  a  channel tron 
electron  multiplier.  Due  to  the  difficulty  of  mounting  the  mass  spec¬ 
trometer  probe  on  one  of  the  moving  arms,  it  was  mounted  in  two  fixed 
positions,  one  facing  the  beam  directly,  the  other  forming  an  angle  of 
144°  from  the  source  to  the  target  to  the  mass  spectrometer  ionizer 
(Figure  1). 

For  studies  of  the  beam,  the  range  of  the  instrument  was  extended 
from  300  AMU  to  about  700  AMU  by  reducing  the  rf  frequency  of  the 
quadrupole  from  1.70  MHz  to  1.11  MHz.  Some  loss  of  attainable  resolution 
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occurred  when  this  was  done.  Further  consequences  of  extending  the 
mass  range  are  discussed  in  more  detail  in  Section  3.1.1. 

For  some  experiments,  time-of-fl ight  measurements  were  made  replac¬ 
ing  the  original  ioni2er  with  the  mass  spectrometer’s  ionizer/quadrupole/ 
electron  multiplier  assembly.  This  allowed  simultaneous  mass  and  velocity 
resolution.  For  these  experiments,  the  output  of  the  electron  multiplier 
in  the  mass  spectrometer  was  fed  directly  through  the  usual  time-of-fl ight 
electronics,  while  the  mass  spectrometer  was  locked  onto  a  particular 
mass  peak. 

2.5.  TIME-OF-FIIGHT  SIGNAL  PROCESSING 

The  signal  path  is  diagrammed  in  Figure  7.  The  electron  multiplier 
served  as  a  low  noise  high  gain  preamplifier.  The  signal  then  passed 
through  a  current  preamplifier  (gain  10^  V/A)  and  an  additional  amplifier 
(gain  50),  if  necessary,  to  achieve  a  signal  of  a  few  volts.  The  ampli¬ 
fied  signal  was  generally  very  noisy,  and  some  signal  averaging  was 
necessary.  All  the  data  reported  here  were  processed  through  Princeton 
Applied  Research  model  TDH9  Waveform  Eductor  which  is  a  100  channel 
analogue  signal  averager.  (The  Eductor  memory  consisted  of  a  bank  of 
precision  capacitors.)  Comparison  with  a  demonstration  model  of  a  new 
digital  signal  averager  made  by  Nicolet  Instruments  (model  1170)  showed 
the  Waveform  Eductor  to  be  a  very  effective  device  for  processing  the 
signals  we  generated.  Since  it  used  the  entire  signal  (no  sample-and- 
hold),  it  was  actually  somewhat  more  efficient  than  the  digital  instru¬ 
ment.  The  digital  device  was  better  for  very  noisy  signals,  of  course, 
because  it  could  continue  to  average  indefinitely  or  at  least  until 
available  memory  is  full.  Also  the  A/D  conversion  was  sufficiently 
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fast  that  one  could  achieve  about  a  ten  fold  improvement  in  time  resolu¬ 
tion  over  the  Waveform  Eductor.  This  had  interesting  consequences  for 
the  signal  from  a  103  channel  pseudorandom  chopper.  The  signal,  after 
sufficient  averaging,  still  looked  slightly  noisy,  because  the  effect 
of  the  shutter  function  (defined  by  an  individual  chopper  slot  passing 
the  chopper  collimator)  was  resolved:  there  were  103  little  bumps  super¬ 
imposed  on  the  signal! 

The  trigger  signal  for  the  Waveform  Eductor  was  supplied  by  a 
photodiode-phototransistor  which  sensed  the  passage  of  a  long  slot  which 
occurred  once  for  each  repetition  of  the  chopping  sequence.  A  dead 
space  of  3  channels  was  set  so  that  the  Waveform  Eductor  channel  width 
matched  the  chopper  channel  width.  This  dead  space  was  usually  set 
to  span  the  time  of  the  phototransistor  signal  pulse  which  was  picked 
up  slightly  on  the  signal  line.  The  missing  channels  were  interpolated. 
For  most  signals,  interpolation  did  not  introduce  noticeable  error,  since 
the  deconvolution  process  for  pseudorandom  chopping  used  all  channels 
to  reconstruct  each  channel  of  the  deconvoluted  signal. 

The  averaged  signal  was  displayed  on  an  oscilloscope.  When  a 
signal  of  interest  was  ready,  it  was  read  through  a  buffer  amplifier 
(to  protect  the  A/0  converter  from  voltage  spikes)  and  A/D  converter 
into  an  on-line  PDP-11 /03  computer.  This  A/0  converter  was  not  fast 
enough  to  handle  the  signal  directly,  which  is  why  the  Waveform  Eductor 
was  necessary.  A  typical  period  (time  for  103  channels)  was  1-2  msec. 

For  read-in  to  the  computer,  a  period  of  about  .1  sec.  was  used.  Once 
the  signal  was  in  the  computer,  a  great  variety  of  manipulations  were 
possible.  These  are  described  in  Section  2.6. 
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2.6  SOFTWARE 


A  series  of  routines  were  developed  to  process  the  data.  They  are 
given  in  Appendix  3  along  with  user  instructions.  All  were  written  in 
FORTRAN  and  run  under  the  RT -IT  version  3  operating  system.  System  hard¬ 
ware  included  dual  floppy  disks,  a  programmable  clock,  a  VT55  graphic 
display  terminal  with  electrostatic  hard  copy  unit,  and  a  Hewlett-Packard 
model  7225A  digital  plotter.  Special  subroutines  ware  available  from 
the  manufacturer  or  were  written  for  these  devices. 

Subroutines  were  written  to  be  self-contained  and  accessible  in 
any  order.  Basic  data  handling  routines  included  routines  to  read  in 
new  data  and  store  it  on  disk,  read  in  old  data  off  disk,  and  to  con¬ 
solidate  disk  files  at  the  end  of  an  experimental  run  (each  trace  was 
initially  stored  in  a  separate  file  to  provide  crash-resistance  for  the 
data).  The  data  could  be  plotted  directly  as  read  in.  All  plotting 
could  be  done  either  on  the  graphic  display  terminal  or  the  digital 
plotter.  Hard  copies  of  the  display  on  the  terminal  were  the  most  com¬ 
mon  means  of  collecting  paper  copies  of  graphs  since  it  was  faster  than 
using  the  plotter.  Good  plotter  copies  could  always  be  generated  as 
needed,  since  all  data  were  on  file. 

The  data  processing  routines  are  as  follows:  first,  of  course,  a 
deconvolution  routine  (52]  recovered  a  time-of-fl ight  spectrum  from  the 
pseudorandom  chopped  signal.  As  long  as  one  channel  width  time  resolu¬ 
tion  was  satisfactory,  no  further  corrections  were  necessary.  See 
Reference  [52]  and  Appendix  A  for  a  discussion  of  shutter  function 
effects  and  the  effects  of  a  finite  length  ionizer.  A  smoothing 
routine  was  available  using  a  low-pass  digital  filtering  scheme  [59].  In 
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general,  this  routine  was  not  found  to  be  useful.  The  noise  accompanying 
the  signal  was  due  primarily  to  statistical  fluctuation  in  the  back¬ 
ground  gas  density.  As  such,  it  was  white  noise.  A  low  pass  filter  cut 
out  the  high  frequency  components,  but  left  low  frequency  ones.  The 
human  eye/brain  handled  the  unfiltered  version  better  since  the  low 
frequency  noise  might  look  like  signal  when  it  alone  was  present. 

Various  theoretical  olots  were  relevant  to  the  measurements  and 
appropriate  fitting  routines  were  devised  for  comparing  the  data  to 
them.  Most  of  these  routines  fitted  a  Maxwellian  or  translating  Max¬ 
wellian  or  the  sum  of  two  such  velocity  distributions  to  the  data.  To 
avoid  the  hazards  of  non-linear  least-squares  fitting,  only  amplitudes 
were  used  as  fitting  parameters.  Other  parameters  were  taken  as  fixed, 
and  if  they  were  not  independently  measurable,  they  were  determined 
manually  by  a  guess  and  try  method  (manual  non-linear  least-squares 
fitting).  Such  guessed  parameters  included  Mach  number  and  temperature 
for  translating  Maxwellians;*  for  static  Maxwellians,  all  parameters 
were  known.  Fits  could  be  made  to  all  or  part  of  an  experimental  curve 
as  appropriate.  Routines  were  also  available  to  subtract  out  parts  of 
a  signal.  For  example,  it  was  often  desirable  to  subtract  out  a  Max¬ 
wellian  component  to  be  better  able  to  study  the  remaining  features. 

In  practice,  the  most  useful  information  could  be  extracted  by 
these  fitting  routines,  or  by  an  even  simpler  routine  which  just  identi¬ 
fies  peak  positions  and  amplitudes.  This  is  in  accord  with  common  prac¬ 
tice.  See  Appendix  C  for  a  more  detailed  discussion. 

*  For  "good"  data  (signal/noise  > 50) ,  experience  with  fitting  suggested 
that  these  parameters  could  be  determined  with  some  confidence  (e.g.±5S 
for  temperature)  by  this  technique.  Note,  however,  the  caution  at  the  end 
of  Appendix  C. 
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SECTION  III 


DATA 

A  lot  of  data  were  collected  for  this  study.  There  were,  for  exam¬ 
ple,  some  2000  time-of-fl ight  spectra  taken  in  the  course  of  some  100 
experimental  runs.  It  is,  therefore,  not  practical  to  present  all  the 
raw  data.  They  are  stored  on  disk  and  in  a  file  at  the  lab.  (Fluid 
Dynamics  Research  Lab,  Rm.  37-4A2,  MIT)  Instead,  the  information  has 
been  consolidated  and  combined  as  much  as  possible.  Hopefully,  this 
has  made  the  material  more  digestible,  and  no  important  data  have  been 
omitted. 

The  data  are  presented  and  discussed  in  four  main  divisions.  First 
there  is  a  discussion  of  instrumental  calibration  and  of  spurious  sig¬ 
nals  which  were  eliminated  and  are  not  present  in  any  of  the  remaining 
data.  Then  there  is  a  discussion  of  the  properties  of  the  beam.  A 
few  results  of  some  studies  of  evaporation  follow,  and  finally  the 
scattering  data  are  presented. 

3.1.  CALIBRATIONS  AND  ELIMINATION  OF  ARTIFACTS 

3.1.1.  MASS  SPECTROMETER 

The  mass  calibration  was  factory  set  for  the  1-300  AMU  nominal 
range  of  the  instrument.  Checks  against  known  gases  (He,  Ng,  H2O,  O2, 
COg,  Ar,  Kr)  verified  this  calibration,  at  least  for  the  low  AMU  range. 
Above  84  AMU  (Kr),  and  for  the  instrument  as  modified  for  higher  mass 
range,  the  clusters  in  the  water  beam  provided  a  very  convenient  cali¬ 
bration  series.  There  are  peaks  in  the  mass  spectrum  of  a  clustered 
water  beam  at  18n+l  AMU  intervals  corresponding  to  protonated  water 
clusters. 
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Mass  peak  amplitude  calibration  is  a  more  difficult  matter.  For 
the  most  part  the  problem  was  ignored,  and  quantitative  measurements 
were  not  made.  To  make  a  satisfactory  calibration  requires  known 
amounts  of  the  species  to  be  measured  and  approximately  the  same  vacuum 
and  flow  conditions.  Short  of  this,  several  corrections  can  be  made, 
which  should  allow  relative  peak  amplitude  comparison  at  least  to  within 
about  a  factor  of  two. 

First,  ionization  efficiency  is  a  function  of  AMU.  Data  provided 
by  the  manufacturer  [55]  taken  from  work  by  Flaim  and  Ownby  [56]  suggest 
as  a  good  approximation: 

<2-'> 

where  £  is  the  ionization  efficiency  relative  to  N2,  and  Z  is  the  sum 
of  the  atomic  numbers  (number  of  electrons)  of  the  atoms  in  the  ion 
in  question. 

Second,  the  gain  of  the  electron  multiplier  depends  on  species  in 
an  irregular  manner  depending  on  mass,  charge,  reactivity,  and  accel¬ 
erating  voltage.  Roughly,  one  finds  that  there  is  an  m"1^  dependence[55] 
where  m  is  the  mass  of  the  ion  or  ion  fragment. 

Finally,  transmission  through  the  quadrupole  is  also  a  function  of 
m.  For  the  instrument,  as  delivered,  the  manufacturer  reports  [55] 
that  transmission  is  near  lOOt  for  less  than  40  AMU  and  drops  off  ap¬ 
proximately  one  decade  for  every  150  AMU  thereafter.  This  is  presumably 
correct  for  some  standard  setting  of  the  resolution  control .  There  is 
a  tradeoff  between  quadrupole  transmission  and  resolution.  No  estimate 
of  the  magnitude  of  the  effect  of  this  tradeoff  has  been  made. 
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As  the  r f  frequency  is  lowered  to  extend  the  mass  range,  the  transmis¬ 
sion  goes  up  roughly  proportionally  to  the  mass  range.  (This  last 
observation  is  based  on  some  qualitative  observations  and  may  be  inac¬ 
curate.  The  adjustments  in  the  resolution  setting  which  go  along  with 
extending  the  mass  range  may  be  more  important.)  As  yet,  we  have  not 
been  able  to  apply  corrections  satisfactorily  to  obtain  relative  ampli¬ 
tudes  over  a  large  mass  range. 

3.1.2.  TIME-OF-FLIGHT  APPARATUS  WITHOUT  MASS  FILTERING 

Many  problems  were  encountered  before  satisfactory  operation  of  the 
time-of-flight  apparatus  was  achieved.  Most  of  these  are  associated 
with  the  time  resolution  and  sensitivity  achieved  with  the  present 
equipment,  and  were  not  noticeable  with  most  of  our  previous  equipment. 
Problems  were  found  associated  with  each  of  the  components  of  the  system: 
chopper,  ionizer,  and  electron  multiplier.  We  consider  each  of  these 
problems  and  their  solutions  separately. 

When  the  photoetched  pseudorandom  chopper  disk  was  first  installed, 
we  observed  consistent  spurious  "bumps"  on  the  time-of-flight  spectrum. 
These  turned  out  to  be  due  to  a  drafting  error  by  the  company  which  made 
the  disks.  One  line  representing  the  edge  of  a  slot  was  drawn  1/2  chan¬ 
nel  off  from  where  it  should  have  been.  Simulating  the  effect  of  this 
error  on  a  computer  indeed  produced  bumps  of  the  kind  we  saw.  The  error 
was  corrected. 

A  more  elusive  problem  was  a  spurious  signal  which  appeared  to  be 
due  to  particles  which  passed  through  some  portion  of  the  chopper  and 
struck  at  least  one  other  surface  before  entering  the  ionizer.  We  also 
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suspected  that  beam  particles  might  strike  the  edge  of  a  chopper  slot 
and  continue  on  to  produce  a  conerent  signal.  The  evidence  *or  this 
comes  from  observation  of  variations  in  the  spurious  signals  as  a  func¬ 
tion  of  chopper  rotation  direction.  Note  that  the  pseudorandom  chopper 
provides  a  substantially  larger  total  area  of  slot  edge  per  cycle  than 
does  a  single  slot  chopper.  These  spurious  signals  manifested  them¬ 
selves  as  an  assortment  of  bizarre  waveforms  (i.e.  they  could  not  readi¬ 
ly  be  explained  as  time-of-fl ight  spectra)  which  could  be  generated 
even  with  gross  misalignments  of  the  components  of  the  apparatus  (for 
example,  with  the  chopper  at  90°  from  the  target-ionizer  line).  A 
set  of  properly  placed  baffles  solved  these  problems.  Disks  were  placed 
shielding  both  sides  of  the  chopper  disk  except  where  the  beam  was  to 
pass  through.  Baffles  attached  to  the  ionizer  restricted  the  possible 
multi  bounce  paths  from  the  chopper  to  the  ionizer.  There  is,  of  course, 
the  potential  hazard  that  all  this  extra  baffling  might  introduce  prob¬ 
lems  of  its  own.  In  particular,  the  baffles  next  to  the  chopper  might, 
in  effect,  momentarily  trap  blocked  beam  particles  which  could  then 
pass  through  the  chopper  when  next  it  opened.  This  would  produce  an 
unpredictable  coherent  noise  source  which  would  not  be  averaged  out. 

We  can  only  say  that  there  is  no  evidence  in  any  data  for  an  artifact 
of  this  sort.  With  all  the  baffling  in  place  signals  were  obtained 
when  and  only  when  they  were  expected,  and  a  good  calibration  match 
could  be  obtained  with  a  Maxwellian  stream  (see  below). 

The  ionizer's  job  is  to  form  ions  from  the  neutral  beam  and  pass 
them  on  to  the  electron  multiplier  in  a  time  interval  negligible  com¬ 
pared  to  the  flight  time  between  chopper  and  ionizer.  Not  all  ionizer 
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designs  successfully  perform  this  task.  A  box  ionizer  (similar  to  many 
widely  used  ionizers)  which  we  used  initially  appeared  to  introduce 
delays  of  tens  of  microseconds.  When  that  ionizer  was  replaced  with  a 
tapered  cylindrical  ionizer  which  provided  a  draw  out  potential  of  a 
few  volts  throughout  the  ionizing  region,  the  delays  disappeared.  This 
problem  reappeared  when  mass-filtered  time-of-fl ight  spectrometry  was 
attempted  as  noted  in  Section  3.1.3. 

Vibration  of  the  ionizer  can  cause  an  undersired  signal  by  changing 
the  filament-grid  spacing.  Most  sources  of  vibration  are  not  at  the 
chopper  frequency,  so  that  as  long  as  they  do  not  produce  too  large 
a  signal,  they  average  out  nicely.  The  chopper  itself,  however,  was 
initially  a  source  of  vibration  at  the  chopper  frequency.  With  the 
chopper  mounted  on  (hung  from)  springs,  and  the  chopper  and  ionizer 
on  independent  arms,  this  problem  was  eliminated. 

The  major  problem  caused  by  the  electron  multiplier  was  due  to  the 
fact  that  its  gain  had  both  short  and  long-term  drifts  which  made  it 
difficult  to  obtain  precise  amplitude  information.  One  problem  did 
surface  briefly  when  the  channeltron  electron  multiplier  was  first  in¬ 
stalled.  The  time-of-fl ight  peaks  inexplicably  appeared  to  split  into 
double  peaks.  The  problem  was  apparently  again  one  of  non-negl igible 
ion  flight  time— in  this  case,  probably,  a  short  and  a  long  time  both 
being  possible.  A  grid  across  the  cone-shaped  opening  of  the  multiplier 
set  at  the  multiplier  high  voltage  (typically  -1500  V  to  -2000  V)  solved 
the  problem,  presumably  by  providing  a  stronger  field  to  attract  the 
ions  to  the  multiplier. 
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With  all  these  p'roolems  'roned  out  a  careful  calibration  of  the 
instrument  was  possible  (and  necessary  in  light  of  the  large  number 
of  things  which  can  evidently  go  wrong  with  a  time-of-fl ight  system). 
First,  of  course,  a  check  was  made  to  make  sure  that  signals  were  a- 
chieved  only  under  conditions  wnere  they  were  expected.  One  extraneous 
signal  source  was  never  fully  eliminated:  a  small  electrical  pickup 
of  the  phototrigger  pulse.  For  the  most  part,  the  dead  time  in  the 
Waveform  Eductor  sweep  could  be  set  to  span  this  pulse,  but  occasionally 
one  or  two  data  channels  might  have  to  be  discarded  and  interpolated. 

Second,  the  time  origin  had  to  be  located  precisely.  This  pro¬ 
cedure  has  been  previously  described  in  some  detail  [52,53]  but  as  the 
most  recent  description  was  never  published,  it  is  included  in  this 
document  as  Appendix  A.  Note  especially  Section  A4. 

Once  the  time  origin  was  located,  it  was  possible  to  make  a  direct 
calibration  comparison  using  a  beam  whose  velocity  distribution  is  known. 
The  only  such  distribution  available  is  the  Maxwellian.  A  reliable 
source  of  a  Maxwellian  stream  is  an  effusive  flow  at  high  Knudsen  num¬ 
ber  using  an  inert  gas.  Such  a  calibration  was  performed  using  He,  Ne, 

A r,  and  Kr.  Data  are  presented  in  References  [52,53].  The  reported 
time  origin  discrepancy  [52]  has  been  eliminated.  The  calibration 
has  been  generally  confirmed  by  some  of  the  experimental  data  itself 
which  include  studies  of  two  additional  sources  of  near  Maxwellian 
streams,  evaporation  from  a  surface,  and  backscattering  from  a  surface 
(which  in  the  present  experiments  consisted  only  of  trapped-reemitted 
particles).  While  these  are  not  satisfactory  for  primary  calibration, 
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they  do  provide  a  convenient  check  on  the  continued  oroDer  function¬ 
ing  of  the  apparatus  with  time. 

3.1.3.  TIME-OF-FLIGHT  APPARATUS  WITH  MASS  FILTERING 

The  mass-fil tered  version  of  the  time-of-flight  apparatus  (described 
in  Sections  2.4.1  and  2.4.3)  was  identical  to  the  unfiltered  version 
except  in  the  region  between  the  ionizer  and  electron  multiplier.  For 
the  quadrupole  mass  filter,  ions  are  accelerated  only  to  about  15  V.  The 
filtering  action  is  not  directly  dependent  on  this  voltage,  but  for 
proper  functioning,  the  ion  energy  must  be  kept  low  enough  so  that 
the  ions  are  in  the  finite  length  quadrupole  field  for  at  least  several 
rf  cycles.  In  principle,  the  ion  travel  time  should  still  be  negligi¬ 
ble:  at  15  V,  the  ions  travel  about  1000  times  faster  than  the  neutral 
beam  and  traverse  the  quadrupole  in  about  1  ys. 

Several  tests  were  performed  which  allow  us  to  assess  the  reli¬ 
ability  of  the  mass-resol ved  time-of-flight  data.  First,  the  system 
was  checked  to  see  if  any  ionizer  parameters  (electron  energy  and  cur¬ 
rent,  ion  energy,  and  focus  voltage)  affected  a  time-of-flight  spectrum 
in  any  way  other  than  the  expected  amplitude  changes.  Moderate  changes 
(up  to  about  50  %)  in  these  parameters  from  the  normal  settings  recom¬ 
mended  and  set  by  the  manufacturer  had  no  measurable  effect.  Larger 
changes  caused  waveform  distortions  and  overall  delays,  but  the  normal 
settings,  gave  reliable  reproducible  results. 

Next,  comparisons  were  made  between  time-of-flight  spectra  with  and 
without  mass-filtering.  An  argon  beam  provided  a  signal  which  was  al¬ 
most  entirely  at  40  AMU.  Since  the  flight  paths  for  the  two  versions 
were  different,  the  comparison  was  made  by  fitting  the  same  translating 
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Maxwellian  to  both  spectra.  A  satisfactory  match  was  achieved,  but  only 
by  introducing  a  delay  of  330  us  to  the  fitted  cur^  for  the  mass-fil¬ 
tered  spectra  (cf.  the  chopper  period  of  ?.96  ms  and  peak  position  at 
1.09  ms).  A  similar  calibration  using  water  evaoorating  from  a  surface 
(normal  emission  direction)  and  the  mass  spectrometer  positioned  to 
look  at  a  scattered  stream  (the  beam  passes  through  the  ionizer  in  a 
different  direction  in  this  case)  yielded  an  apparent  delay  of  90  us. 

The  reason  for  the  discrepancy  has  not  been  determined,  and  additional 
tests  need  to  be  made.  If  the  difference  is  due  entirely  to  the  change 
in  the  orientation  of  the  beam  relative  to  the  ionizer,  then  the  delay 
can  be  measured  and  the  appropriate  correc^'on  can  be  applied  to  all 
subsequent  data.  The  possibility  of  mass-dependence  has  not  yet  been 
ruled  out,  however.  Ideally,  of  course,  the  delay  should  be  eliminated. 
The  use  of  an  ionizer  of  the  type  we  built  for  our  time-of-flight 
apparatus  without  mass-filtering  might  eliminate  the  delay,  but  one 
might  still  expect  problems,  since  in  our  design,  there  is  an  accelerating 
potential  of  1.5-3  kV  drawing  the  ions  from  the  ionizer  to  the  electron 
multiplier,  whereas  in  the  mass  spectrometer  there  is  only  a  15  V 
potential  drawing  ions  into  the  quadrupole  field. 

The  position  we  have  taken  tentatively  in  the  interpretation  of 
mass-filtered  time-of-flight  data  is  that  for  a  given  experimental 
configuration,  there  is  a  constant  delay  independent  of  species  (mass) 
and  velocity.  The  magnitude  and  constancy  of  the  delay  is  still  un¬ 
certain  and  this  should  be  borne  in  mind  in  interpreting  the  data;  how¬ 
ever,  the  shape  of  the  distribution  should  be  reliable. 

In  the  category  of  artifacts,  there  were  some  rather  interesting 
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observations  made  with  mass-filtered  time-of-fl ight  soectrometry.  It 

was  oossible  to  obtain  quite  satisfactory  time-of-fl ight  spectra  for 

species  present  in  very  small  quantities,  since  the  statistical  noise 

is  only  that  associated  with  density  fluctuations  of  the  species  being 

observed.  It  was,  for  example,  quite  easy  to  obtain  a  time-of-fl ight 
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spectrum  A r  and  Ar  in  an  argon  beam.  Unfortunately,  we  could  also 

obtain  clean  signals  at  some  mass  numbers  which  could  not  be  associated 
with  anything  in  the  beam(see  Figure  8)  and  were  attributed  to  species 
present  in  the  background  gas.  The  time-of-fl ight  spectra  for  these 
background  species  were  different  from  those  of  beam  species.  In  fact, 
the  spectra  cannot  be  interpreted  as  time-of-fl ight  spectra  because 
they  dipped  below  zero  in  intensity.  The  dip  below  zero  is  real  (not 
an  artifact  of  overlapping  spectra);  extending  the  chopper  period 
substantially  showed  the  steady-state  background  level  to  be  above  the 
dip.  A  demonstration  of  this  phenomenon  was  made  using  a  mixed  Ar/He 
beam.  Normal  spectra  were  found  at  4  and  40  AMU,  while  anomalous  spectra 
were  found  at  18,  28,  32,  and  44  AMU  (see  Figure  9).  No  signal  was 
achieved  for  any  mass  not  present  in  the  mass  spectrum  of  the  beam  or 
background. 

It  is  not  possible  to  explain  these  signals  as  due  simply  to 
chopped  background  gas.  The  number  of  particles  which  just  happened  to 
be  passing  in  the  right  direction  to  go  through  the  chopper  to  the  ion¬ 
izer  was  too  small.  Experimentally,  this  was  readily  verified  by  turning 
off  or  blocking  the  beam;  in  which  case  no  time-of-fl ight  signal  could 
be  observed.  In  order  to  explain  the  signals,  one  must  postulate  a 
mechanism  involving  entrainment,  or  more  precisely,  single  collisions 
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between  beam  and  background  particles.  To  account  for  significant 
signal  intensity,  such  collisions  must  occur  after  the  chopper,  and, 
most  likely,  in  the  immediate  vicinity  of  the  ionizer.  As  long  as  only 
a  single  collision  is  involved,  it  is  reasonable  to  suopose  that  corre¬ 
lation  with  the  chopper  is  well  maintained.  Collisions  could  also 
knock  particles  out  of  the  ionizer  which  might  otherwise  be  ionized, 
and  this  could  explain  the  dip  below  the  background  level.  A  detailed 
modeling  of  this  process  has  not  been  attempted,  but  it  seems  to  be 
the  only  plausible  explanation  for  the  observed  features.  At  present 
we  know  of  no  way  of  extracting  useful  information  from  these  signals. 
Since  they  can  be  readily  identified  by  their  signature,  there  appears 
to  be  no  problem  identifying  spectra  which  are  associated  with  the  beam. 
It  is,  however,  interesting  to  note  that  the  signals  from  background 
species  must  also  contribute  to  the  signal  one  obtains  without  mass- 
filtering.  Fortunately,  it  is  unlikely  that  a  situation  could  be  devised 
where  these  spurious  signals  would  contribute  measurably  to  the  total 
signal.  At  least  for  the  mechanism  suggested  above,  there  must  be  many 
more  chopped  beam  particles  than  affected  background  particles  reaching 
the  ionizer. 

3.2.  SOURCE  STUDIES 

The  primary  focus  of  the  present  study  was  on  the  scattering  of  a 
water  beam  from  surfaces.  First,  however,  the  characteristics  of  the 
source,  and  more  particularly,  the  water  beam  had  to  be  understood.  The 
initial  study  of  the  source  was  with  argon. 
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3.2.1.  ARGON 


The  relatively  low  pumping  speed  for  argon  in  the  collimating  cham¬ 
ber  limited  source  pressures  to  values  below  25  kPa ;  at  -25  kPa,  the 

_2 

pressure  in  the  collimating  chamber  rose  to  about  10  Pa,  resulting  in 
significant  attenuation  of  the  beam  due  to  scattering.  The  velocity 
distribution  of  the  Ar  beams  could  be  fitted  well  to  a  translating 
Maxwellian,  when  the  beams  were  generated  at  source  pressures  below 
about  18  kPa  and  nozzle  temperatures  around  room  temperature  (Figure  10a). 
At  higher  source  pressures,  argon  dimerization  resulted  in  the  production 
of  two  peaks  in  the  time-of-fl ight  spectrum  (Figure  10b).  Dimers  have 
been  observed  before  by  others  [e.g.28]  using  mass  spectrometry ,  and 
we  also  subsequently  verified  their  presence  in  our  beam  by  mass 
spectrometry.  At  nozzle  temperatures  >  100°  C,  the  dimers  were  not 
observed. 

3.2.2.  WATER  -  GENERAL  CHARACTERISTICS  AND  VELOCITY  ANALYSIS 

The  high  pumping  speed  available  for  water  allowed  the  use  of 
source  pressures  to  over  1  MPa.  In  fact,  no  limit  on  source  intensity 
was  ever  reached  in  our  studies.  The  velocity  distribution  (without 
mass-filtering)  as  a  function  of  source  pressure  is  shown  in  Figure  11. 

A  typical  mass  spectrum  is  shown  in  Figure  12.  Stagnation  gauge  in¬ 
tensity  measurements  as  a  function  of  source  pressure  are  shown  in 
Figure  13.  A  plot  of  peak  amplitude  in  time-of-fl ight  spectra  versus 
stagnation  gauge  intensity  measurements  is  given  in  Figure  14.  At 
very  low  pressures  the  velocity  distribution  approaches  a  Maxwellian.  At 
high  pressures,  the  velocity  distribution  could  not  be  fitted  even  to  a 
translating  Maxwellian.  In  particular,  a  substantial  low-velocity 
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tail  aopears,  which  we  attribute  to  clusters  (see  Section  2.2.4).  at 
high  pressures,  mass  spectrometry  indicated  that  the  beam  was  highly 
clustered.  Substantial  superheating  (heating  the  nozzle  above  the 
boiler  temperature)  was  required  to  eliminate  all  clusters.  For  exam¬ 
ple,  for  a  boiler  temperature  of  376°  K,  clusters  were  formed  for 
nozzle  temperatures  below  about  55Qa  K. 

3.2.3.  MASS  ANALYSIS 

Referring  again  to  the  mass  spectrum  in  Figure  12,  several  inter¬ 
esting  features  should  be  noted.  Mass  peaks  were  located  at  18n  +  1  AMU 
(and  also  1,  17,  18  and  20  AMU)  in  the  range  1  s  n  s  37  which  is,  at 
present,  the  extended  range  of  the  mass  spectrometer.  Presumably  the 
pattern  would  continue  beyond  n  *  37  if  the  range  of  the  instrument 
were  extended  further.  Electron  diffraction  studies  of  beams  formed 
under  similar  conditions  [35-37]  have  indicated  clusters  of  n  *  1000 
or  more.  There  are  two  possible  explanations  for  peaks  at  18n  +  1 
instead  of  18n.  Either  an  OH'  was  removed  on  ionization  of  a  cluster 
or  an  H+  was  added.  The  mechanism  of  proton  addition  is  favored  by 
mass  spectrometri sts  [34,57]  as  an  explanation  of  our  mass  spectra. 

It  is  a  common  process  in  high  density  ionization  processes.  The  mean 
free  path  of  neutral  species  in  the  ionization  region  was  on  the  order 
of  meters.  (The  "total  pressure"  measured  in  the  ionization  region  was 
never  >  10  Pa.)  On  this  basis  collision  probabilities  would  be  ex¬ 
pected  to  be  very  low.  There  is,  however,  a  significant  attractive 
force  between  a  proton  and  a  neutral  water  molecule  due  to  the  dipole 
moment  of  the  molecule,  which  increases  the  effective  "size"  of  a 
molecule  and  decreases  the  mean  free  path  for  a  proton  -  molecule 
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collision  by  two  to  three  orders  of  magnitude.  Cn  this  basis  the 
mechanism  of  proton  addition  appears  plausible. 

It  is  interesting  to  note  that  whatever  the  ionization  process  for 
the  clusters  is,  it  appears  to  be  "clean"  for  the  low  mass  range  studied. 
We  have  been  unable  to  find  any  evidence  of  fragmentation  of  clusters 
on  ionization;  we  have  observed  no  changes  in  the  mass  spectrum  over  a 
range  of  electron  energies,  currents,  and  ion  energies.  This  is  in 
contrast  to  mass  spectrometry  of  large  molecules  which  typically 
fragment  into  many  pieces.  Half-order  peaks  appeared  for  ISh  s  n  s  3 T>$ 
these  had  a  threshhold  of  ~4Q  V  for  the  electron  energy,  suggesting 
that  they  were  due  to  double  ionization.  If  this  is  the  case,  clusters 
of  up  to  75  molecules  have  been  detected.  An  alternative  explanation 
[34]  is  the  attachment  of  contaminant  molecules  to  the  clusters.  Unless 
nitrogen  or  oxygen  or  other  permanent  gases  are  the  contaminants,  this 
explanation  seems  unlikely,  especially  in  light  of  the  prolonged  steam 
cleaning  of  the  source  assembly  through  several  months  of  use  prior 
to  the  mass  analysis. 

Lin  also  observed  the  half-order  peaks  in  his  mass  spectrometer 
data  [32]  for  the  range  1%  s  n  s  30H-  He  attributed  these  peaks  to 
the  attachment  of  a  molecule  of  silicone  pump  oil  to  the  clusters.  He 
did  not  identify  his  pump  oil,  but  claimed  a  molecular  weight  of  333 
(*  18H  x  18).  Data  from  manufacturers  on  various  diffusion  pump  oils 
give  "average"  molecular  weights  of  between  400  and  500.  A  range  of 
values  should  be  present  in  any  given  oil  sample,  so  a  value  of  exactly 
333  Is  far-fetched.  Double  ionization  remains  the  most  plausible 
explanation  for  the  half-order  peaks. 
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For  clusters  of  n  >  3,  one  looks  for  a  signal  at  13n  ♦  1  AMU.  It 
is  less  clear  where  one  should  expect  the  signal  from  a  dimer  or 
trimer,  or  to  what  original  particles  to  attribute  the  signal  at  mass  19. 
Three  sources  may  contribute  to  m  3  19:  protonated  water  monomers  and 
fragmented  dimers  (already  mentioned)  and  isotopic  species.  Ordinary 
terrestrial  water  is  known  to  consist  of  99.732  0.202 

0.042  V70.  and  0.032  ]HD160  [2].  Thus  0.072  is  mass  19  and  0.22 
is  mass  20.  With  a  high  intensity  water  source  these  species  are  quite 
detectable.  Since  the  amplitude  of  the  mass  19  signal  is  appreciably 
greater  than  0.072  of  the  mass  18  signal,  the  contribution  of  isotopes 
alone  is  not  an  adequate  explanation;  probably  the  other  two  sources  of 
m  *  19  give  rise  to  the  observed  signal. 

The  general  amplitude  trend  in  the  mass  spectrum  (Figure  12)  is  a 
monotonic  decrease  with  increasing  AMU.  The  actual  relative  amplitudes 
of  clusters  of  different  size  are  not  the  same  as  those  shown  in  Figure  12 
because  of  the  amplitude  calibration  difficulties  (see  Section  3.1.1). 
There  are,  however,  three  anomalies  in  the  plot  which  should  be  noted. 
First  is  the  relatively  large  peak  at  n  *  21  and  the  relatively  small 
one  at  n  3  22.  The  amplitude  increases  from  n  3  22  to  n  3  25.  Finally 
there  is  a  small  dip  in  amplitude  at  n  3  30.  (See  Section  4.1  for 
further  discussion.) 

3.2.4.  MASS  FILTERED  VELOCITY  ANALYSIS 

Time-of-flight  spectra  associated  with  specific  cluster  sizes  in 
the  beam  are  shown  in  Figure  15.  See  also  Figure  11  which  gives  the 
spectrum  without  mass-filtering.  For  n  >  3  (m  a  73  AMU),  the  spectrum 
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for  each  cluster  size  could  be  fitted  very  nicely  to  a  translating 
Maxwellian  (Figure  15).  In  order  to  perform  such  a  fitting,  it  was,  of 
course,  necessary  to  select  a  time  delay  (see  Section  3.1.3).  The 
delay  was  assumed  constant  for  all  cluster  sizes.  To  select  a  delay, 
the  peak  position  of  the  larger  clusters  (n  £  15)  was  matched  to  the 
peak  position  of  the  spectrum  without  mass-filtering.  This  peak  position 
also  coincided  with  the  position  of  the  larger  of  the  two  peaks  in  the 
monomer  distribution,  so  it  was  thought  that  the  combined  effect  of  a 
large  number  of  monomers  and  different  sized  clusters  at  this  particular 
time-of-fl ight  would  produce  a  peak  at  the  same  position  when  no  mass- 
filtering  was  used.  The  delay  chosen  in  this  way  was  157  us.  In 
any  case,  an  error  in  the  choice  of  time  delay  affects  the  magnitude 
of  the  temperatures  and  velocities  determined  from  the  fitting,  but 
not  the  qualitative  trends. 

Three  parameters  were  obtained  from  this  fitting  process:  ampli¬ 
tude,  velocity,  and  translational  temperature  (along  the  beam  direction; 
there  could  conceivably  be  a  different  temperature  associated  with 
thermal  motion  orthogonal  to  the  beam  or  within  a  cluster).  These 
are  plotted  in  Figures  16,  17  and  18.  Only  a  few  amplitudes  can  be 
compared,  since  to  maximize  signal-to-noise  in  time-of-fl ight  spectra, 
the  electron  multiplier  gain  was  changed  occasionally.  The  corrections 
cited  In  Section  3.1.1  apply  equally  to  both  mass  spectrometer  data  and 
mass-filtered  time-of-fl ight  data,  and  hence  the  amplitudes  in  Figure  16 
can  be  compared  directly  to  the  amplitudes  in  Figure  12.  Note,  however, 
that  the  time-of-fl ight  apparatus  without  mass  filtering  is  not  hampered 
by  quadrupole  transmission  losses.  The  low  velocity  tail  in  Figure  11, 
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which  is  evidently  due  to  clusters  beyond  the  range  we  have  been  able 
to  study  with  the  mass  spectrometer ,  will  not  be  as  prominent  in  mass- 
filtered  analysis. 

The  translational  velocity  data  shown  in  Figure  17  is  striking  in 
that  most  clusters  moved  with  nearly  the  same  speed.  Note  that  this 
speed  is  substantially  faster  than  that  which  could  be  achieved  in 
an  unclustered  expansion  by  converting  all  the  enthalpy  in  the  source 
into  directed  motion  as  indicated  in  the  figure  (see  also  Section  4.1). 

The  temperature  data  as  plotted  in  Figure  18  assumes  a  monomer  mass 
for  all  cluster  sizes.  Using  the  mass  of  the  clusters  in  the  fitting, 
one  would  find  that  the  temperature  is  proportional  to  the  cluster  size 
for  the  range  of  sizes  where  the  temperature  plotted  in  Figure  18  is 
constant!  See  Section  4.1  for  a  discussion  of  this  phenomenon.  Note  also 
that  the  temperatures  determined  in  this  way  are  somewhat  lower  than 
the  100-150°  K  suggested  by  Stein  and  Armstrong[35]  for  their  beams. 

No  explanation  is  available  for  the  discrepancy  between  the  data  taken 
in  the  two  separate  runs.  The  source  conditions  were  nearly  identical, 
and  the  experiments  must  be  repeated  and  extended  to  resolve  the  dis¬ 
crepancy. 

Only  a  small  amount  of  data  has  been  taken  thus  far  as  a  function 
of  source  pressure.  The  effects  on  the  time-of -flight  spectrum  without 
mass  filtering  were  cited  In  Section  3.2.2  and  are  shown  in  Figure  11. 

We  attribute  the  appearance  of  a  low  velocity  tail  to  the  formation  of 
large  (n  >  37)  clusters  in  the  beam.  The  double  peak  in  the  time-of-flight 
spectra  for  masses  17,  18,  19,  20  and  to  a  lesser  extent  37  and  55  also 
correlates  with  the  formation  of  large  clusters.  For  example,  raising 
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the  source  temperature  only  20°  caused  substantial  changes  as  shown  in 
Figure  19.  Apparently,  the  slower  peak  is  in  someway  associated  with 
the  formation  of  large  clusters.  See  Section  4.1  for  further  discussion. 


3.3.  EVAPORATION 

Some  interesting  information  was  obtained  by  looking  at  the  material 
emitted  by  an  evaporating  ice  surface  with  no  impinging  beam.  The  study 
of  evaporation  was  never  a  major  goal  of  the  project,  so  the  data  are 
somewhat  incomplete.  There  are  several  interesting  features  in  the 
data,  though,  so  we  present  what  we  now  know. 

A  summary  of  experimental  conditions  studied  is  given  in  Table  1. 


Table  1.  Evaporation  Oata. 


Surface 

Temperatures 

Findings 

Single  Crystal  Ice 

214-224°K 

For  low  end  of  temperature  range, 
emission  is  monomer  cosine/Maxwelli¬ 
an  at  a  temperature  10-20®  below 
surface  temperature.  At  higher  tem¬ 
peratures,  emission  is  non-Maxwell i- 
an  and  there  may  be  dimers  emitted. 

Polycrystall ine  Ice 

229°K 

Emission  is  monomer  Cosine/Maxwel 1 i - 
an  at  ~209°K. 

Deposit  on  Covellite 

204°  K 

Higher  amplitude  toward  tangent,  velo 
city  Maxwellian  at  tangent,  slower  at 
normal . 

Oeposit  on  Silicon 

210®K,  212°K 

Maxwellian  at  surface  temperature. 

Oeposit  on  Platinum 

207°K,  212°K 

Maxwellian  at  surface  temperature. 
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Two  bulk  ice  surfaces  were  used,  one  cut  f’-on  single  crystal  glacial 
ice,  and  one  from  a  polycrystalline  samole.  Also  studied  were  evapo¬ 
ration  of  accumulated  ice  layers  on  covellite,  silicon,  ana  platinum. 

It  was  only  possible  to  study  a  fairly  narrow  temperature  range  'ap¬ 
proximately  210°  K  <  T  <  225°  K) ,  because  the  vapor  pressure  of  an  ice 
surface  rises  very  rapidly  with  temperature. 

Roughly  speaking,  the  emission  from-  all  evaporating  surfaces  was 
Maxwell  ian/cosine  (i.e.  the  velocity  distribution  was  Maxwellian  at 
the  surface  temperature  and  the  spatial  distribution  was  cosine). 

There  were,  however,  some  notable  exceptions,  and  these  are  the  obser¬ 
vations  of  particular  interest. 

The  mass-filtered  time-of-flight  spectra  allowed  mass  analysis  of 

the  evaporating  particles  separately  from  the  background  [see  Section  3.1.3). 

Under  all  circumstances,  signals  could  be  found  at  mass  13  and  17.  At 

higher  target  temperatures  a  signal  appeared  at  mass  19.  No  signal  was 

ever  conclusively  identified  at  mass  37  or  higher.  The  signal  at  mass  19 

is  most  likely  due  to  dimers  (see  Section  3.2.3).  Proton  addition  is 

not  likely  to  be  important,  since  the  pressure  at  the  ionizer  is  on 
-5 

the  order  of  10  Pa.  Isotopic  contributions  are  not  ruled  out  entirely, 
because  relative  amplitudes  of  mass  18  and  19  signals  were  not  obtained. 
Dimers  may  be  formed  shortly  after  leaving  the  surface,  rather  than  be¬ 
fore,  since  the  mean  free  path  in  the  vicinity  of  the  surface  may  be 
quite  small.  For  comparison,  the  mean  free  path  at  the  equilibrium 
vapor  pressure  of  ice  at  225°  K  is  about  1  mm. 

In  almost  all  cases,  there  were  no  measurable  deviations  from  a 
cosine  spatial  distribution.  The  exception  was  the  evaporation  of  a 
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substantial  frost-like  accumulation  on  covellite.  "he  amplitude  of 
the  emitted  stream  was  significantly  higher  toward  the  tangent  than  at 
the  normal  (see  Figure  20).  This  may  simply  reflect  the  needle-like 
crystal  habit  (parallel  needles  sticking  out  from  the  substrate)  which, 
in  effect,  provided  a  large  area  of  surface  with  very  short  linear  di¬ 
mension  oriented  perpendicular  to  the  target.  Cosine  emission  from  such 
needles  may  well  give  the  unexpected  amplitudes  observed. 

Only  one  surface  showed  a  measurable  variation  in  the  velocity 
distribution  as  a  function  of  angle  (see  Figure  21).  This  was  again  the 
frost  on  covellite.  A  good  fit  to  a  Maxwellian  was  achieved  for  angles 
near  the  tangent.  For  near-normal  angles  there  was  a  significant  devi¬ 
ation  toward  slower  velocities.  If  this  is  interpreted  as  a  temperature 
change  in  a  Maxwellian  distribution,  the  normally  emitted  particles 
are  some  20°  C  colder  than  the  target.  Again,  the  variation  may  re¬ 
flect  the  crystal  habit. 

Evaporation  of  accumulation  on  silicon  and  platinum  gave  a  good 
fit  to  a  Maxwellian  at  the  surface  temperature.  Evaporation  from  bulk 
ice  samples  both  single  crystal  and  polycrystalline,  showed  a  deviation 
toward  slow  velocities  for  all  angles,  which  again  can  be  interpreted 
as  a  30°  C  cooling.  The  immediate  question  is  whether  or  not  the  surface 
temperature  was  in  fact  being  measured  correctly.  In  most  cases  the 
thermocouple  was  clamped  to  the  surface.  The  calibration  of  both  the 
thermocouple  and  time-of-fl ight  equipment  were  reliable  (at  least  to 
within  about  ±3°  C),  and  were  checked  together  by  backscattering  ex¬ 
periments  (see  Section  3. 4. 2. 2.1).  Even  when  the  thermocouple  was 
embedded  in  the  ice,  temperature  gradients  were  not  likely  to  be  a 
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problem.  At  typica’  evaporation  rates,  the  power  being  carried  away 
from  the  surface  is  only  sufficient  to  maintain  a  temoerature  gradient 
of  a  fract  on  of  a  degree  across  the  ice  sample. 

An  execution  to  this  behavior  was  noted  for  a  single  crystal  ice 
sample  at  221°  K  and  is  shown  in  Figure  22.  The  velocity  distribution 
now  had  a  distinctly  non-Maxwell ian  shape.  The  average  velocity  was 
somewhat  faster  than  the  average  of  a  surface  temperature  Maxwellian. 
Whether  this  behavior  was  due  to  surface  phenomena,  or  collisions  in  the 
gas  phase  immediately  adjacent  to  the  surface  is  not  clear  and  cannot 
readily  be  y,esolved  by  time-of-fl ight  techniques. 

3.4.  SCATTERING  DATA 

3.4.1.  Early  Studies 

In  this  category  we  include  all  data  prior  to  the  installation  of 
the  new  source  system  in  May,  1979  (see  Section  2.2).  Experiments 
involved  scatter  of  water  and  carbon  dioxide  beams  from  surfaces  of 
covellite  and  aluminum.  Argon  scattering  was  attempted,  but  the  signal 
levels  achieved  were  inadequate.  The  water  data  are  fully  superseded 
by  more  recent  data.  Source  pressures  of  less  than  10  kPa  were  used, 
and  only  the  scattering  patterns  which  we  now  associate  primarily  with 
monomer  scattering  appeared. 

The  carbon  dioxide  scattering  data  are  of  some  interest,  because 
they  are  unique.  Measurements  were  limited  to  a  room  temperature  nozzle 
and  a  source  pressure  of  ~5  kPa.  Signal  to  noise  levels  were  not  great 
(about  15:1  after  averaging),  but  some  tentative  conclusions  were  possi¬ 
ble.  The  scattering  was  primarily  diffuse.  No  deviations  from  a 
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Maxwellian  velocity  distribution  at  the  surface  temperature  were  aer'-'- 
nitely  identified.  Such  deviations  as  were  suspected  ’were  similar  to 
those  observed  with  water.  Deviations  from  a  cosine  spatial  distribution 
also  showed  patterns  like  water:  there  was  a  greater  intensity  at  the 
target  tangent  than  at  the  target  normal.  It  would  be  interesting  to 
repeat  these  experiments  with  the  higher  intensity  source,  but  to  achieve 
comparable  source  conditions  would  require  a  chilled  nozzle. 

3.4.2.  THE  WATER  SCATTERING  DATA 

3. 4. 2.1.  EXPERIMENTAL  CONDITIONS  TESTED 

We  describe  here  the  results  of  studies  made  using  a  water  beam 
formed  in  the  final  source  assembly  described  in  Section  2.2.  The 
primary  detection  system  used  was  the  time-of-fl ight  spectrometer  (with¬ 
out  mass  filtering).  Experimental  parameters  varied  were  source  pres¬ 
sure  and  temperature,  target  material,  temperature,  and  history,  and 
incident  and  scattering  angles. 

Source  pressure  was  varied  up  to  about  1  MPa.  Nozzle  temperature 
varied  little  on  an  absolute  scale.  The  range  was  from  the  boiler  temper¬ 
ature  of  370-470°  K  up  to  about  620°K.  The  main  effect  of  these  variations 
was  on  the  degree  of  clustering  of  the  beam  (see  Section  3.2). 

Six  targets  were  tested:  A1 ,  Pt,  Si,  Glass,  CuS,  and  H2O.  The 
differences  in  scattering  from  these  surfaces  are  summarized  in  Table  2 
and  described  in  more  detail  in  Section  3. 4. 2. 2.  Descriptions  of  the 
surfaces  themselves  are  given  in  Section  2.3.2.  Incident  and  scattering 
angle  effects  are  also  described  in  detail  in  Section  3. 4. 2. 2. 
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Table  2-  Surface  Temperature  Range  eor  Observed  features. 


Surface 

Maxwel 1 lan 
Peak 

fast 

Peak 

Specular 

Slow  Peak 

Non-Specu' ar 
Slow  Peak 

Cleaved  Si ( 11 1 ) 

d 

\ 

i 

! 

210°!(<Tss300oK 

2io°:-:<'ssuCC',,< 

Pol i shed  ?? 

polycrystal  1  ine 

: 

Not  present 

210°K<Tss5C0  - 

Polished  CuS 

. 

(1000) 

Ts  >27  0°  K 

Ts>210°K 

210°K<TS<3C0°.<3 

210°'<<T  r330r< 

Polished  A1 

1 

! 

i 

Present  at  2254KC 

polycrystal  1  ine 

i 

I 

?C 

Glass 

i 

210°KsT$53CC9K 

210°.K<T  S5G0°< 

H23  (1000) 

1 

f 

N 

f 

Present11 

?d 

a  No  upper  threshold  found 
b  A  weaker  feature  than  on  other  surfaces 
c  Except  as  noted,  conditions  not  checked 

d  Very  difficult  to  identify  because  of  evaporation  from  the  surface 


3. 4. 2. 2.  FEATURES  OBSERVED 

Four  distinct  types  of  scattering  were  detected.  To  facilitate 
discussion,  we  label  these  in  terms  of  their  distinguishing  charac¬ 
teristics  in  the  time-of-fl ight  spectra:  the  "fast  peak",  the  "Max¬ 
wellian  peak",  the  "specularly  directed  slow  peak",  and  the  "non-specu- 
lar  slow  peak".  Labels  which  are  descriptive  of  the  mechanism  of  scat¬ 
tering  must  await  understanding  of  the  mechanism. 

Each  of  the  four  features  cited  has  been  independently  resolved 
under  some  experimental  conditions.  Under  most  conditions  at  least 
two  features  overlapped  to  the  extent  that  there  was  no  valley  between 
them.  This  was  particularly  the  case  with  the  fast  peak  and  Maxwellian 
peak,  but  one  example  where  they  are  resolved  is  shown  in  Figure  23. 
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Also  note  the  spectrum  shown  in  "igure  25b,  where  tne  Maxwellian  pea* 
is  undetectable.  In  cases  where  the  two  peaks  are  unresolved,  it 
seems  reasonable  to  assume  on  tht  basis  of  backscattering  observations 
(where  only  the  Maxwellian  component  survives  [see  Figure  24])  that 
the  "Maxwellian  peak"  in  the  unresolved  spectrum  is  indeed  a  Maxwel¬ 
lian  at  the  surface  temperature  and  cosine  spatially,  hence  one  can 
isolate  the  fast  peak  with  some  confidence  by  subtracting  out  a  Max¬ 
wellian  at  the  surface  temperature.  The  amplitude  of  this  Maxwellian 
cannot  be  determined  with  satisfactory  precision  from  backscattering 
data,  because  there  are  substantial  drifts  in  the  gain  of  the  electron 
multiplier  with  time.  Instead,  a(  amplitude  is  typically  chosen  such 
that  the  remaining  peaks  have  plausible  shapes  and  there  is  no  dip  be¬ 
low  the  zero  amplitude  line. 

The  properties  of  the  four  peaks  are  described  individually  in  the 
following  sections.  Discussion  of  possible  scattering  mechanisms  is 
deferred  to  Sections  4.3  and  4.4. 

3. 4. 2. 2.1.  THE  MAXWELLIAN  PEAK 

As  discussed  above,  our  experiments  indicate  that  this  feature  is  Max¬ 
well  ian/cosine  (monomer)  at  the  surface  temperature  (see  Figure  24).  It  was 
found  under  almost  all  experimental  conditions,  though  it  was  sometimes 
dwarfed  by  either  the  fast  peak  or  one  of  the  slow  peaks.  Small  deviations 
from  Maxwellian  were  difficult  to  detect  because  of  the  presence  of  the  oth¬ 
er  peaks  in  a  scattering  pattern.  Such  deviations,  if  present,  would  not, 
of  course,  alter  the  analysis  of  the  remaining  peaks  appreciably.  If 
one  assumes  a  trapping/desorption  model  for  the  Maxwellian  peak,  then 
the  data  should  look  much  like  the  evaporation  of  an  accumulated  layer 


of  ice  (at  least  in  the  limit  where  the  layer  is  very  thin).  The  evacora 
tion  data  of  Section  3.3  thus  further  support  the  conclusion  that  tne 
Maxwellian  peak  is,  in  fact,  Maxwellian  at  the  su-^ace  temperature. 

The  particles  in  the  incident  beam  giving  r*s.  to  the  Maxwellian 
peak  could  be  monomers  or  clusters  or  both.  Mass  spectrometer  measure¬ 
ments  indicate  that  any  particles  which  are  trapped  at  the  surface,  re¬ 
gardless  of  their  original  size,  evaporate  as  monomers.  Some  measure 
of  the  time  scale  of  this  accommodation  process  can  be  made  by  comparing 
data  obtained  by  chopping  the  incident  beam  to  data  obtained  by  chopping 
the  scattered  stream.  The  two  kinds  of  spectra  are  identical  (taking 
into  account  the  changed  flight  path  and  the  fact  tnat  tne  velocity 
distribution  is  that  of  the  incident  beam  until  the  particles  reach  the 
target).  Evidently  the  residence  time  for  the  trapped  particles  is  less 
than  a  few  microseconds  (the  smallest  time  shift  observable).  Actually, 

this  result  should  be  expected;  since  molecular  vibrations  have  periods 
-13 

on  the  order  of  10  seconds,  a  microsecond  is  a  long  time.  A  practical 
experiment  for  measuring  residence  time  must  be  able  to  measure  such 
shorter  times. 

There  was  one  experimental  condition  where  a  notable  difference  was 
observed  between  the  spectra  obtained  by  chopping  the  incident  and  scat¬ 
tered  streams.  The  pair  of  spectra  shown  in  Figure  25  are  for  a  water 
beam  scattering  from  a  silicon  (111)  surface  at  212°  K.  In  this  case, 
the  target  was  cold  enough  that  most  of  the  beam  stuck  indefinitely  to 
the  surface.  The  few  particles  reemitted  as  a  Maxwellian  stream  after 
sticking  lost  correlation  with  the  incident  beam  (i.e.  have  resided 
more  than  several  hundred  microseconds  on  the  surface).  No  intermediate 
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target  temoerature  range  giving  a  measurable  residence  time  could  be 
found.  The  behavior  appeared  to  be  bimodal .  Either  particles  stuck 
for  a  short  time  before  remission,  or  they  stuck  for  some  indeterminate 
(and  random)  long  time. 

The  deviations  from  Maxwellian  behavior  noted  in  the  early  results 
(Section  3.4.1)  can  all  be  attributed  to  a  manifestation  of  a  fast 
peak  which  overlaps  the  fast  end  of  the  Maxwellian  peak.  The  signal-to- 
noise  ratio  in  those  data  was  not  high  enough  to  separate  the  fast  peak 
satisfactorily.  Hindsight  and  a  bit  of  imagination  show  the  peak  to 
be  present. 


3. 4. 2. 2. 2.  THE  FAST  PEAK 


The  behavior  of  the  fast  peak  for  typical  operating  conditions  is 
★ 

shown  in  Figures  26  and  27.  The  major  features  are  a  small  but  measura¬ 
ble  decrease  in  speed  from  that  of  the  incident  beam  and  a  scattering 
pattern  in  which  intensity  increases  from  the  surface  normal  to  the 
tangent  (no  well  defined  lobe),  regardless  of  incident  angle  (0°  s 
s  75°).  Furthermore,  the  fast  peak  is  approximately  constant  in  ampli¬ 
tude  for  constant  “turning  angle",' i.e.  for  9^  +  9S  *  constant. 


*  In  this  and  all  subsequent  plots  of  peak  amplitude,  the  Maxwellian 
peak  amplitude  is  used  for  normalization.  The  detector  collimation 
is  such  that  for  angles  less  than  about  75°,  the  detector  sees  an  in¬ 
finite  target.  The  Maxwellian  cosine  distribution  is  then  constant 
with  scattering  angle.  The  flux  per  unit  area  on  the  surface  is  a 
function  of  incident  angle,  so  the  Maxwellian  peak  intensity  is  not 
constant  for  changes  in  incident  angle.  By  dividing  by  the  Maxwellian 
amplitude,  this  variation  and  any  experimental  drifts  are  eliminated. 
Note  that  since  peak  amplitudes  and  not  integrals  are  used,  the  ratios 
do  not  represent  the  relative  numbers  of  particles  even  for  constant 
ionization  efficiency. 
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"he  velocity  loss  also  varies  somewhat  with  angle  and  also  appears 
constant  for  constant  turning  angle.  This  peak  was  observed  for  all 
targets  at  temperatures  above  the  point  where  virtually  all  beam  parti¬ 
cles  stick  to  the  surface  (around  210°  K  for  typical  beam  intensities). 
Even  the  amplitude  and  velocity  of  the  peak  appear  to  be  independent 
of  target  conditions,  at  least  for  the  cases  which  have  been  checked, 
and  this  peak  is  the  last  to  disappear  as  the  target  temperature  is 
lowered.  Even  when  nearly  everything  is  sticking  a  few  particles  con¬ 
tinue  to  bounce  off  with  little  loss  of  energy. 

3. 4. 2. 2. 3.  THE  SPECULARLY  DIRECTED  SLOW  PEAK 

A  striking  feature  in  the  scattering  data  is  the  appearance  under 
a  fairly  narrow  range  of  experimental  conditions  of  a  specularly  direct¬ 
ed  slow  peak.  This  peak  may  dominate  a  time-of-fl ight  spectrum  and  vary 
widely  in  peak  position  and  amplitude  as  shown  by  the  examples  in 
Figure  28.  If  attributed  to  water  monomers,  its  average  velocity 
would  be  lower  than  a  Maxwellian  at  any  surface  temperature  in  the 
system  (the  coldest  being  77°  K) .  The  intensity  pattern  is  a  lobe  at 
the  specular  angle  (at  least  to  within  +5°)  with  a  half  width  of  20° 

(at  least  for  cleaved  Si  (111)  for  which  a  detailed  survey  was  made). 
Figure  29  shows  this  as  a  plot  of  the  slow  peak  amplitude  versus  scat¬ 
tering  angle.  A  plot  of  stagnation  gauge  flux  measurements  for  similar 
conditions  is  shown  In  Figure  30.  (Since  the  slow  peak  can  be  made  to 
dominate  the  scattering,  the  stagnation  gauge  primarily  measures  the 
slow  peak  amplitude.)  Note  that  the  stagnation  gauge  and  flow  through 
Ionizer  may  detect  clusters  of  different  efficiencies.  The  peak  shape 
is  non-Maxwell  Ian  and,  in  fact,  resembles  fairly  closely  the  incident 
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beam  distribution  scaled  in  time  by  a  factor  of  2  to  5  (see  rigure  31}. 

The  peak  varied  in  amplitude  but  not  in  position  as  a  function  of 
scattering  angle  as  expected  for  lobular  scattering.  There  were  pro¬ 
nounced  changes  in  both  amplitude  and  peak  position  as  a  function  of 
incident  angle  (Figures  29  and  32).  The  peak  has  not  been  detected 
for  incident  angles  less  than  50°.  The  arrangement  of  the  apoaratus 
limited  satisfactory  data  to  incident  angles  less  than  ~75°. 

As  indicated  in  Table  2,  some  evidence  for  the  specularly  directed 
slow  peak  was  found  with  all  surfaces  tested  except  platinum,  and  for 
target  temperatures  in  the  range  210°  K  <  T  <  300°K.  Below  210°  K, 
all  particles  in  the  incident  beam  stick  to  the  target  (though  as 
the  temperature  is  lowered,  the  slow  peak  disappeared  at  a  slightly 
higher  temperature  than  the  fast  peak  and  Maxwellian  peak).  The 
strongest  slow  peak  amplitude  was  achieved  in  the  range  240°  K  <  T 
<  245°  K  (at  least  for  Si  (111));  as  the  temperature  was  raised  above 
~300°  K,  the  peak  disappeared.  The  effect  of  target  temperature  on 
the  slow  peak  position  (velocity)  was  relatively  small  (Figure  33). 
However,  in  the  available  temperature  range,  temperature  could  only 
vary  by  about  15  %  on  an  absolute  scale,  so  major  effects  were  not  to 
be  expected. 

The  behavior  of  the  specularly  directed  slow  peak  as  a  function 
of  nozzle  temperature  and  source  pressure  provides  the  strongest  evi¬ 
dence  for  the  importance  of  large  clusters  in  the  beam.  The  effect 
can  be  described  in  two  ways.  For  constant  source  pressure  of  350  kPa, 
the  slow  peak  disappears  when  the  nozzle  temperature  was  only  40-50°  C 
above  the  boiler  temperature.  For  constant  nozzle  temperature,  there 


was  a  threshold  source  pressure  (bailer  temperature) ,  below  which  the  slew 
peak  could  not  be  found  [see  Table  3).  Because  of  the  rapid  rise  in  vapor 


Table  3.  Source  Pressure  Thresholds  for  the  Appearance  of  the  Specularly 
Directed  Slow  Peak 


Source 

Temperature  (T  ) 

Threshold 

Pressure  (PQ)a 

Relative  . 

Beam  Flux0 

455°  K 

350  kPa 

1.0 

485 

550 

1.4 

535 

800 

1.5 

605 

1300 

1.9 

a  No  consistent  criterion  was  used  to  define  threshold-  Values  good  to 
within  about  ±15%. 

b  Conversion  from  stagnation  gauge  pressure  versus  source  pressure  measur- 
ments.  Probably  gives  a  good  indication  of  relative  mass  flux  (see  Sec¬ 
tion  3.2.2  and  Figure  13). 


pressure  of  water  with  boiler  temperature,  the  source  pressure  threshold 
for  the  onset  of  the  slow  peak  increased  approximately  as  the  square  of 
the  nozzle  temperature,  and  the  degree  of  superheating  necessary  to  elimi¬ 
nate  the  slow  peak  increased  also:  at  1.3  MPa,  a  superheating  of  about 
140°  C  is  needed.  Note  also  that  for  constant  nozzle  temperature,  the 
fast  peak  and  Maxwellian  peak  amplitudes  rose  approximately  linearly  with 
source  Intensity,  while  the  specularly  directed  slow  peak  amplitude  rose 
much  more  rapidly  (see  Figure  33). 

Finally,  there  is  an  interesting  time  dependence  in  the  behavior  of 
the  slow  peak  amplitude.  (Only  the  amplitude  was  affected  in  observations 


to  date.)  If  a  target  was  heated  to  ~4GQ0  K  to  drive  off  adsorbed  water, 
cooled  to  2^0°  k,  then  exposed  suddenly  to  a  glancing  incidence  team 


(a.  >  50°) ,  it  took  minutes  to  hours  for  tne  slow  peak  to  reach  its  full 
amplitude  from  a  near  zero  initial  amplitude.  The  process  can  be  accele¬ 
rated  by  using  a  normally  incident  beam,  for  which  about  15  minutes  suf¬ 
fices  to  achieve  steady-state  surface  conditions.  Apparently  the  slow 
peak  depends  on  the  presence  of  a  surface  water  or  ice  layer  which  forms 
slowly  in  a  temperature  range  where  a  macroscopic  layer  does  not  accumu¬ 
late.  (The  vapor  pressure  of  bulk  ice  at  240°  K  causes  bulk  ice  to  evap¬ 
orate  much  faster  than  could  be  resupplied  by  the  beam.)  No  surface 
monitoring  instruments  besides  the  beam  (such  as  IEED,  Auger  Electron 
Spectroscopy,  or  a  quartz  crystal  microbalance)  were  available,  so  one 
can  only  speculate  on  the  conditon  of  the  surface.  Certainly  no  macro¬ 
scopic  (visible  to  the  eye)  ice  layers  were  present. 

3. 4. 2. 2. 4.  THE  NON-SPECULAR  SLOW  PEAK 

As  source  intensities  were  Increased  above  those  used  for  most  of 
the  studies  of  the  specularly  directed  slow  peak  (PQ  *  600  kPa,  Tq 
»  450°  K),  a  new  peak  was  detected  in  the  scattered  streams.  The  main 
features  of  this  peak  were  a  velocity  comparable  to  that  of  the  specularly 
directed  slow  peak  and  a  tangential  reemission  angle  for  all  incident 
angles  (55°  s  8^  s  80°)(see  Figure  34).  The  peak  shape  of  the  time-of- 
flight  spectrum  was  non-Maxwell  Ian  but  broader  than  that  of  the  specularly 
directed  slow  peak  (see  Figure  35).  Under  conditions  where  both  slow 
peaks  were  present,  it  was  difficult  or  impossib...  to  separate  them.  For¬ 
tunately,  there  were  conditions  where  each  could  be  studied  separately. 

Exact  intensity  thresholds  for  the  appearance  of  the  non-specular 
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slow  peak  were  not  established,  but  source  pressures  of  at  least  1.3  MPa 
were  required,  suggesting  that  this  non-specular  slow  peak  is  associ¬ 
ated  with  larger  clusters. 

The  peak  position  as  a  function  of  incident  and  scattering  angles 
C for  the  range  55°  to  80°  where  detection  was  possible}  is  shown  in 
Figure  36.  The  variation  in  amplitude  is  quite  large  and  similar  to 
that  of  the  specularly  directed  slow  peak.  Note,  that  unlike  the  specu¬ 
larly  directed  slow  peak,  there  are  amplitude  variations  as  a  function 
of  both  incident  and  scattering  angles,  analogous  to  those  of  the  fast 
peak. 

The  non-specular  slow  peak  was  found  in  scattering  from  silicon, 
platinum,  glass,  and  covellite.  Aluminum  and  ice  were  not  checked. 

There  were  striking  variations  in  its  behavior  as  a  function  of  target 
temperature  (see  Table  2).  On  silicon  both  slow  peaks  were  detected 
only  below  280°  K.  On  covellite,  the  non-specular  slow  peak  persists 
to  about  330°  K;  on  glass,  the  threshold  was  about  470°  K;  and  on 
platinum  the  slow  peak  showed  no  sign  of  extinction  at  585°  K  (no  higher 
temperatures  were  attempted).  On  glass  there  was  an  observable  hystere¬ 
sis:  after  heating  the  target  until  the  slow  peak  vanished,  it  was 
necessary  to  cool  it  to  near  273° K  before  the  slow  peak  reappeared. 

The  implication  is  that  there  is  again  a  surface  coverage  effect,  though 
the  rank  ordering  of  the  threshold  temperatures  is  mystifying. 

Finally,  it  is  interesting  to  note  that  the  non-specular  slow  peak 
(and  presumably  also  the  specularly  directed  slow  peak,  though  this  was 
not  checked)  vanished  completely  when  the  inside  of  the  nozzle  became 
partially  blocked.  The  source  pressure  ratio  and  nozzle  temperature 
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remained  as  before,  but  the  reduced  intensity  evidently  reduced  the 
clustering  in  the  expansion  from  the  nozzle. 

3. 4.2. 3.  SCATTERING  FROM  ICE 

Nearly  all  beam  particles  stuck  to  the  surface  for  target  temper¬ 
atures  where  the  ice  surface  was  not  evaporating  excessively.  Even  when 
chopping  the  incident  beam,  high  evaporation  rates  destroyed  the  signal, 
because  the  local  density  became  high  enough  for  significant  gas-phase 
collisions.  With  some  careful  searching,  it  was  possible  to  find  a 
temperature  where  a  noisy  signal  was  achieved  and  the  fast  peak  was 
visible.  Chopping  the  scattered  beam,  of  course,  showed  mostly  simple 
evaporation.  With  a  bit  of  imagination  it  was  possible  to  identify  a 
slow  peak  under  appropriate  conditions,  but  the  evaporating  particles 
domi nated. 

It  does  not  appear  likely  that  it  will  ever  be  possible  to  observe 
the  scattering  of  water  from  ice.  Possibly  a  higher  intensity  source 
would  allow  the  use  of  a  lower  target  temperature,  but  then  it  would 
be  difficult  to  avoid  high. densi ties  caused  by  the  beam  itself.  Clear¬ 
ly  the  ice  layers  believed  to  be  on  the  higher  temperature  targets  must 
be  more  tightly  bound  to  the  underlying  substrate  than  are  the  surface 
molecules  of  bulk  ice. 

3. 4.2. 4.  MASS-FILTERED  VELOCITY  ANALYSIS 

The  mass-filtered  version  of  the  time-of -flight  equipment  is  not 
yet  set  up  to  do  detailed  spatial  scans.  It  can  function  in  two  fixed 
positions,  one  along  the  incident  beam  axis  and  the  other  at  an  angle 
of  144°  from  the  source  target  line,  thus  limiting  scattering  studies 
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to  a  narrow  range  of  angles.  Calibration  of  tnis  equipment  [see  Section 
3.1.3)  suggested  that  it  is  generally  reliable  but  requires  correction 
for  delay.  Results  obtained  from  observations  of  evaporation  reinforced 
confidence  in  the  system. 

Results  of  scattering  from  silicon  under  conditions  where  the  non- 
specular  slow  peak  was  present  are  shown  in  Figure  37.  Signals  were 
found  for  masses  18,  19,  37  and  no  higher.  This  is  consistent  with 
our  hypothesis  that  the  slow  peaks  are  due  to  clusters  of  masses  beyond 
the  range  of  the  mass  spectrometer.  Unfortunately  the  time-of-fl ight 
spectra  obtained  at  low  mass  (18,  19  AMU)  (Figure  37)  resemble  neither 
the  incident  beam  (Figure  15)  nor  the  spectra  obtained  without  mass¬ 
filtering  (Figure  38).  Correcting  for  delay  (see  Section  3.1.3),  one 
finds  that  the  first  peak  in  Figure  37  probably  corresponds  to  the  fast 
Maxwellian  peak  in  Figure  38.  The  second  peak  must  be  a  manifestation 
of  the  slow  peak  in  Figure  38,  though  the  peak  shape  is  distinctly 
different.  As  shown  in  Figure  37,  this  second  peak  did  appear  when  and 
only  when  the  source  conditions  were  such  that  the  non-specular  slow  peak 
appeared.  Without  more  data  at  a  higher  mass  range,  it  is  difficult  to 
interpret  the  monomer  time-of-fl ight  spectra  in  Figure  37.  The  presence 
of  a  slow  peak  in  the  monomer  time-of-fl ight  spectrum  does  not  necessarily 
contradict  the  hypothesis  that  the  slow  peak  is  due  to  clusters.  Monomers 
could  be  detached  from  clusters  in  the  ionizer. 
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SECTION  IV 


DISCUSSION  AND  GUIDELINES  FOR  FUTURE  WORK 
4.1.  HOMOGENEOUS  NUCLEATION  IN  THE  SEAM 

There  is  a  good  deal  of  interesting  information  about  clusters  in  jet 
expansions  both  in  the  mass  spectrometer  and  time-of-fl ight  measurements 
of  the  beam  and  in  the  scattering  data.  Unfortunately ,  we  are  still  at 
the  stage  where  we  are  using  the  beam  to  probe  the  surface  and  the  surface 
to  probe  the  beam  at  the  same  time,  and  detailed  interpretation  of  the  data 
is  difficult. 

It  is  fairly  safe  to  say  that  the  beam  typically  consists  of  clusters 
of  all  sizes  from  monomer  out  to  at  least  a  few  hundred  molecules.  Elec¬ 
tron  diffraction  measurements  of  beams  generated  similarly  to  ours  [35-37] 
indicate  that  clusters  of  several  thousand  molecules  may  be  present  under 
some  circumstances.  From  the  magnitude  of  the  slow  peaks  we  observe  in 
scattering,  and  from  the  behavior  of  stagnation  gauge  intensity  and  time- 
of-flight  peak  amplitude,  it  appears  that,  at  least  for  higher  incident 
beam  intensities,  a  substantial  portion  of  the  mass  of  the  beam  may  be  in 
large  clusters.  It  is  only  possible  to  construct  an  approximate  cluster 
size  distribution  with  any  data  that  now  exist  on  cluster  beams  (see 
Section  3.1.1).  One  can  still  draw  some  conclusions  about  the  behavior 
of  clusters  in  two  size  ranges:  (1)  clusters  less  than  50-100  molecules, 
where  most  molecules  reside  on  the  surface  and  few  in  the  bulk  (consider, 
for  example,  a  simple  cubic  structure  of  5x5x5*125  molecules:  98  are 
surface  molecules,  and  only  27  are  bulk),  and  (2)  clusters  near  the  cri¬ 
tical  nucleus  size  for  homogeneous  nucleation,  thought  to  be  around  35- 
1000  molecules  for  water  [9,10]  . 

*  The  notion  of  a  critical  nucleus  is  that  clusters  of  lesser  size  tend  to 
shrink  on  the  average,  while  clusters  of  larger  size  tend  to  grow.  Theor- 
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The  ^formation  about  the  smaller  clusters  is  derived  *rom  mass  scec- 
ra  (Figure  12,  for  exemple;  and  mass-fi i tered  velocity  analysis  (ri:vre  ’ 


No  clusters  within  the  range  of  the  mass  spectrometer  (n  <  75)  of 
more  than  2  or  perhaps  3  molecules  survive  any  collision  with  a  surface* 
They  are  apparently  fully  destroyed  and  reemitted  as  monomers.  They  do, 
however,  survive,  at  least  in  part,  the  collision  with  a  70  eV  ionizing 
electron  or  perhaps  a  free  proton.  If  electron  impact  is  the  ionization 
mechanism,  at  least  an  OH  must  be  detached  leaving  a  free  proton  in  the 
cluster.  The  attachment  of  a  proton  traveling  at  low  velocity  relative 
to  the  cluster  (glancing  collision:  the  protons,  being  formed  from  beam 
particles,  move  with  the  beam)  could  also  be  the  ionization  mechanism. 

In  either  case,  the  proton  might  stabilize  a  cluster  by  interaction  be¬ 
tween  it  and  the  water  dipole  moments.  The  charged  clusters  clearly 
do  survive  2  kV  acceleration  in  the  rf  fields  of  the  quadrupole  to  form 
the  mass  spectra  we  observe  (else  we  would  see  no  clusters). 

One  charge  center  may  stabilize  a  cluster,  but  two  may  destroy  it. 
The  repulsion  between  the  two  charge  centers  may  tear  a  small  cluster 
apart.  This  might  explain  the  absence  of  double-ionization  peaks  (see 


etical  values  for  the  size  of  the  critical  nucleus  calculated  thermodyna¬ 
mically  by  considering  the  balance  between  surface  and  bulk  energies  of 
the  growing  cluster  in  a  quasi equi 1 ibri urn  system  (a  bad  approach  for  very 
small  clusters).  The  actual  value  depends  on  temperature  and  the  "super- 
saturation  ratio"  (the  ratio  of  the  actual  pressure  to  the  pressure  at 
which  solid  and  vapor  would  be  in  equilibrium  at  the  system  temperature) . 
One  must  use  some  care  in  extending  the  notion  of  critical  nucleus  to  the 
highly  non-equilibrium  conditions  in  a  nozzle  beam  and  to  beam-surface 
collisions. 

*  More  precisely,  we  have  not  observed  any  clusters  in  the  scattered 
stream.  While  we  have  not  surveyed  all  possible  e,  and  9  ,  nor  looked 
for  scattered  particles  out  of  the  principal  plane,  we  have  looked  at 
near  tangential  9,  and  9  where  we  would  expect  cluster  survival  to  be  . 
most  favored.  (The  normal  component  of  momentum  is  smaller  for  large 
9 j ,  and  one  would  expect  the  normal  component  to  be  the  determining  factor 
fAr  cluster  survival /destruction.) 
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Section  3.2.3)  in  the  mass  spectrum  below  n  *  13%.  Either  clusters 
larger  than  n  3  37  are  more  stable,  or  the  charge  centers  are  more  effec¬ 
tively  shielded  from  one  another. 

All  small  clusters  travel  with  nearly  the  same  velocity.  The  simolest 
explanation  for  this  observation  is  that  the  larger  clusters  are  swept 
along  by  the  monomers.  This  explanation  fails,  because  it  also  predicts 
that  all  clusters  should  have  the  same  temperature,  which  is  contrary  to 
our  results  (Figure  18).  The  clusters  could  be  formed  by  glancing  col¬ 
lisions  between  molecules  and  clusters  having  very  nearly  zero  relative 
velocity  (both  directed  and  thermal  translational). 

A  plausible  model  can  be  constructed  from  the  data  on  the  variations 
in  the  velocity  distribution  of  small  clusters  of  different  n  values. 
Consider  first  the  energetics  of  cluster  formation.  We  start  with  mono¬ 
mers  at  PQ  and  TQ  in  the  source.  Their  available  energy  (ignoring  vis¬ 
cous  losses  in  the  nozzle)  is  the  enthalpy  of  the  source  gas,  cpTQ. 

We  have  detailed  mass-filtered  time-of-fl ight  data  for  clusters  in  the 
beam  for  only  one  source  condition:  Tp  =  458°  K.  At  this  temperature 
rotational  modes  in  water  are  fully  excited  and  vibrational  modes  are 
not.  The  internal  energy  per  molecule,  u,  is  thus  3/2kTfl  +  3/2kTp  3  3kTQ 
(translational  and  rotational).  The  enthalpy  per  molecule,  h=u+Pv 
(3u+kTQ  for  an  ideal  gas)  is  thus  h34kTQ.  Tabulated  values  for  cp[59], 
actually  give  h  3  4.17  kTp  at  TQ  3  458°K. 

We  begin  by  assuming  that  each  cluster  of  n  water  molecules  is  formed 
from  a  group  of  monomers  each  of  mass  m  which  interact  only  among  them¬ 
selves.  The  energy  balance  for  this  process  can  be  written  as 
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(4.1) 


? 

n4.17kT  a  '>s(nm)u  +  £„  +  £.  +  £r 

o  tr  i  f 

2 

wnere  *s(nm)u  is  the  kinetic  energy  of  the  directed  motion  of  the  cluster, 
Etr  12  the  additional  energy  associated  with  thermal  translational  motion, 

E.j  includes  all  other  thermal  energies  (rotation  and  vibration),  and  Z, 
is  the  energy  of  formation  of  the  cluster. 

We  measure  7Q,  u,  Etr  and  n  (see  Figures  16-18),  while  E.  and  Ef 
are  unknown.  Note  that  since  velocity  distributions  for  small  clusters 
(3  <  n  <  50)  are  all  nearly  Maxwellian  and  have  half-widths  which  are  nearly 
independent  of  n  (Figure  15),  the  apparent  translational  temperature, 

Ttr,  is  proportional  to  n;  i.e.,  the  fitting  of  the  same  translating 
Maxwellian  velocity  distribution,  v2e'M^v'u^2//2kTtr,  to  all  the  data  of 
the  type  shown  in  Figure  15  implies  that  the  exponential  is  constant. 

Since  M  =  nm,  we  must  have  T  s  nT  ,  where  T  is  a  temperature  we  shall 

Ui  c  c 

refer  to  as  the  "cluster  temperature".  We  clearly  do  not  have  an  equi- 
partition  of  thermal  energy  in  the  cluster  beam.  We  have  no  way  of  pre¬ 
dicting,  then,  how  much  energy,  E^  is  contained  in  internal  modes.  One 
can  speculate  that  internal  modes  will  be  populated  according  to  a 
temperature  T  (nT  would  destroy  a  sizeable  cluster),  but  one  needs  to 

v  t 

postulate  a  cluster  structure  to  count  available  modes,  anyway,  and  we 
are  not  yet  in  a  position  to  do  so.  One  way  to  explain  the  apparent 
variation  in  Ttr  with  cluster  size  is  to  postulate  that  clusters  are  formed 
only  from  monomers  moving  together  with  virtually  no  relative  motion, 
i.e.,  monomers  whose  directed  velocities  and  thermal  translational 
velocity  are  within  the  same  narrow  range.  The  thermal  motion  of  a  cluster 
would  then  be  determined  by  a  coherent  addition  of  the  monomer  motion 
rather  than  the  more  typical  incoherent  (random)  addition.  This  implies  a 
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very  weak  binding  energy  for  the  clusters  ( insufficient  attractive  force 
to  pull  together  molecules  with  significant  relative  motion),  a  notion 
which  is  supported  by  the  observed  vulnerabil ity  of  clusters  to  impact 
(see  Section  3. 4. 2. 4).  It  might  be  more  appropriate  to  think  of  the  width 
of  the  velocity  distribution  as  arising  from  a  statistical  spread  in  the 
directed  velocity  rather  than  from  residual  thermal  motion. 

Accordingly,  E  *  34k(nTcK  and  we  have 
£.+£*  ~ 

-~-4-  *  4.17kTQ-  3/2mu  -  3/2kT.  .  (4.2) 

T  ,  u , and  T  can  be  obtained  from  the  data,  hence  everything  on  the  right 
hand  side  is  known.  We  can  thus  calculate  a  sort  of  potential  energy 
per  molecule  (consisting  of  internal  energy  and  bond  energy)  as  a  function 
of  n. 

In  order  to  determine  u  and  T  ,  we  must  set  a  time  origin  for  the 

U 

mass-filtered  time-of-flight  spectra.  Insofar  as  the  absolute  time  origin 
is  uncertain,  we  must  be  cautious  in  assigning  absolute  values  to  u 
and  T  .  It  is  reasonable  at  this  point  to  postulate  that  the  mass  spectro¬ 
meter  introduces  a  constant  delay  to  the  velocity  spectra  which  is  greater 
than  zero  (see  Section  3.1.3).  On  inspection  of  the  spectra  of  Figure  15,  it 
seems  reasonable  to  match  the  peak  position  of  the  larger  peak  in  the 
monomer  time-of-flight  spectrum  to  the  peak  of  the  spectrum  without 
mass-filtering  (Figure  15a).  This  is  also  the  peak  position  of  the  time-of- 
flight  spectra  for  n  >  15,  so  the  largest  part  of  the  mass  of  the  beam 
probably  has  this  peak  velocity.  This  match  gives  a  delay  of  157  us,  which 
is  both  the  best  estimate  and  a  maximum  reasonable  value  (all  other  peaks 
appear  at  shorter  time). 
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E,*Ef  £f 

Table  4  gives  values  for  .  For  comparison,  note  that  —  for 

bulk  ice  is  about  -580  meV/molecule.  The  main  effect  of  an  error  in  the 

time  origin  determination  would  be  just  to  shift  all  values  of  (E^+E^/n 

by  a  constant  amount.  The  trend  would  still  be  in  the  same  direction: 

larger  clusters  more  weakly  bound  than  smaller  clusters  (for  n  <  40). 

Particularly  disturbing,  though,  is  the  negative  value  for  £.  +  E^ 
for  n  *  1 .  for  a  monomer  is  zero  by  definition.  The  extra  energy 
for  monomers  is  probably  a  result  of  collisions  between  monomers  and  the 
larger  clusters,  which  we  have  neglected  to  consider  thus  far.  The 
formation  of  the  bond  involved  in  the  growth  of  the  larger  clusters 
may  well  require  the  presence  of  a  third  body  to  remove  some  excess 
energy,  so  it  is  not  reasonable  to  neglect  such  processes. 

The  amount  of  energy  exchanged  between  monomers  and  clusters  of  a 
particular  size,  n,  cannot  be  found.  Now,  instead  of  writing  an  energy 
balance  for  a  single  cluster,  the  best  we  can  do  is  to  write  an  energy 
balance  for  the  system  as  a  whole.  Equation  4.2  must  be  rewritten  as 
follows: 

A„  -  V*A  V2m[u(n)]2 
n  n 

-3/2kT,ZAn  (4'3) 

c  n  n 

where  A  Is  the  number  of  clusters  of  size  n  and  £  nA„  is  the  total  num- 
n  n  n 

ber  of  molecules.  This  expression  would  give  us  a  weighted  average  of 
the  potential  energy  per  molecule,  provided  we  could  determine  the  An's 
with  reasonable  accuracy.  One  may  be  able  to  determine  the  An's  with 
sufficient  accuracy  by  using  the  corrected  peak  heights  from  the  mass 
spectrum.  To  obtain  useful  information  in  this  way  requires  more  data 


Tk 

(—>  n 


4.17kT 


X 


67 


n 

n 

1 

-202  meV/molec 

2 

-177 

3 

-146 

4 

-12Q 

5 

-97 

6 

-97 

7 

-30 

8 

-80 

9 

-80 

10 

-62 

11 

-62 

12 

-65 

13 

-65 

14 

-65 

15 

-51 

16 

-51 

a  Calculated  from  equation  4.2 
using  data  of  Figure  15. 

cf.  binding  energy  of  water  in 
bulk  ice  »  580  meV/molec 
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with  different  distributions  for  A  (different  T  and  P  ).  Using  such 

n  oo 

data  it  might  be  possible  to  make  some  interesting  statements  about 
c  +£ 

“i  f  as  a  function  of  cluster  size.  Of  course,  it  would  be  very  useful 
n 

to  resolve  the  time  origin  problems  and  extend  the  useful  mass  range  of 
the  mass  spectrometer. 


Perhaps  a  more  fruitful  approach  is  to  look  directly  at  the  Afl.  In 
the  range  of  cluster  sizes  (below  the  critical  nucleus  size)  which  have 
had  sufficient  opportunity  to  grow  and  diminish  by  collision  (this  might 
include  all  but  the  largest  clusters  in  the  beam  at  any  particular  source 
conditions),  one  might  expect  the  number,  An,  of  clusters  of  size  n  to 
be  distributed  as  a  function  of  the  energy  of  formation,  Ef,  according  to 

An  *  Aie"EfAT  •  (4.4) 

Probably,  the  appropriate  temperature  is  T  *  T  .  One  can  thus  obtain 
Ef(n)  directly  from  the  Ap.  Furthermore,  one  can  attempt  to  fit  E^(n)  to  a 
function  of  the  form 


3  *  -*»* 


+  Bn 


(4.5) 


where  A  represents  the  surface  energy  per  molecule  of  a  cluster,  and  8 
represents  bulk  energy.  This  would  be  the  first  experimental  determination 
of  A  and  8.  Note  that  one  also  obtains  a  measure  of  critical  nucleus 

O  A 

size  from  A  and  8:  nc*(-3$)  •  Such  a  fitting  is  shown  in  Figure  12.  We 
obtain  A»1.56,  8-. 29,  and  n  »46.  Note  that  these  energies  are  an  order  of 

t 

magnitude  smaller  than  those  obtained  by  assuming  bulk  crystalline  values. 
[10]  See  Appendix  E  for  a  discussion  of  the  effect  of  mass  spectrometer 
amplitude  correction  factors  on  this  fitting.  It  would  be  desirable  to 
collect  data  for  a  variety  of  source  conditions,  of  course,  and  we  are 
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still  in  the  process  of  refining  this  approach. 

The  anomalies  in  the  mass  spectrum  must  also  contain  potentially 
useful  information  about  the  smaller  clusters.  For  n<4G,  say,  there  is 
no  reason  to  expect  clusters  with  bulk  ice  structure.  Indeed,  one  theory 
[6,10]  suggests  a  variety  of  other  more  likely  (or  energetically  favorable) 
structures.  The  differences  in  free  energy  between  these  structures  and 
those  of  bulk  ice  are  not  large,  and  alternate  structures  may,  in  fact, 
coexist.  The  feature  at  n*21  has  now  been  observed  by  several  groups 
[32-34].  One  proposed  structure  [59]  is  a  regular  dodecahedron  (20 
vertices)  with  a  water  molecule  trapped  in  the  center.  Only  a  small 
distortion  of  the  tetrahedral  bond  angle  of  hexagonal  ice  is  necessary 
(  108°  instead  of  109°28').  The  apparent  dip  in  the  mass  spectrum  for 
clusters  of  22-24  molecules  can  also  be  explained  by  postulating  that 
at  least  some  of  these  structures  are  unstable  compared  to  the  simple 
clathrate  and  lose  molecules  to  achieve  n*21 .  Our  model  of  weakly  bound 
clusters  is  not  consistent  with  this  picture,  however,  and  we  do  not 
yet  have  a  satisfactory  alternative. 

The  feature  at  n=30  has  not  been  previously  reported,  but  is  quite 
reproducible.  We  have  no  good  explanation  for  its  presence. 

For  homogeneous  nucleatlon,  the  main  interest  in  cluster  structure 
is  to  establish  the  range  of  contributing  configurations  and  the  paths 
among  them.  From  the  point  of  view  of  nucleation  and  growth  of  macro¬ 
scopic  crystals,  the  most  critical  events  take  place  at  cluster  sizes 
close  to  the  critical  nucleus  size.  With  clusters  of  smaller  size, 
it  is  only  Important  to  know  that  there  are  no  important  barriers  to 
stable  nucleation  or  phenomena  which  affect  the  macroscopic  nucleation 
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rata.  The  20  or  21  molecule  clathrate  structure  is  potentially  such  a 
oarrier,  but  present  experimental  evidence  suggests  that  the  effect  is 
small.  Perhaps  the  clathrate  structure  cannot  even  become  a  critical 
nucleus,  but  must  first  break  up  and  re-form  in  another  configuration, 
in  which  case,  it  is  not  important  unless  it  ties  uo  a  significant  amount 
of  mass. 

Detailed  information  about  large  clusters  is  even  more  difficult 
to  obtain.  Specially  designed  mass  spectrometers  have  approached  the 
ten  thousand  molecule  mark  for  and  [34,39,43],  but  mass  spectrometry 

in  this  size  range  is  difficult.  Our  instrument  has  been  extended  to 
(^0)37  (or  (1120)75,  if  half-order  peaks  are  interpreted  as  due  to  double 
ionization).  Electron  diffraction  measurements  [35-37]  are  useful  for 
the  larger  n  values,  but  are  limited  in  the  sort  of  information  they  provide. 
It  appears  that  nearly  all  large  water  clusters  formed  in  nozzle  expansions 
have  the  cubic  structure  of  ice  Ic  (from  bulk  studies  [60],  this  implies 
an  effective  temperature  at  formation  of  about  100°-150°  K.).  However, 
only  rough  estimates  of  average  cluster  size  can  be  obtained,  and  it 
is  not  possible  to  track  the  growth  of  a  cluster. 

Our  data  provides  two  indirect  clues  about  the  behavior  of  large 
clusters.  First  the  monomer  velocity  distribution  undergoes  a  dramatic 
change  (becomes  bimodal)  when  a  certain  source  intensity  (size  of  cluster) 
is  reached.  At  present,  this  phenomenon  is  poorly  understood,  but  one 
can  speculate  that  some  significant  event  takes  place  when  some  large 
cluster  size  is  achieved.  Similarly,  since  the  slow  peaks  in  the  scatter¬ 
ing  data  appear  to  involve  the  scattering  of  large  clusters,  while  small 
clusters  simply  do  not  survive  collision  with  a  surface,  there  seem  to 
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De  two  transitions  in  the  scattering  behavior  as  a  function  of  cluster 
size.  Some  or  all  of  these  transitions  could  be  related  to  the  critical 
nucleus  size.  Clusters  larger  than  the  critical  nucleus  size  are  more 
stable  with  respect  to  gas  phase  collisions,  and  one  might  expect  them 
to  be  more  stable  with  respect  to  impact  as  well.  There  are  some  ex¬ 
citing  possibilities  for  studying  critical  nuclei,  if  these  speculations 
can  be  substantiated.  Some  means  must  be  found  for  monitoring  cluster 
size;  extending  the  range  of  mass  spectrometry  is  probably  still  the 
most  promising  approach. 

4.2  EVAPORATION  ANO  DETAILEO  BALANCING. 

The  usual  theoretical  argument  for  the  inference  that  the  velocity 
and  spatial  distribution  of  particles  evaporating  from  a  surface  is 
Maxwellian  is  as  follows;  Consider  a  solid  in  equilibrium  with  its  vapor. 
The  vapor  may  consist  of  clusters  as  well  as  monomers.  (Most  discussions 
of  evaporation  treat  the  vapor  phase  as  monomers.  At  least  one  other 
group  [61]  has  definitely  identified  small  clusters  in  evaporation  from 
metals.)  In  equilibrium  the  vapor  and  solid  are  at  the  same  temperature. 
The  vapor  has  Maxwellian  velocity  distributions  associated  with  each 
cluster  size.  The  surface  is  neither  growing  nor  shrinking,  so  detailed 
balancing  requires  that  any  particles  which  stick  to  the  surface  must 
be  exactly  balanced  by  particles  leaving  the  surface.  Furthermore  for 
equilibrium  to  be  maintained,  every  dynamic  process  must  be  accompanied 
by  its  inverse.  Monomers  and  clusters  In  the  gas  phase  have  a  Maxwellian 
distribution.  If  these  are  incident  on  the  surface,  they  must  leave  the 
surface  in  a  cosine  spatial  distribution  with  a  Maxwellian  velocity 
distribution  at  the  surface  temperature.  Since  the  incident  stream  is 
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independent  of  the  evaporating  stream,  if  the  vapor  is  removed,  and 
particles  no  longer  strike  the  surface,  the  evaporating  particles  should 
still  retain  the  equilibrium  properties.  How  then  can  one  explain  observed 
deviations  from  this  behavior? 

The  answer  is  that  the  detailed  balancing  argument  is  not  strictly 
valid.  In  equilibrium,  the  argument  is  sound,  but  the  extension  to  non¬ 
equilibrium  is  not  necessarily  valid.  For  example,  as  soon  as  collisions 
between  evaporating  molecules  become  likely  (when  the  stream  density 
becomes  high  enough),  equilibrium  arguments  fail.  If  there  are  gas 
phase  collisions  only  between  molecules  leaving  the  surface  and  not  be¬ 
tween  these  molecules  and  those  of  an  incoming  vapor,  there  is  no  reason 
to  expect  equilibrium  arguments  to  work,  and  a  specific  kinetic  model 
must  be  constructed.  This  problem  could  apply  to  all  the  data  for  higher 
target  temperatures  (Section  3.3,  Figure  22).  It  does  not  necessarily 
imply  that  the  higher  temperature  data  are  of  no  interest,  however.  Rapid 
non-equilibrium  evaporation  is  a  common  occurrence,  and  its  study  is  also 
interesting. 

4.3.  MONOMER  SCATTERING 

Me  consider  here  the  possibilities  for  the  contribution  of  monomers 
to  all  the  features  of  the  time-of-flight  scattering  data  (Section  3.4. 2,2). 
The  Maxwellian  peak  was  discussed  in  Section  3. 4. 2. 2.1.  A  plausible 
explanation  was  presented  involving  processes  of  trapping/ desorption.  We 
attribute  the  fast  peak  to  the  scattering  of  monomers,  but  we  must  also 
consider  possible  mechanisms  whereby  monomers  might  give  rise  to  the 
slow  peaks. 
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4.3.1.  THE  FAST  PEAK 


The  properties  of  the  fast  peak  were  presented  in  detail  in  Section 
3. 4. 2. 2. 2.  We  believe  this  scattering  feature  to  be  the  same  as  that 
reported  by  Hurst  et  al  [62]  for  heavy  inert  gas  scattering  from  .metals. 

They  used  only  grazing  incident  angles,  and  interpreted  the  scattering 
as  specular.  Our  observations  show  some  rather  unexpected  properties. 
Roughly  speaking,  the  scattering  appears  to  be  almost  elastic  in  velocity 
space,  but  inelastic  in  coordinate  space.  Figure  39  shews  incident  and 
scattered  velocity  components.  Note  that  neither  normal  nor  tangential 
momentum  is  conserved.  There  are,  in  fact,  regions  where  there  is  a 
gain  of  either  normal  or  tangential  momentum. .  The  total  scattered 
velocity  is  always  less  than  the  incident  velocity,  so  there  is  a  net 
loss  of  total  momentum.  There  is,  of  course,  no  problem  of  total  momentum 
conservation:  the  extra  momentum  is  either  gained  or  lost  by  the  surface. 
However,  most  simple  interaction  models  that  have  been  developed  to  explain 
classical  scattering  in  the  thermal  regime[63]  cannot  explain  any  gain  of 
momentum  components  (unless  the  surface  is  much  hotter  than  the  beam). 

Most  models  require  tangential  momentum  conservation  and  allow  for  either 
gain  or  loss  of  normal  momentum. 

A  simple  explanation  which  addresses  itself  to  all  the  features  of 
the  fast  peak  is  as  follows:  Referring  to  Figure  26,  note  that,  to 
within  experimental  accuracy  at  least,  the  (normalized)  amplitude  of  the 
fast  peak  is  constant  for  constant  "turning  angle"  (9>9$).  For  example, 
the  amplitude  at  es*75°  for  a  beam  Incident  at  9^*65°.  An  amplitude  pattern 
of  this  sort  could  be  generated  by  specular  scattering  from  a  rough  surface. 
All  surfaces  used  in  this  study  (see  Table  2)  can  be  expected  to  be  rough 


on  the  scale  of  a  single  water  molecule.  Even  an  optically  smootn  sur¬ 
face  need  only  be  smooth  down  to  a  scale  of  about  1  COO  X. 

For  any  given  incident  angle  (relative  to  the  average  surface  normal), 
a  full  range  of  local  surface  orientations  could  present  themsel  os  to 
an  incoming  molecule  from  the  beam.  Molecules  which  scatter  specularly 
relative  to  the  local  surface  normal  would  be  found  at  any  scattering 
angle  relative  to  the  average  surface  normal.  One  need  only  add  the 
hypothesis  that  molecules  which  scatter  specularly  do  so  as  a  function 
of  (local)  incident  angle,  and  a  complete  match  to  the  amplitude  data 

is  obtained.  A  Oebye-Waller  attenuation  may  be  adequate  to  formalize 

★ 

this  relationship. 

To  explain  the  small  loss  in  velocity,  then,  we  need  only  invoke 
a  small  energy  exchange  as,  for  example,  from  the  hard  cube  model.  Then 
normal  velocity/momentum  loss  manifests  itself  as  greater  velocity  loss 
small  turning  angles. 

The  data  are  not  sufficiently  accurate  to  warrant  detailed  numerical 

comparison.  There  are  a  couple  of  predictions  which  could  be  tested, 

for  which  data  have  not  yet  been  taken.  Scattering  outside  the  principal 

plane  should  also  show  the  fast  peak.  The  peak  amplitude  should  be 

constant  for  a  constant  angle  between  the  vectors  associated  with  the 

incident  and  scattering  direction.  (Note,  that  this  is  no  longer  equal  to 

the  sum  of  the  incident  and  scattering  angles.)  If  a  surface  can  be 

polished  smooth  to  the  scale  of  a  water  molecule,  real  specular  scattering 

should  occur.  This  may  not,  in  practice,  be  possible,  especially  when 

a  surface  layer  of  ice  is  present  (such  a  layer  would  probably  be  patchy 

*  A  study  of  this  possibility  is  under  way  with  Nicolas  Garcia,  but  the 
results  are  not  yet  ready  for  presentation. 


and  rough) . 


4.3.2.  THE  SLOW  PEAKS 

We  have  ascribed  the  slow  peaks  in  the  time-of-fl ight  spectra  to  the 
scattering  of  clusters  from  the  surface.  There  is  good  evidence  that 
this  hypothesis  is  reasonable  (see  Sections  3. 4. 2. 2. 3  and  3. 4. 2. 2. 4),  but 
without  positive  identification  of  clusters  in  the  slow  peaks,  one  should 
not  rule  out  the  possibility  of  monomer  scattering. 

The  scattering  cannot  involve  any  sort  of  thermal  equilibration: 

if  the  particles  are  monomers,  neither  the  target  nor  any  surface  in  the 

entire  system  is  cold  enough  to  produce  the  observed  velocities.  Hence 

for  monomers,  the  scattering  would  have  to  be  direct  scattering  with  an 

attendant  kinetic  energy  loss.  The  energy  can  either  be  lost  to  the 

surface  (via  surface  phonon  creation)  or  to  internal  modes  of  the  water 

molecule  (either  rotation  or  vibration).  The  observed  energy  loss  is 
-20 

about  2  x  10  J/molec.  There  is  a  strong  peak  in  the  bulk  phonon 

-20 

spectrum  of  ice  at  about  1.3  x  10  J/molec[2].  Of  the  possible  vi¬ 
brational  modes,  the  only  one  close  is  the  OH  stretching  mode  at  about 

-20  -20 
6.5  x  10  J/molec. [2]  Rotational  modes  are  excited  at  about  1.7  x  10 

J/molec. (2]  Both  translation-rotation  transfer  and  phonon  creation 
thus  appear  plausible.  Where  these  explanations  fail  most  seriously  is 
tithe  behavior  of  the  slow  peaks  as  a  function  of  nozzle  temperature 
for  constant  source  pressure.  There  is  a  sharp  cut-off  in  the  slow  peak 
amplitude  for  fairly  small  (<10%)  increases  in  temperature  (see  Section 
3. 4. 2. 2. 3).  This  is  reasonable  if  there  is  a  resonance,  but  a  correspond¬ 
ing  low  energy  cut-off  cannot  be  found.  An  explanation  in  terms  of  clusters 
appears  much  more  promising. 


4.4.  CLUSTER  SCATTERING 

Small  dusters  (less  than  75  molecules)  apparently  never  survive 
collision  with  a  surface  (see  Section  3. 4. 2. 4).  Above  some  critical 
size,  possibly  in  the  neighborhood  of  1  GOO  molecules,  it  is  plausible  to 
assume  that  some  fraction  of  the  clusters  survive,  at  least  in  part,  and 
give  rise  to  the  slow  peaks.  The  non-specular  slow  peak  apparently 
requires  even  larger  clusters  than  the  specular  slow  peak  (see  Sections 
3. 4. 2. 2. 3  and  3. 4. 2. 2. 4) 

Constructing  a  convincing  model  of  cluster  scattering  to  explain 
the  observed  features  is  a  formidable  task.  There  are  many  possible 
mechanisms  to  consider.  The  cluster  may  scatter  intact,  fragment,  or 
lose  a  few  surface  molecules.  It  may  change  structure  substantially  and/or 
have  internal  modes  excited.  Clusters  and  monomers  present  on  the  surface 
may  play  an  important  role.  One  can  imagine  a  sputtering  process  where 
the  incoming  cluster  knocks  material  loose  from  the  surface.  In  any  case, 
the  nature  of  the  surface  is,  at  present,  largely  unknown.  There  may 
be  a  layer  of  water/ice  one  or  more  molecules  thick.  Such  a  layer  might 
be  significantly  deformed  on  collision  with  a  cluster  even  if  sputtering 
does  not  take  place.  At  the  present  limited  state  of  knowledge  about  the 
system,  any  model  will  be  necessarily  somewhat  ad  hoc,  and  it  would  be 
difficult  to  decide  on  reasonable  approximations  to  create  a  numerical 
model.  There  are,  nevertheless ,  some  things  which  can  be  said  in  general 
terms  about  the  scattering. 

First,  it  should  be  noted  that  the  scattering  of  water  clusters  from 
surfaces  appears  to  be  quite  different  from  the  scattering  of  nitrogen, 
helium,  and  hydrogen  clusters  (see  Section  1.3) [44].  Neither  of  the  two 
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types  of  scattering  observed  with  these  other  species  has  features  similar 
to  those  of  the  slow  oeaks  observed  for  water.  The  model  proposed  to 
explain  one  of  the  scattering  patterns  for  hydrogen  and  helium  [44]  involved 
the  concept  of  evaporative  recoil:  some  molecules  from  the  surface  of 
the  cluster  evaporate  and  leave  the  surface;  conservation  of  momentum 
requires  that  the  cluster  gain  momentum  in  the  opposite  direction.  The 
problem  with  this  explanation  is  that  the  "evaporation"  from  the  cluster 
must  take  place  during  contact  with  the  target,  and  the  evaporating  mole¬ 
cules  probably  end  up  on  the  target  and  do  not  necessarily  carry  away 
any  momentum.  In  any  case,  such  a  model  does  not  predict  the  features 
observed  in  the  water  scattering  data. 

Consider  what  we  do  know  about  the  scattering:  The  incoming  clusters 
are  cold  (probably  50° -150°  K  )  and  moving  very  fast  for  their  size 
(nearly  as  fast  as  the  monomers).  We  need  to  explain  two  kinds  of  behavior. 
The  clusters  lose  much  of  their  initial  velocity  and  either  come  out  at 
the  specular  angle  or  near  the  tangent.  These  two  cases  are  considered 
separately. 

The  appearance  of  a  scattering  maximum  at  the  specular  angle  is 
startling  in  any  "classical"  scattering  system.  In  addition  we  have  just 
exolained  the  absence  of  specularly  directed  scattering  for  the  fast  peak 
Section  4.3.1)  by  postulating  a  rough  surface.  A  rough  surface  for 
v/w'  scattering  could  appear  smooth  for  cluster  scattering,  though, 

*  *n*  scale  of  roughness  were  in  between  the  size  of  a  monomer  and  the 
••  i  c’jst er.  We  have  also  considered  a  model  where  there  is  total 
'  in  the  locally  normal  direction  on  a  rough  surface.  Such 

^  r-.ws  ;-oduce  a  near-specular  reemission  angle,  but  it  does  not 
myway)  a  lobe  of  the  sort  observed,  and  does  not 


account  for  sufficient  loss  of  total  momentum. 


Returning  to  the  smooth  surface  (for  cluster  scattering)  hypothesis, 
we  have  considered  a  couple  of  other  models.  A  rough  analogy  to  the 
observations  is  the  bounce  of  a  very  soft  rubber  ball.  The  ball  is  so 
weakly  held  together  that  it  fragments  for  incident  angles  less  than  50°. 

For  larger  incident  angles,  it  survives  collision  but  loses  much  of  its 
translational  kinetic  energy.  The  loss  mechanism  might  be  conversion  to 
thermal  motion  via  vibrational  deformation  of  the  cluster  and/or  the  surface, 
or  perhaps  the  energy  could  go  into  breaking  some  of  the  bonds  of  the 
cluster  resulting  in  some  loss  of  mass.  Such  a  process  might  explain 
the  intensity  variations  with  incident  angle.  From  the  mass  spectro¬ 
metry  data,  it  appears  that  a  wide  range  of  cluster  sizes  must  contribute 
to  the  slow  peaks.  Perhaps  clusters  of  smaller  and  smaller  size  survive 
collision  with  the  surface  as  incident  angle  is  increased.  Detailed 
mass  analysis  with  a  mass  spectrometer  of  sufficient  range  might  verify 
such  a  behavior.  In  any  case,  the  equal  loss  of  normal  and  tangential 
momentum  must  be  accounted  for. 

The  non-specular  slow  peak  evidently  has  quite  a  different  scattering 
mechanism.  Almost  all  of  the  normal  momentum  of  the  incident  cluster  is 
lost.  A  large  fraction  (50-80%)  of  the  tangential  momentum  is  lost 
as  well.  A  rough  surface  explanation  is  of  no  use  (comforting  if  one 
believes  the  non-specular  slow  peak  to  involve  the  largest  clusters).  Again 
a  detailed  mechanism  can  be  only  guessed  at.  One  might  envision  a  cluster 
almost  sticking  and  then  rolling  or  bouncing/hopping  along  the  surface 
until  It  just  happens  to  get  enough  energy  to  escape  the  attractive 
potential  well  and  goes  skimming  off  near  the  surface.  It  is  certainly 
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remarkable  that  these  clusters  can  survive  collision  with  a  hot  surface. 
The  variation  of  threshold  temperatures  among  targets  is  quite  baffling. 
What  changes  can  be  taking  place  on  the  surface  at  threshold.  A  water 
monolayer  is  the  most  obvious  explanation,  but  what  role  does  it  play? 

The  immediate  directions  for  further  study  of  cluster  scattering 
include  the  following:  Mass  analysis  (or  more  precisely,  mass-filtered 
velocity  analysis)  of  the  slow  peaks  needs  to  be  accomplished.  This 
will  involve  a  very  high  mass  range  mass  spectrometer ,  but  as  long  as 
individual  mass  peaks  can  be  resolved  (and  for  time-of-fl ight  work  even  if 
they  cannot),  mass  resolution  can  be  sacrificed  for  range.  Some  other 
means  must  be  found  for  monitoring  the  condition  of  the  surface.  Helium 
scattering  or  LEED  would  probably  be  most  useful.  Auger  Electron  Spectro¬ 
scopy  or  a  quartz  crystal  microbalance  are  other  possibilities.  There 
is  probably  useful  information  to  be  gained  from  other  surfaces  as  well; 
we  still  do  not  know  what  surface  features  are  important. 


APPENDIX  A.  TIME-OF-FLIGHT  ANALYSIS  3Y  THE  PSEUDORANDOM  CHOPPING  TECHNIQUE: 

CALIBRATION  AMO  APPLICATION  TO  MASS  SPECTROMETRY 

This  appendix  contains  a  copy  of  a  paper  [53]  presented  at  the  11th 
Rarefied  Gas  Oynamics  Symposium  (July,  1973,  Cannes,  France).  It  was 
not  included  in  the  published  proceedings  (due  to  a  disagreement  over 
whether  it  contained  any  "new"  material),  and  is  therefore  reproduced 
in  full  here.  It  contains  its  own  figures  and  references. 
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Time-of-Fl ight  Analysis  by  the  Pseudorandom  Chopping  Technique: 

Calibration  and  Application  to  Mass  Spectrometry 

David  D.  Dreyfuss,*  Robert  B.  Doak,+  Harold  Y.  Wachman^ 

A  time-of-flight  apparatus  utilizing  a  cross  correlation  chopper 
for  beam  modulation  has  been  calibrated  against  a  Maxwellian  stream. 

A  consistent  and  reliable  scheme  for  determining  the  time  origin,  to 
±5  y  sec  for  any  arbitrary  time-of-flight  wave-form  has  been  devised. 
It  is  also  shown  that  the  need  for  deconvolution  over  shutter  function 
is  virtually  eliminated,  since  with  this  scheme  it  is  possible  to  work 
in  principle  with  an  arbitrarily  narrow  width  of  unit  slot.  A 
relatively  simple  and  straightforward  analytical  procedure  for  decon- 
voluting  the  modulated  signal  is  described. 
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1.  Introduction 

In  1966  we  reported  [1,  2]  from  this  laboratory,  on  the  calibration 
and  performance  of  a  time-of-fl ight  (TOP)  apparatus  (fig.  1)  in  which  the 
shutter  consisted  of  a  single  rotating  disk  having  four  narrow  slots  of 
equal  widths  located  at  90°  intervals  around  the  rim.  The  detector  con¬ 
sisted  of  an  electron  beam  ionizer  which  (ideally)  produced  a  sheet  of 
electrons,  which  ionized  the  molecular  stream  at  a  plane.  The  ions  were 
drawn  out  of  the  stream  and  collected  on  an  electron  multiplier. 

There  are  several  limitations  to  this  scheme.  The  time  resolution  of 
the  apparatus  is  determined  by  the  ratio  of  the  open  to  closed  arc  segments 
around  the  disk  (the  chopping  period,  which  is  the  sum  of  these  two  lengths, 
is  fixed  by  the  need  to  prevent  overlap  of  fast  molecules  from  one  pulse 
and  slow  molecules  of  the  previous  pulse).  To  obtain  good  time  resolution, 
one  uses  a  narrow  slot.  It  is,  in  principle,  possible  to  deconvolute  from 
the  effects  of  a  wide  slot  (shutter  function),  but  the  mathematical  manip¬ 
ulations  involved  are  very  sensitive  to  the  magnitude  of  the  noise  in  the 
data.  In  practice,  if  the  shutter  function  is  broad  enough  to  distort  the 
TOF  waveform,  then  the  best  that  deconvolution  provides  are  the  lowest 
order  moments  of  the  velocity  distribution  represented.  Usually,  in  the 
past,  it  was  found  necessary  to  sacrifice  information  on  velocity,  and  use 
a  wide  enough  slit  to  transmit  a  sensible  signal  in  each  pulse. 

The  ionization  detector  has  its  own  limitations.  An  electron  beam 
ionizer  is  a  universal  detector,  hence  is  usable  for  all  molecular  species 
In  the  test  stream.  As  a  consequence,  however,  background  molecules  are 
detected  as  well  (in  particular,  with  permanent  gases,  all  those  stream 
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molecules  which  strike  closed  portions  of  the  disk  are  detected  as  back¬ 
ground).  Statistical  fluctuations  in  this  background  are  usually  the 
major  source  of  noise  in  the  TOF  data.  Other  detector  problems  include 
the  finite  extent  of  the  ionization  region  along  the  beam  direction  and 
non-zero  ion  travel  time  between  the  instant  of  ionization  and  detection. 
The  extent  of  the  ionization  region  does  not  present  a  serious  problem 
because  it  has  only  a  second  order  effect  on  the  TOF  distribution  (the 
finite  shutter  function  is  a  first  order  effect).  We  have  shown  [3] 
that  a  considerable  flight  path  through  the  ionization  region  (up  to 
about  25%  of  the  total  flight  path)  can  be  used  before  appreciable  dis¬ 
tortion  in  waveform  occurs,  and  we  have  taken  advantage  of  this  to 
improve  the  efficiency  of  ionization.  As  to  the  effects  of  the  ion 
travel  time,  it  can  be  made  insignificant  with  sufficient  draw-out 
potentials  or  alternatively  it  can  be  accounted  for  in  interpreting  the 
data. 


2.  Pseudorandom  Chopping 

An  alternative  approach  to  molecular  beam  chopping,  which  avoids 
many  of  the  problems  cited  above  is  to  use  a  pseudorandom  chopping 
scheme  (fig.  2).  This  technique  was  originally  developed  in  connection 
with  thermal  neutron  studies  [4-7],  and  has  only  recently  been  used  with 
neutral  molecular  beams  [8,  9].  In  comparison  to  the  scheme  described 
above  a  largernumber  of  pulses  of  molecules  are  produced  over  each 
revolution  of  the  disk.  As  a  consequence  several  pulses  are  produced 
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within  each  chopping  period,  which  overlap  at  the  detector.  Sy  spacing 
the  pulses  "randomly"  within  the  period,  it  is  possible  by  a  simple 
arithmetic  procedure  to  undo  the  overlap  produced  [3,  9].  It  is  possible, 
in  fact,  to  attain  a  duty  cycle  of  5055  for  such  a  chopper  independent  of 
individual  slot  width.  With  this  chopper,  the  time  resolution  is  determined 
by  the  length  of  the  pseudorandom  chopping  sequence,  where  "length"  counts 
the  total  number  of  unit  slot  widths  (note  that  several  adjacent  locations 
may  be  open,  thus  forming  a  single  larger  slot  which  is  treated  as  a  set 
of  unit  slots).  The  deconvolution  procedure  recovers  the  waveform  which 
would  be  obtained  from  a  chopper  with  a  single  unit  slot  per  period.  While 
it  is  conceivable,  in  principle,  to  improve  time  resolution  arbitrarily  by 
taking  longer  and  longer  sequences,  in  practice,  limitations  in  electronic 
response  time  and  the  need  to  have  unit  slot  width  compatible  with 
molecular  beam  width  precludes  this  extension.  Avoiding  these  problems  by 
using  larger  diameter  disks  will  cause  other  experimental  difficulties. 
(Given  an  electronic  chopping  scheme  of  some  sort  (e.g.  chopping  laser 
excitation  of  a  beam),  this  last  limitation  may  not  apply,  since  chopping 
could  be  entirely  independent  of  collimation,  and  even  better  time 
resolution  might  be  obtainable.) 

3.  Calibration 

Although  the  pseudorandom  chopping  technique  has  been  demonstrated 
for  molecular  beams,  to  our  knowledge,  no  calibration  of  a  system  using 
the  technique  has  been  reported.  We  have  performed  such  a  calibration 
against  inert  gas  Maxwellian  beams.  A  schematic  diagram  of  the  apparatus 
Is  shown  in  figure  3.  A  room  temperature  gas  reservoir  at  a  pressure  on 
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the  order  of  10  Pa  supplies  a  beam  through  a  small  hole  in  a  thin  stainless 
steel  wall.  In  operation,  conditions  were  such  that  adequate  signals  were 
achieved  for  a  source  Knudsen  number,  Kn  ^  1.  Flow  at  this  value  of  Kn  is 
perhaps  not  fully  molecular  (one  would  like  to  use  Kn  =  100),  but  it  appears 
adequate  for  velocity  distribution  measurements  from  comparison  of  velocity 
distributions  for  a  series  of  decreasing  source  pressures  which  indicates 
no  measurable  change  in  distribution  when  values  of  Kn  ~  .1  are  used. 

Some  examples  of  calibration  runs  are  shown  in  figure  4.  To  compare 
the  data  with  the  known  velocity  distribution  of  the  gas  in  the  source, 
a  least  square  fit  of  a  Maxwellian  velocity  distribution  to  the  data  was 
made.  Source  parameters  such  as  molecular  weight  and  temperature  which  are 
known  or  independently  measurable  were  held  fixed.  Fitting  parameters  were 
a  baseline  position  (baseline  information  is  lost  experimentally  because  the 
signal  Is  superposed  on  a  large  DC  background),  and  an  overall  amplitude 
factor.  The  distribution  is  linear  in  these  quantities,  so  the  fitting  is 
straight- forward  and  well-defined.  It  avoids  the  hazards  of  peak  height 
matching  in  that  it  weights  all  data  points  equally,  and  can  be  used  for 
comparison  with  moderately  noisy  data.  Also  it  is  readily  extended  to  gas 
mixtures,  which  should  have  a  velocity  distribution  consisting  of  a 
weighted  sum  of  two  different  Maxwell ians. 

4.  Time  Origin  Considerations 

One  parameter  important  to  comparison  of  the  data  with  Maxwellian 
theory,  which  has  not  so  far  been  mentioned,  is  the  location  of  the  time 
origin  in  the  TOP  waveform.  Physically,  this  is  the  time  at  which  all 
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molecules  represented  by  the  TOF  waveform  passed  the  chopper.  Because  of 
the  sharp  leading  edge  of  the  Maxwellian  TOF  waveform,  small  differences  in 
time  origin  (less  than  one  channel  (slot)  width  or  15  usee)  are  readily 
detected.  In  making  a  calibration  it  is  especially  important  that  the  time 
origin  be  fixed  independently,  and  not  determined  by  matching  peak  locations 
or  the  leading  edge  or  by  some  similar  scheme. 

In  principle,  the  time  origin  is  easy  to  determine.  Using  a  photo¬ 
trigger  reference,  whose  position  relative  to  the  beam  position  is  measured, 
time  origin  can  be  determined  in  terms  of  some  constant  number  of  channels. 
Time'  delays  associated  with  either  ion  drawout  time  or  electron  signal 
processing,  are  essentially  constant,  and  may  also  be  assessed.  It  is 
difficult,  however,  to  measure  these  delays  independently  with  sufficient 
precision.  It  seems  best  to  measure  them  for  the  composite  apparatus.  Of 
the  several  schemes  attempted,  the  following  appeared  the  most  satisfactory. 

Using  the  complete  system  with  molecular  beam  and  all  processing 
electronics  in  place,  electronic  delays  were  determined  from  background  gas 
detection  without  the  chopper.  For  this  process  the  electron  accelerating 
grid  in  the  Ionizer  was  pulsed  electronically,  and  the  time  delay  from  the 
electronic  pulse  to  signal  rise  was  measured.  To  determine  geometric 
timing  relationships,  a  series  of  TOF  spectra  were  measured  for  a  constant 
Intense  signal.  (The  stream  was  not  necessarily  Maxwellian:  we  used  a 
fairly  high  pressure  HgO  beam  for  which  we  could  obtain  an  excellent  signal - 
to-nolse  ratio.)  These  TOF  waveforms  were  produced  by  rotating  the  chopper 
both  "forward"  and  "backwards"  (the  relative  location  of  the  photo-trigger 
and  molecular  beam  define  directions).  Then  the  relative  location  of  some 
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known  feature  (e.g.  peak  position)  for  the  two  resulting  waveforms  were 
compared  to  give  a  measure  of  the  geometric  timing  relationship  which  is 
independent  of  electronic  delays.  Improved  precision  in  the  measurement 
can  be  obtained  by  making  the  measurement  at  several  chopping  frequencies. 
We  have  been  able  to  locate  the  time  origin  to  within  about  ±1/3  channel 
(typically  ±  a  few  microseconds  for  our  operating  chopping  period  of 
1-2  msec  and  pseudorandom  sequence  of  length  103)  using  these  techniques. 

5.  Applications  to  Mass  Spectrometry 

One  possible  application  of  the  improved  resolution  TOF  system  we 
have  developed  is  to  mass  spectrometry  (see  fig.  5).  As  mentioned  above 
in  connection  with  the  fitting  of  theoretical  curves  to  the  data,  it  is 
possible  to  obtain  a  measure  of  the  relative  amounts  of  gases  in  a  gas 
mixture  by  comparing  amplitude  of  the  fitting  Maxwell ians.  The  advantage 
of  the  technique  is  that  it  is  non-destructive,  in  the  sense  that  ioniza¬ 
tion  takes  place  only  after  the  necessary  information  (time  of  flight 
between  two  fixed  points)  has  been  obtained.  There  are  difficulties, 
however,  because  the  velocity  distributions  for  the  species  to  be  separated 
must  be  known  independently.  This  information  is  available,  a  priori  only 
for  free  molecular  flows.  If  Maxwellian  velocity  distributions  must  be 
used,  only  a  few  fairly  widely  spaced  molecular  weights  can  reasonably  be 
separated  since  the  separate  TOF  peaks  are  broad  and  overlapping.  In  some 
experimental  situations,  for  example,  detection  of  dimers  and  higher  order 
clusters  in  near  condensing  systems,  the  advantage  of  non-destructive 
sensing  may  outweigh  the  disadvantages  cited. 
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FIGURE  1 

Typical  TOF  Apparatus 


FIGURE  2 

Pseudorandom  Chopper 
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FIGURE  3 
The  Apparatus 
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FIGURE  4 


calibration  Examples 
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APPENOIX  B.  DATA  PROCESSING  PROGRAMS  AND  DOCUMENTATION 


This  appendix  contains  user  instructions  for  the  data  handling 
software,  as  well  as  copies  of  the  programs.  Computer  hardware  is  described 
in  Section  2.6.  Programs  were  written  in  FORTRAN  IV.  Non-standard 
subroutines  include  a  collection  provided  by  the  manufacturer  for  support 
of  specific  hardware  (clock,  A/D,  and  graphic  display  terminal).  An 
additional  set  of  subroutines  were  specially  written  for  use  with  the 
plotter.  Documentation  for  these  routines  is  available  (and  is  kept 
with  the  system  software  documentation)  but  is  not  reproduced  here. 
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DOCUMENTATION  FOR  EXPERI. 

Includes  documentation  for  ell  routines  present  as  of  Matt  1980. 

Programs  and  documentation  written  by  David  Dreyfus*. 

The  program  is  accessed  btt  Putting  disk  labeled  EXPERIMENT  in 
drive  0  and  a  disk  for  data  in  drive  1.  Routines  are  available  to  collect 
and  store  data*  as  well  as  to  perform  a  variettt  of  data  processms  tasks. 

A  summary  and  list  are  available  btt  running  the  Program  (type  R  EXPERI) 
and  typing  HELP. 

All  routines  ars  accessible  in  any  order.  When  the  program  is 
started  and  whenever  a  -'utine  *s  finished*  a  prompt  character  (  >  ) 
indicates  that  a  new  routine  can  be  selected.  Simply  type  the  name  of  the 
routine  to  start  it.  Generally*  data  is  preserved  until  a  routine 
specifically  changes  or  replaces  it.  To  return  to  the  RT11  monitor*  type 
QUIT.  If  an  error  causes  the  program  to  abort*  it  should  be  restarted  using 
the  monitor  REENTER  command  which  will  give  the  EXPERI  prompt  character  <  >  )* 
and  leave  the  data  unchanged  (i.e.  if  a  spectrum  has  been  read  in  from  disk* 
it  will  still  be  there).  Many  errors  are  detected  and  dealt  with  by  the 
program.  Routines  with  two  word  names  can  be  written  as  either  one  or  two 
words  (with  or  without  a  space). 

Assumed  hardware  is  a  103  channel  pseudorandom  chopped  time-of-f light 
system  with  a  100  channel  Waveform  Eductor  providing  signal  averaging. 

1)  Oata  reading  and  file  maintenance. 

DATA  -  Reads  in  a  spectrum.  Queries  as  follows: 

NEW  or  OLD?  -  Select  NEW  to  accept  a  new  spectrum  via  the  A/D 

converter  and  store  it  on  disk*  or  OLD  to  retrieve  a  spectrum 
already  on  disk. 

File  name:  -  If  OLD  was  selected  this  is  where  the  file  name  is 

entered*  assumed  to  be  on  drive  1  and  have  a  .DAT  extension. 

Run  no.:  -  Spectra  are  labeled  by  the  run  number  entered  here.  A 

new  spectrum  is  stored  in  a  file  DXi:RUn.DAT*  where  n  is  the 
run  number.  This  puts  each  new  spectrum  in  a  separate  file 
to  provide  crash  resistance  for  the  data.  These  files  are 
generally  combined  into  one  file  at  the  end  of  a  data 
gathering  session  using  CONCAT.  Old  spectra  are  thus 
identified  by  file  name  and  run  number.  Old  spectra  from  the 
same  day's  run  would  be  accessed  by  file  name  RUn  and  run 
number  n. 

CONCAT  —  Concatenates  temporary  data  filee  into  one  larger  file.  Queries 

ae  follows: 

Output  file:  -  A  6  letter  (maximum)  name  for  the  output  file. 

Standard  practice  is  to  use  the  date  in  the  form  02MAY0* 
where  the  first  2  digits  are  the  day  of  the  month*  the  next 
3  characters  give  the  month*  and  the  final  disit  is  the  final 
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digit  of  th*  stir  (e.g.  0  for  199C-- 

Are  films  consecutively  numbered  starting  with  l7  -  Answer  YES  or  NO 

as  appropriate  Usually  th*  answer  should  b*  YES. 

How  mans  files?  -  Give  nuabsr  of  files  to  b*  concatenated. 

If  th*  filas  arm  consecutively  njr  era d»  the  concatenation  proceeds. 

If  not.  each  file  must  ba  iivan  i  .hen  asked  for)  by  name.  The  process 
•aw  be  slow,  since  a  lot  of  disk,  operations  arm  reauired. 

2)  Plotting  routines. 

Many  of  the  data  erocassins  routines  automatically  produce  Plots,  but 
in  addition  there  are  the  following  routines  for  controllins  Plots! 

PLOTON 

PLOTOFF  - —  These  routines  simply  enable  and  disable  Plotting  on  the  HP  72C3A 
Platter.  The  plotting  is  sufficiently  sl.u  that  standard  practice  is 
to  use  the  plotter  only  to  Produce  report  Quality  reproductions.  For 
routine  Plotting  the  display  screen  is  adeauate.  and  the  terminal 
copier  can  provide  Quick  hard  copies.  The  default  condition  when  the 
program  is  started  is  with  the  plotter  disabled. 

ERASE 

DISPLAY  -  These  routines  can  suppress  or  reactivate  th*  graphic  display  on 

th*  terminal  (primarily  to  avoid  confusion  with  text).  Stored  data 
is  not  affected,  and  subseauent  Plotting  proceeds  regardless  of  how 
these  routines  are  used. 

HISTOGRAM  -  Allows  activation  or  suppression  of  histogram  display  on 

terminal  (filling  in  space  below  graph).  Queries  are  obvious. 

PLQTBATA  -  Plots  raw  data  (prior  to  deconvolution).  Generally  used  to  check 

visually  that  data  was  read  in  correctly,  and  to  check  for  necessary 
interpolation. 

3)  Data  manipulation  routines. 

These  routines  change  a  spectrum  in  various  ways.  Th*  raw  data  is 
always  maintained  on  disk.  In  addition  th*  current  spectrum  can 
also  be  recovered  by  reinvoking  BECQNV  (it  is  not  necessary  to  reread 
the  raw  data  unless  a  smaller  interpolation  is  desired). 

DECONU  -  Recovers  a  time-of-f 1 isht  spectrum  from  th*  rseudo random-chopped 

raw  data.  Queries  as  follows: 

FQRUARD  or  REVERSE  chopping'7  -  Normally  th*  chopper  is  run  in  one 

direction  labeled  ’FORWARD*.  If  for  some  reason  th*  chopper 
is  run  in  th*  other  direction.  REVERSE  can  be  selected. 

The  channels  around  th*  region  to  be  interpolated  are  printed,  and 
then  th*  muerw! 

Low  and  high  channels  for  interpolation!  -  If  all  100  channels  of 

th*  Waveform  Eductor  signal  are  good,  th*  standard  answer 
is  100. 10A.  It  is  sometimes  desirable  to  extend  the  interpo¬ 
lation  region  by  on*  or  two  channels  in  either  direction. 
Interpolation  is  simply  linear  between  the  endpoints. 

A  plot  of  th*  deconvoluted  spectrum  is  automatically  produced. 


SUBMAX 

SUBTRANS  -  Subtract  a  manually  chosen  Maxwellian  (SUBMAX)  or  translating 

Maxwellian  (SUBTRANS)  fro*  a  deconvoluted  spectrum.  This  procedure 
is  invoked  when  a  known  Maxwellian  component  is  obscuring  other 
features  of  interest.  Queries  as  follows: 
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Flight  path  (cm.): 

Period  (msec.): 

Eductor  delay  (usec.>: 

Temperature  (centigrade) : 

Mach  number:  C SUB TRANS 

Molecular  weight: 

Maxwellian  amplitude: 

Eductor  delay  is  the  delay  fro*  the  center  of  t ha  trigger  rulse  to  the 
baainnina  of  the  first  channel  as  sat  on  the  Uaveform  Eductor. 
Maxwellian  amplitude  is  the  peak  amplitude. 

Aftar  tha  ravisad  spectrum  is  plotted  comes  the  auers! 

Another  try?  -  A  response  of  NO  or  <CR>  ends  tha  routine.  A 

response  of  YES  allows  tha  selection  of  a  new  amplitude  in 
case  the  subtraction  does  not  appear  reasonable.  The  auery 
•Maxwellian  amplitude:*  is  repeated*  and  a  new  plot  is 
produced.  The  cycle  can  be  repeated  indefinitely. 

SMOOTH  — -  smooths  data  bw  low-pass  digital  filterin*  (sa#  Kaisar  and  Reed* 
R.S.I.  48:1447  (1977)  for  datails  of  tha  tachniaue).  Lambda  is 
set  to  80  dB»  and  delta  is  sat  to  about  2.5  Fourier  components*  which 
sets  NP  to  103*  divine  good  accuracy  and  reasonable  run-time.  Queries 
are  as  follows: 

How  many  Fourier  components  to  be  kept?  -  For  typical  time-of-f light 

spectra*  around  10-15  components  represent  tha  usable 
information.  For  a  103  channel  sianal*  52  components  contains 
all  the  information  available.  This  filtering  process  does 
affectively  eliminate  all  of  the  hiah-f reouency  noise*  but 
since  tha  low  freouencv  noise  is  generally  of  about  the  same 
amplitude*  eliminating  Just  the  high  freouency  noise  tends  to 
create  features  which  aren't  really  present.  The  human  brain 
seems  capable  of  doing  a  more  realistic  smoothins  operation 
on  the  unsmoothed  data  by  fittind  a  curve  of  approximately 
the  right  shape. 

The  smoothed  data  is  Plotted  over  the  oriainal  spectrum.  Additional 
Queries  allow  another  choice  of  number  of  components  and 
ask  whether  the  smoothed  version  should  replace  the  oriainal. 

4)  Information  extracting  and  fitting  routines. 

PEAKS  -  Locates  peak  positions  and  amplitudes  and  calculates  average 

velocity  for  any  time-of-f light  spectrum.  Queries  are  as  follows: 

Flight  path  (cm.):  ! 

Period  (msec.):  Responses  as  for  SUBMAX 

Eductor  delay  (usee.)!  ! 

How  many  channels  are  "zero*?  -  Humber  of  channels  at  beginning 

of  spectrum  which  can  be  taken  as  zero  for  determining  a 
baseline  (enter  as  integer  Cno  decimal  roint3). 

Peaks  found  are  marked  on  display*  and  positions  and  amplitudes  are 
printed  along  with  the  average  velocity. 

MAXFIT  -  Fits  one  or  two  Maxwellians  of  the  form  U**2*EXP(-V**2>  to  a 

spectrum.  Plots  fit.  Queries  are  as  follows: 

How  many  species?  — -  1  or  2*  Determines  whether  1  or  2  Maxwellians 
are  to  be  fit  (integer). 

Flight  path  (ca.>:  ! 

Period  (msec.): 

Temperature  (centigrade):  Responses  as  for  SUBMAX  and  PEAKS 

Molecular  weimhtCsl: 


Responses  are  obvious.  All 
numbers  should  be  entered  as 
real  (with  decimal  point). 
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How  many  channel*  ere  ,:»ro"' 

Uhat  fraction  do  you  want  to  fit  (starting  from  left)''  ---  Real  number 
between  0.  and  1.  indicating  how  much  of  the  spectrum  is  to  be 
used  for  fitting  (allows  fitting  to  a  tail  or  fitting  to  one 
eortion  of  spectrum  if  desired) •  This  option  is  available 
only  for  a  single  Maxwellian  fit. 

Fitting  is  by  least  sauares.  Fitting  parameters  are  just  the  ampli¬ 
tudes  of  the  1  or  2  Maxwellians  specified.  The  two  Maxwellian*  are 
distinguished  only  by  their  molecular  weight  as  might  be  appropriate 
for  a  gas  mixture. 

TRANSFIT  — -  Similar  to  MAXFITi  but  one  of  the  Maxwellian*  can  be  a 

translating  Maxwellian:  V#S2*£XP(-(V-U)**2> .  Queries  as  follows! 

How  many  species?  -  If  1  is  selected  a  fit  is  made  for  a  trans¬ 

lating  Maxwellian.  2  gives  on*  of  each  type. 

Other  Queries  as  for  MAXFIT  with  the  following  changes!  Only  one 
molecular  weight  can  be  specified.  Two  temperatures  are  allowed’  and. 
of  course>  a  Mach  number  is  needed  for  the  translating  Maxwellian. 

The  Mach  number  is  the  ratio  of  U  to  tha  average  thermal  velocity! 

SORT (8kT/<pi )*) . 

SAVE  -  Stores  a  spectrum  for  use  by  DOUBLEFIT. 

DOUBLEFIT  -  Fits  the  sum  a  a  spectrum  stored  by  SAVE  and  a  Maxwellian. 

Queries  as  for  MAXFIT!  with  a  warning  that  SAVE  must  be  invoked 
firsti  and  a  Question  as  to  whether  a  spectrum  has  been  stored. 
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PROGRAM  EXPERI 

Monitor  rrasra*  for  control 


of  data  acouisition  and  rrocassins 


REAL  DECONWf 309) .MAXI .MINI 
INTEGER  DA  TA  <  206 ) » GRAPH (312) 

BYTE  COMMAN <21 ) »T0ATE<9) » TTIME<  8 ) >RUNNG<3) .ERR 

The  folloMin*  values  are  passed  around  aeons  subroutines. 

COMMON  /RAUD/  DATA  .MIN .MAX »  /DECD/  DECGNV. MAXI .MINI , ICHNL1  .GTO 

+./PL0T/  GRAPH.  /RIJMUM/  RUNNO 

CALL  DATE  <  TOATE ) 

CALL  TIME^TTIME) 

T*10*<TTIME< 1)-48)+TTIME(2)-4S 
IF  (T  ,LE.  12)  GOTO  3 
IT-(T-12)/10 


TTIME<1)*48+IT 

TTIHE(2>»36+T-IT*10 

IF  (TTIME(l)  .EQ.  48)  TTIME<1)*32 

CALL  PLTOFF 

CALL  PL0T33< 13 » 72 . ) 

CALL  FLOTSSdO..) 

TYPE  10.TDATE, (TTIME(I) .1*1.3) 

FORMAT ( '  Data  acouisition  i  rrocesains  s*st»*  version  i: 
+  last  chansed  4-MAR-80.V'  Todaa  i*  '»9A1»'  and  it  is  ' 
+.3A1.'  o' 'clock..  V/) 

Entar  rroart  character. 

TYPE  30 

FORMAT  ('  >'.»> 

CALL  GETSTR (3. COMM AN >20) 

CALL  SCOPY ( COMMAN . COMMAN .16) 

CALL  TRIM (COMMAN) 

Identify  coaaand  and  invoke  appropriate  routine. 

CALL  SCOMP  < ' HELP ' »  COMMAN . I VAL ) 

IF  (IWAL  .NE.  0)  GOTO  40 
CALL  HELP 
GOTO  20 

CALL  SCOMP < 'DATA' .COMMAN. IWAL) 

IF  ( IWAL  .NE.  0)  GOTO  30 
CALL  READER 
GOTO  20 

CALL  SCOMP ( 'DECONW' .COMMAN . I VAL > 

IF  (IWAL  .NE.  0)  GOTO  60 
CALL  DEKONW 
GOTO  20 

CALL  SCOMP < 'MAXFIT' .COMMAN. IWAL) 

IF  (IWAL  .NE.  0)  GOTO  70 
CALL  MAXFIT 
GOTO  20 

CALL  SC0MP< 'PLOTDATA' .COMMAN. IWAL) 

CALL  SC0MP('PL0T  DATA' .COMMAN. JVAL) 

IF  ( IWAL  .NE.  0  .AMD.  JWAL  .NE.  0)  GOTO  80 
CALL  PLDATA 
GOTO  20 

CALL  SCOMP < 'CONCAT' .COMMAN. IWAL) 

IF  (IWAL  .NE.  0)  GOTO  90 
CALL  JOIN 
GOTO  20 

CALL  SCOMP  < ' ERASE ' . COMMAN . I VAL ) 

IF  (IWAL  .NE.  0>  GOTO  100 
CALL  ERASE 
GOTO  20 

CALL  SCOMP ( 'DISPLAY' .COMMAN. IWAL) 

IF  (IWAL  .NE.  0)  GOTO  110 
CALL  CSPLAY 
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GOTO  20 

UO  CAUL  SCOMP(  'HISTOGRAM'  .COMMAN. IVAL) 

IF  ( I UAL  .NE.  0)  GOTO  120 
CALL.  HISTO 
GOTO  20 

120  CALL  SCOMP( 'OUIT' .COMMAN, IVAL) 

IF  (IVAL  . NE .  0)  GOTO  130 
CALL  EXIT 

130  CALL  SCOMP (  'SAVE'  .COWMAN.  IVAL) 

IF  ( I UAL  . NE.  0)  GOTO  140 
CALL  SAVE 
GOTO  20 

140  CALL  SCOMP( 'DOUBLE  FIT ', COMMAN, IVAL > 

CALL  SCOMP( ' DOUBLEFIT ' » COMM AN »  JVAL ) 

IF  ( IOAL  .NE.  0  .AND.  .JVAL  .NE.  0)  GOTO  ISO 
CALL  DBLFIT 
GOTO  20 

130  CALL  SCOMP ( ' TRANSFIT ' .COMMAN , I UAL) 

IF  < IVAL  .NE.  0)  GOTO  160 
CALL  TRAP IT 
GOTO  20 

160  CALL  SCOMP ( 'SUBMAX' .COMMAN, IVAL) 

IF  (IVAL  .NE.  0)  GOTO  170 
CALL  SUBMAXI  1) 

GOTO  20 

170  CALL  SCOMP  (  'SU8TRANS  '  .COMMAN*  I  UAL) 

IF  ( IVAL  .NE.  0)  GOTO  130 
CALL  SUBMAX (2) 

GOTO  20 

130  CALL  SCOMP  < ' PEAKS ' ,  COMMAN . IVAL ) 

IF  (IVAL  .NE.  0)  GOTO  190 
CALL  PEAKS 
GOTO  20 

190  CALL  SCOMP ( 'SMOOTH' * COMMAN, IVAL) 

IF  < IVAL  .NE.  0)  GOTO  200 
CALL  SMOOTH 
GOTO  20 

200  CALL  SCOMP ( 'PLOTON ' * COMMAN  * IVAL ) 

CALL  SCOMP ('PLGT  ON ' * COMMAN  *  JVAL ) 

IF  < IVAL  .NE.  0  .AND.  JVAL  .NE.  0)  GOTO  210 
CALL  PLTON 
GOTO  20 

210  CALL  SCOMP ( 'PLOTOFF ' .COMMAN , IVAL > 

CALL  SCOMP ( ' PLOT  OFF ' , COMMAN  * JVAL ) 

IF  (IVAL  .NE.  0  .AND.  JVAL  .NE.  0)  GOTO  220 
CALL  PLTOFF 
GOTO  20 

220  CONTINUE 

TYPE  510 

510  FORMAT ( '+Whet*' ) 

GOTO  20 

C  TO  ADO  A  NEW  ROUTINE*.  WRITE  ROUTINE*  ADD  4  LINES  AS  ABOVE ,  AND 

C  ADD  NEW  INSTRUCTION  TO  LIST  IN  HELP  SUBROUTINE. 

END 

SUBROUTINE  REPL( INSTRN. RNUM) 

C  Extends  capabilities  of  PLTRER  to  handle  larder  numbers. 

BYTE  INSTRN(Sl) 

IF  .RNUM  .LT.  10000.)  GOTO  10 
CALL  CONCAT < INSTRN. '  ♦♦'  » INSTRN) 

CALL  PLTREK INSTRN* INT(RNUM/10000. ) > 

RNUW«RNUM— l 0000 . * INT ( RNUM/ 10000 .  ) 

IF  (RNUM  .GC.  1000)  GOTO  10 

CALL  INSERT < ' OS* '. INSTRN, INDEX ( INSTRN* '*#')) 

IF  (RNUM  .QE.  100)  GOTO  10 

CALL  INSERT < '0**' .INSTRN* INDEX (INSTRN. '*♦' ) > 
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IF  (RNUM  .GE.  10)  GOTO  10 

CAUL  INSERT ( '0 ♦  »'  . INSTRN. INDEX ( INSTRN.  '**')> 

10  CALL  PUTREK INSTRN. IMKRNUM+. 3) > 

CALL  CONCATt INSTRN. ' . ' .INSTRN) 

RETURN 

END 

subroutine  help 

TYPE  10 

10  FORMAT  ('  Current  valid  command*  are:'/ 

+'  HELP ' . TOO .' P roduce*  this  list.'/ 

+  '  DATA '» T20 »' Reads  in  a  tieu  or  old  spectrum  and  stores  it.'/ 

+'  PLOTDATA ' . TOO . 'P lo ts  ran  data.'/ 

+'  DECONU ' . T20 . ' Deconvolu tes  data  and  plots  it.'/ 

+  '  SUBMAX '» T20 »' Subt racts  a  manually  chosen  Maxwellian  from  data  1 
+  reelots . ' / 

+  '  SUBTRANS' »T20» 'Same  as  SUBMAX.  but  using  translating 
+  Maxwellian.  '  / 

+  '  SMOOTH '» T20 >' Per f a rms  low-pass  filtering  on  deconvoluted  data.'/ 
+  '  MAXFIT ' . T20 . 'Fi ts  one  or  two  Maxwei 1 lans  to  deconvoluted  data  t 
+  plots' /T23 .' over  data  Plot.'/ 

+  '  SAVE' »T20. 'Stores  a  deconvoluted  trace  for  DOUBLE  FIT.'/ 

+  '  DOUBLE  FIT' .T20» 'Fits  the  sum  of  a  Maxwellian  and  a  stored 
+  waveform'/T23» 'and  plots  over  data  Plot.'/ 

+'  TRANSFIT ' . T20 » ' Fi ts  a  translating  Maxwellian  and  plots  it  or 
+  fits  the'/T23» 'sum  of  that  and  a  regular  Maxwellian  and 
+  plots  it. '/ 

+'  PEAKS' »T20» 'Locates  and  gives  resolved  peak  amplitudes;  finds 
+  averase'/T23» 'velocity. '/ 

+'  ERASE' »T20. 'Stops  all  graphic  display.'/ 

>'  DISPLAY' »T20» 'Reactivates  graphic  display.'/ 

+'  HISTOGRAM' »T20. 'Disables  or  re-enables  histogram.'/ 

+'  C0NCAT' »T20» 'Condenses  all  temporary  data  files  into 
+  one  f ile . ' / 

+'  PLOT  ON' .T20. 'Activates  7223A  plotter.'/ 

+'  PLOT  OFF ' .T20 .' Deactivates  7223A  Plotter.'/ 

+'  QUIT' »T20» 'Returns  system  to  RT-ll  monitor.'/ 

+'  For  more  help,  type  ‘TYPE  SYJEXPERI.DOC*  from  RT-ll  monitor.') 

RETURN 

END 

SUBROUTINE  ERASE 
CALL  PL0T33(2.0.1) 

RETURN 

END 

SUBROUTINE  OSPLAY 
CALL  PL0T33<  2 .1.0) 

CALL  PL0T33< 13.72. ) 

CALL  PLQT33< 10» » ) 

RETURN 

END 

SUBROUTINE  PLDATA 

INTEGER  0ATA(2OA) »GRAPH(312) 

COMMON  /RAUD/  DATA. MIN. MAX.  /PLOT/  GRAPH 
ICHNL1-60  (set  starting  point  for  Plot. 

CALL  PL0T33  <2.312. ) 

CALL  PLT33U.0.) 

CALL  PLT33  <2.137. 374 ) 

CALL  PLQT33( 13.72. ) 

CAU.  PLQT33(10. ,  1 
DO  10. IX*0 .303*3 

XY»INT (FLOAT (DATA < ICHNL1 +IX/3)-MlN) *233. / ( MAX-MIN)  ) 

DO  10.  1-1.3 
10  GRAPH( XX+X )-IY 
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I Y*INT( FLOAT (DATA (  I CHNL1  r  1 02  )  — M  IN  )  *233  . / ( MAX-M IN  >  > 

GRAPH  ( 31 1 ) *  I Y 
GRAPH<312)*IY 
CALL  PL T33<  7,0,0) 

CALL  PLT33 (3.-312,  GRAPH ) 

CALL  PLT33 (2*1 . ) 

RETURN 

END 

SUBROUTINE  HI3TC 
BYTE  ANSWER ( 10 ) *  ERR 
ERR*. FALSE. 

3  TYPE  10 

10  FORMAT  ('♦Turn  ON  or  OFF'*  ',*> 

CALL  GETSTR ( 3 »  ANSWER  >  9  , ERR ) 

CALL  SCOMP ( 'ON' .ANSWER, IVAL) 

IF  (IVAL  .EO.  0)  GOTO  30 
CALL  SCOMP ( ' OFF ' , ANSWER  » IVAL ) 

IF  (IVAL  .EQ.  0)  GOTO  70 
TYPE  20 

20  FORMAT  (  '♦What?'/) 

GOTO  3 

30  CALL  PL0T33  (2,8,0) 

GOTO  100 

70  CALL  PL0T33 ( 2  »0  » 3  ) 

100  CALL  PL0T33( 13,72,  ) 

CAL_  PL0T33 ( 10  ,  ,  ) 

RETURN 
*  END 

SUBROUTINE  DEKONV 
BYTE  TDATE<  81 ) .RUNNO ( 3 ) 

COMMON  /RUNUM/  RUNNO 
REAL  MAXI. MINI. 0EC0NV(309> 

INTEGER  DATA <  206  > ,HIGH.GRAPH(312) 

BYTE  S( 103) , IDIREC IS) »ANSWER< 12) , ERR . IVAL , JUAL 
COMMON  /RAUD/DATA, /DECD/OECCNV  »MAX1 .MINi.XCHNLl.GTO 
+  ,  /PL OT/  GRAPH 

C  Sat  up  dac  (involution  saauanca,  sat  IDIF:EC  initially  to  'FORWARD'. 

DATA  ERR/ .FALSE ./, 3/7* 1,3* -1,3* 1,-1, -1, 1,-1, 3*1 ,-1,-1, 1,-1, 1,1, 

+3*—  1 , 3*1 ,7*-l ,3*1 »-l » 1 »-l »-l » 1 >  3*”1 , 1 »3*-l , 1 ,-l , 1 »-l > 1 » 1 ,- 1 » 

♦4*1, -1,1, 1,-1, -1,1, -1,1, 1,-1, -1,1» 4*— 1,1,— 1,— 1,1,— 1,1,— 1,3*1, 

♦- 1, 3* 1,-1. 1,1, -1,1, 3* -1/. IDIREC/70 ,79,32,87,63,82,68,0/ 

3  TYPE  10 

10  FORMAT  ('+FORWAR0  or  REVERSE  chopping^  - ,s) 

CALL  GETSTR (3, ANSWER, 10, ERR) 

CALL  SCOMP ( 'FORWARD' .ANSWER, IVAL) 

IF  (IVAL  . NE.  0)  GOTO  20 

ICHNL1*3?  list  channel  to  Plot. 

GT0«60.9  (Geoaetric  tin*  origin  (no  dela*  corrections )  . 

GOTO  30 

20  CALL  SCOMP ( ' REVERSE ' , ANSWER , IVAL  > 

IF  (IVAL  .EO.  0)  GOTO  40 
TYPE  30 

30  FORMAT  ( '♦What*'/) 

GOTO  3 

40  ICHNL1*40  (1st  channel  to  Plot 

GT0*42.1  (Geoaetric  tin*  origin. 

C  Reverse  deconvolution  saouanc*  if  direction  is  changed. 

30  IF  (ISCOMPdDIREC .ANSWER)  .EO.  0)  GOTO  70 

CALL  SCOPY(ANSWER.IDIREC) 

DO  60,1*2,32 
J*S(I) 

S( I >*S( 103-1 ) 

60  3( 103-1 )®J 

C  Choosa  channals  for  interpolation  and  perf or»  linaar  interpolation. 
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non 


70  type  ao 

30  FORMAT  <  '+CHannels  96  thru  103'  ) 

TYPE  90.  (DATA( I ) . I  -96. 103) 

90  FORMAT  (1313) 

TYPE  100 

100  FORMAT  ('  Low  and  hiah  channels  for  interpolation:  '.»> 

ACCEPT  110. LOU. HIGH 
110  FORMAT  (214) 

00  140. I=L0U+1 .HIGH-1 

0ATA(  I  >  =  ( I-LOU  )  *  (  DAT  A  ( HIGH  )  -DAT.4<  LOU  ) ) / ( HIGH-L3U ) +DATA < LOU ) 

IF  (1-103)  130.100.120 
120  DATA( I-103)=OATA( I ) 

GOTO  140 

130  DATA(I+103)=OA7A(I) 

140  CONTINUE 

C  Deconvolute1 

DO  130  1=1.103 
130  DECONV  < I ) =0  . 

MAX1=-10000. 

MIN1  =  10000 . 

DO  170.1=1.103 
DO  160  »  J=1 » 103 

160  0EC0NV( I)=D£C0NV< I >+S( J>*FLOAT(DATA( I+J) ) 

DECONV ( I ) =-D£C0NV ( I ) 

DECONV(I+103)=OECONV(  I) 

DECONV ( I +206 ) =D£C0NV  < I ) 

MAX1*AMAX1 ( MAXI » DEC0NV(  I ) >  1  Deter® me  max  and  min  values  for  plot 

170  MIN1=AMIN1 (MINI • DECONV  til) 

PLOTTING  ROUTINE 

CALL  PL0T33U.0.) 

CALL  PL0T33 (2.512. ) 

CALL  PLT33 (2.9.302) 

CALL  PL0T33( 13.72.  ) 

CALL  PLQT33 ( 10 . .  ) 

DO  200.  IX=0 .303 . 3 

IY=INT( (DEC0NU( ICHNL1+IX/3 ) -MINI >*201  ./ C MAXI -MINI > >+17 
DO  200.  1=1.3 
200  GRAPH ( IX+I)=IY 

IY=INT ( ( DECONV ( ICHNU +102) -MINI) *201 ./(MAXI-MINI) )+17 
GRAPH(3ll  )  =  1Y 
GRAPH(313>=IY 
CALL  PLT33 (7.0.0) 

CALL  PLT33( 3 .-312 » GRAPH) 

CALL  PLT33 (9.44.0) 

CALL  DATE ( TDATE ) 

TDATE ( 10 ) =0 
CALL  PLTTXT <  TDATE ) 

CALL  PLT33 (9.44.1) 

CALL  CONCAT ( 'Run  *' .RUNNO . TDATE ) 

CALL  PLTTXT (TDATE) 

CALL  PLT33 (9.44.2) 

CALL  SCOPY ( 'Full  scale  amplitude:  ♦*'» TDATE) 

CALL  REPL (TDATE* MAX 1-M INI ) 

CALL  PLTTXT( TDATE) 

CALL  PL0T33( 13.72.) 

CALL  PLT33(2.1.) 

RETURN 

END 

SUBROUTINE  MAXFIT 
BYTE  TDATE (81 > .RUNN0(3) 

COMMON  /RUNtJM/  RUNNO 

REAL  MAXUEL. DECONV (309) .MAXI. MINI 

MAXUEL ( V  >  =EXP ( - ( V*V ) ) *V**4 
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C  This  is  a  static  Maxwellian.  The  'J** 4  occurs  because  the  riots  art 

C  tiae-of-f lisht  Cdv»L/t**2  dt*v«*2/L  dt>. 

INTEGER  GRAPH(312> 

COMMON  /DECO/  DECONV. MAXI .MINI . ICHNL1 .3T0 
3  TYPE  10 

10  FORMAT  ('-t-How  many  sxcus7  '»*> 

ACCEPT  20. NSPEC 
20  FORMAT  (13 > 

IF  (NSPEC  .£0.  1  .OR.  NSPEC  .EO.  2)  GOTO  40 
TYPE  30 

30  FORMAT  ('+1  or  2  only!'/) 

GOTO  5 
40  Z£R0*0. 

TYPE  30 

SO  FORMAT  (  '  +F1 lSht  rath  (ca.>:  '.S> 

ACCEPT  60 .PATH 
60  FORMAT  (F9.4) 

TYPE  70 

70  FORMAT  ('(-Period  (i»sec.)5  '»*) 

ACCEPT  60. PERIOD 
TYPE  80 

80  FORMAT  ('+Eductor  delay  (usac.)i  '.S) 

ACCEPT  60. DELAY 

TO*GTO-(OELAY+3. X.103/PERIQD  'locates  tin>e  onain.  3  is  for 
PERI0D-PERI0D/1000.  (electronic  delay. 

IT0*INT(T0+1 . )  '1st  channel  for  fittins. 

IZER0*INT(3.*(T0-ICHNL1 ) )  !No.  of  roints  at  besmnins  of  fitted 

IZ*IZER0  ! curve  to  set  to  zero. 

IF  ( I ZERO  .LT.  0)  IZERO-O 
TYPE  90 

90  FORMAT  < '  +Te«reraturs  ( centisrade ) :  '»*> 

ACCEPT  60. TEMP 
TEMP*TEMP+273.13 
GOTO  (100.120)  NSPEC 
100  TYPE  110 

110  FORMAT  ('+Molecular  ueiaht:  '.S) 

ACCEPT  60. WEIGHT 

TSCALE*SQRT(UEIGHT/TEMP*6.016S9E-9)*PATH  Iscalins  tiise  for 

GOTO  130  ! MAXUEL  velocity. 

120  TYPE  130 

130  FORMAT  ('+Molecular  weUhtsJ  '.$) 

ACCEPT  140.UEIGH1 .UEIGH2 
140  FORMAT  (2F9.4) 

TSCAL1«SQRT( WEIGH1/TEMPS6 . 0 1 639E-9 ) *PATH 
TSCAL2-SQRT ( UEIGH2/TEMP46 . 0 1 639E-9 ) *PATH 
130  TYPE  160 

C  Set  baseline. 

160  FORMAT  ('+How  isany  channels  are  ‘zero’?  '.»> 

ACCEPT  20.NZER0 

DO  170.  I*IT0-NZ£R0 » ITO+NZERO-l 
170  ZER0»ZER0+DEC0NO(I> 

C  Set  rrefittins  roints  on  fitted  curve. 

ZER0-ZER0/NZER0/2 

IY*INT( (ZERO-MINI ) 4201 ./(MAX1-MIN1>)+17 
DO  180.  1*0 . I2ER0+3 
180  GRAPH ( 141 ) *IY 

GOTO  (200.230)  NSPEC 
C  Sinsle  Maxwellian  fit. 

200  S1*0. 

S2*0. 

204  ’YPE  203 

203  FORMAT  ('fWhat  fraction  do  you  want  to  fit  (starting 

rfro»  riaht) y  ',*> 

ACCEPT  206.FRACT 
206  FORMAT  (F6.4) 

IF  (FRACT  .OT.  0.  .AND.  FRACT  .LE.  1.)  GOTO  208 
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207 


TYPE  207 

FORMAT  <'+Huh?  Between  0.  and  1.  rlease*'/) 
GOTO  204 

208  START*INT  < 103*( 1-FRACT) ) 

IF  (START  .LT.  NZERO/2)  START ”NZ£R0/2 
DO  210 *  I-ITO+START, IT0+102 
V*103./( I-TO)/P£RIOD*TSCAL£ 

S1*S1  +  <  DECONY < I ) -ZERO ) 4MAXUEL ( V ) 

210  S2*S2+MAXWEL ( V ) **2 

AMP-S1/S2 

DO  220.  I*IZERC+1 .311 
IF  (I  .LE.  0)  GOTO  220 
V«512 . / < I-IZERO > /PER IODATSCALE 
Y*AMP*MAXUEL ( V > +ZERO 

GRAPH  ( I +4  )  * INT  (.< Y-MIN1 )  *20 1  .  /  <  MAX  1 -MINI ) ) +1 7 
IF  ( GRAPH ( 1+4 )  .LT.  0)  GRAPH ( 1+4 ) =0 
IF  ( GRAPH ( I +4 )  .GT.  235)  GRAPH < 1+4 ) =233 


220 

CONTINUE 

GOTO  300 

C 

Double  Maxui 

230 

S1*0. 

S2-0. 

S3*0. 

S4*0 . 

S3-0. 

DO  240. I”IT0+NZER0/2 » IT0+102 
Vl»103 . / ( I-TO  > /PERIQD*TSCAL1 
V2-103 . / < I-TO ) /PERI0D*TSCAL2 
SI *S1 +MAXUEL ( VI ) **2 
S2*S2+MAXUEL ( VI ) *MAXUEL ( V2 ) 

S3-S3+M AXUEL ( V2 ) **2 

S4*S4+  <  DECONV  < I ) -ZERO ) XMAXUEL ( VI ) 

240  S3»S3+ ( DECONV  < I ) -ZERO ) XMAXUEL ( V2 ) 

DEN-S1XS3-S2X32 
A AMP» ( S4XS3-S2XS3 ) /DEN 
DAMP- ( S1XS3-S4XS2 ) /DEN 
DO  270.  I-IZERO+1.311 
IF  (I  .LE.  0)  GOTO  270 
Vl*312./( I-IZERO) /PER I 0DXTSCAL1 
V2«512. /(I-IZERO) /PERI0DXTSCAL2 
Y-AAMPXMAXUEL  <  VI ) +BAMPXMAXUEL ( V2 ) +ZERO 
GRAPH ( 1+4 )  =  INT  <(Y-MIN1)*201./<NAX1-MIN1>>+17 
IF  < GRAPH ( 1+4)  .LT.  0)  GRAPH(I+4)*0 
IF  < GRAPH ( 1+4 >  .GT.  233)  GRAPH( I+4>»233 
270  CONTINUE 

C  Determine  Standard  Deviation. 

300  VAR*0 . 

DO  330.  I*IT0+NZER0/2» IT0+102 
GOTO  (310.320)  NSPEC 
310  V-103./(I-TO)/P£Rr00*TSCALE 

TERM«AMPXMAXUEL ( V ) 

GOTO  330 

320  V1-103./<I-T0>/PERI0D*TSCAL1 

V2*103./< I-TO ) /PERIQ0XTSCAL2 
TERM-AAMPXMAXUEL ( VI ) +8AMPXMAXUEL ( V2  > 

330  VAR*VAR+ ( DECONV ( I > -ZERO-TERM  >  XX2 

C  PLOTTING  ROUTINE 

CALL  PLQT33< 13.72. ) 

CALL  PL0T33( 10. » > 

CALL  PL0T33 (1.1.) 

CALL  PLT33(2. 147.314) 

CALL  PLT33 (4 » 1 » I Y ) 

CALL  PLT33< 7 .0 .0 ) 

CALL  PLT33 <3.-312. GRAPH ) 

CALL  PL0T53(1.0.)  (Marker  eut  on  eraan  0  since 

CALL  PLT3S(4»IZ.1>  (eraeh  1  earkers  are  flaks' 
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333 


340 


330 


370 


400 


CALL  PLQTC5<  9  <  44 .0 ) 

CALL  OA  r£ <  TDATE ) 

TYPE  333*  (TDATEC I ) *  1*1 *9) 

FORMAT <  '4'  *9A1 ) 

CALL  PL0T33  <9*44*1) 

TYPE  340 

FORMAT  (  '  fRun  *'.*) 

CALL  PUTSTR  (7  *  RUNNO  .  '  4 '  ) 

CALL  PLT33 (9*44*3) 

GOTO  (330.370)  NSPEC 

CALL  SCOPY ( 'Fit  amplitude:  **'. TDATE) 

CALL  REPLCTDATE.AHP*4./EXP(2. > ) 

CALL  PLTTXT ( TDATE ) 

GOTO  400 

CALL  SCOPY < 'Coefficient  for  nass  **:  ♦*'. TDATE) 
CALL  PL TRER <  TDATE  * WEIGH 1 ) 

CALL  REPL(TBATE>AAMP*4./EXP<2. > ) 

CALL  PLTTXT (TDATE) 

CALL  PLT35 ( 9  *  44  *  4 ) 

CALL  SCOPY ( 'Coefficient  for  *a«i  ♦#:  ♦#'. TDATE) 
CALL  PLTRER ( TDATE .WEI GH2 ) 

CALL  REPL ( TDATE . BAMP*4 . /EXP ( 2 . ) ) 

CALL  PLTTXT (TDATE) 

CALL  PLT33 (9*44*3) 

UAR-SGRT (OAR/ ( 102-NZERO/2 ) ) / ( MAXI -MINI ) *100 . 
UAR-INT ( 10 . *VAR  )  * .  1 

CALL  SCOPY ( 'Standard  deviation:  MZ. ' .TDATE) 
CALL  PL TRER (TDATE. WAR) 

CALL  PLTTXT (TDATE) 

CALL  PL0T33( 13  *  72  *  ) 

CALL  PLT35 (2*1.) 

RETURN 

END 


SUBROUTINE  SAVE 

C  Saves  a  srectrue  for  DBLFIT. 

REAL  DECQNU ( 309 ) *SAWED<309> .SAMP. MAXI .MINI 
COMMON  /DECD/  OECONU . MAXI .MINI .  /SAUER/  SAUED-SAMP 
DO  10.  1-1.309 
10  SAWED ( I ) -OECONW  < I > 

SAMP— MAXI “MINI 

RETURN 

END 


SUBROUTINE  DBLFIT 

BYTE  TDATE ( 31 ) • RUNNO ( 3 ) .ERR 

COMMON  /RUNUM/  RUNNO 

REAL  MAXWEL . DECONV ( 309 ) , SAWED ( 309 ). MAXI . MINI 
MAXWEL ( V  > -EXP ( - ( V*W ) )*V**4 
INTEGER  GRAPH <3 12) 

COMMON  /Dr.CD/  DECONW , MAXI .MINI . ICHNL1 »GTO»  /SAUER/  SAUED.SAMP 
TYPE  10 

10  FORMAT  ('4This  routine  fits  the  sue  of  a  Maxwellian  and  a  wavefore 

4  stared  usins'/'  SAUE.  Is  that  wavefore  read*?  '.») 

CALL  GETSTR ( 3 . TDATE . 8 . ERR  > 

CALL  SCOMP('NO'.TDATE.IWAL) 

IF  (IVAL  .EO.  0)  RETURN 
SZERQ-0 . 

ZERO-O . 

TYPE  50 

SO  FORMAT  ( '4Flisht  rath  (ce.)!  '.*) 

ACCEPT  60. PATH 
60  FORMAT  (F9.4) 

TYPE  70 

70  FORMAT  ('+Period  (msec.):  '.*) 

ACCEPT  60.  PERIOD 
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TYPE  80 

80  FORMAT  (  '  YE'Juctor  (n  sac .  >  I  ',♦) 

ACCEPT  60.  DELAY 

T0-GT0-  <  DELAY+3  .  )*.  103/PERIOD 

PERI0D-PERI00/1000. 

ITO-INT (To+l . ) 

IZERQ* INT ( 3  .  *  (  T0-ICHNL1 > ) 

IZ-IZERO 

IF  <  IZERC  .LT.  0)  IZERO-O 
TYPE  90 

90  FORMAT  (  '  >Tein**ratijr@  (cantiaraMa )  :  '»*) 

ACCEPT  60.  TEMP 
TEMP-TEMP+273. 13 
TYPE  110 

110  FORMAT  (  '  +Malecular  waiShtJ  '»*> 

ACCEPT  60. WEIGHT 

TSCALE*S0RT(UEIGHT/TEMP*6.01659E-9)*PATH 
TYPE  160 

160  FORMAT  ('+-How  i»an*  channels  ara  ’zara*?  ',») 

ACCEPT  163>NZER0 
16S  FORMAT  ( 13  > 

DO  170.  I»ITO-NZERO. ITO+NZERO-1 
SZERO*SZERO+SAVED( I ) 

170  ZERO*ZERO+DECONW< I) 

SZERO-SZERO/NZERO/2 

ZER0*Z£R0/NZER0/2 

IY*INT( (ZERO-MIND *201 ./(MAXI-MINI >  )  +  l7 
DO  180.  1*0 » IZERO+3 
180  GRAPH ( X+4 ) *1 Y 

C  FITTING 

51- O. 

52- 0. 

53- 0. 

S**0. 

S3-0. 

DO  260,  I-ITO+NZERO/2 . IT0+102 
V-103./I  I-TO)/PERIOD*TSCALE 

51- S1+MAXUEL ( V  >  **2 

52- S2+MAXUEL ( V ) * ( SAVED ( I ) -SZERO ) 

53- S3+ ( SAVED ( I ) -SZERO )  **2 

54- S4+ ( DECONV ( I ) -ZERO ) 4MAXWEL <  V ) 

260  S3-S3+ <  DECONV  < I ) -ZERO ) *  <  SAVED ( I > -SZERO ) 

DEN-S1*S3-S2*S2 
AAHP*  <  S4*S3-S2*S3 ) /DEN 
8AMP*  <  S1*S3-S4*S2 ) /DEN 
DO  300, I-IZERO+1 .311 
IF  <1  .LE.  0)  GOTO  300 
J-I/3+ICHNL1 
K»I-I/3*3+l 

GOTO  (270,270,270,280.280)  K 

270  TERM-FLOAT (K+2 ) /3 .*( SAVED ( J ) -SAVED  <  J-l ) )+SAVED(J-l ) -SZERO 

GOTO  290 

280  TERM-FLOAT ( K-3  > /3 . * ( SAVED <  J+ 1 ) -SAVED  <  J) > +SAVED ( J ) -SZERO 

290  V-512 . / ( I-IZERO ) /PERIQD*T3CALE 

Y-AAMP*MAXWEL ( V ) +8AMP*TERM+ZER0 
GRAPH ( I+4)»INT( ( Y-MIN1 )*201 ./(MAXI-MINI ) )+17 
IF  ( GRAPH( 1+4 )  .LT.  0)  GRAPH<I+4>-0 
IF  C0RAPH(I+4)  .GT.  233)  GRAPH < 1+4 ) *233 
300  CONTINUE 

VAR-O. 

DO  330, I-IT0+NZER0/2, IT0+102 
V*l03./( I-TO  > /PERIOD* TSC ALE 
TERN*AAMP*MAXUEL ( V ) +BAMP* <  SAVED ( I > -SZERO > 

330  VAR-VAR+ ( DECONV ( I ) -ZERO- TERM ) **2 

C  PLOTTING  ROUTINE 

CALL  PLQT33< 13,72, ) 
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CALL  PL0T33 ( 10. . ) 

CALL  PLT33 (1.1. > 

CALL  PLT33<2,167.344) 

CALL  PLT35 ( 4  « 1 » IY  > 

CALL  PLT33 (7.0.0) 

CALL  PLT35(3. -312, GRAPH) 

CALL  PLT33 (1.0.) 

CALL  PLT33 (6.12.1) 

CALL  PL0T33  <9.44.0) 

CALL  DATE ( TDATE ) 

TYPE  333. TDATE 
333  FORMAT  ('+'.9A1) 

CALL  PL0T33<9.44, 1 ) 

TYPE  340 

340  FORMAT  ( ' fRun  *'.*) 

CALL  PUTSTR  (  7  .  RUNNO  »'+•') 

CALL  PLT33 (9.44,3) 

CALL  SC0PY(  'Coefficient  for  Maxwellian.'  #♦ '  » TDATE) 

CALL  REPL (TDATE. AAMP*4./EXP<2.  )  ) 

CAL!.  PLTTXT( TDATE) 

CALL  PLT33 ( 9 , 44,4) 

CALL  SCOPY( 'Coefficient  for  other  *eak :  *4', TDATE) 

CALL  REPL  <  TDATE , 8AMP4SAMP ) 

CALL  PLTTXT ( TDATE ) 

CALL  PLT33(9,44,5) 

CALL  SCOPY( 'Standard  deviation!  #*X . ' » TDATE ) 

RNUM-INT  < 10*SORT ( VAR/ ( 1 02-NZER0/2 ) ) / ( MAX 1-MINI >  *100 . >  * . I 
CALL  PLTRER(TDATE.RNUM) 

CALL  PLTTXT (TDATE) 

CALL  PL0T33 (13,72.) 

CALL  PLT33 (2.1.) 

RETURN 

END 

SUBROUTINE  TRAP IT 

C  This  routine  functions  essentially  as  MAXFIT,  except  that  the  first 

C  Maxwellian  can  be  a  translating  Maxwellian. 

BYTE  TDATE ( 81 ) « RUNNO ( 3 ) . ERR 
COMMON  /RUNUM/  RUNNO 

REAL  MAXUEL » TRANS, DECONV( 309) , MAXI  .MINI 
MAXUEL ( V ) »EXP ( - ( V*V  ) )*V**4 
TRANS ( V » U ) »EXP  < - < ( V-U )  **2 ) >*V**4 
INTEGER  GRAPH (312) 

COMMON  /DECD/  DECONV , MAXI .MINI , ICHNL1 ,GTO 
3  TYPE  10 

10  FORMAT  ('+How  eanw  seecies?  ',*) 

ACCEPT  20.NSPEC 
20  FORMAT  (13) 

IF  (NSPEC  .EQ.  1  .OR.  NSPEC  .EQ.  2)  GOTO  40 
TYPE  30 

30  FORMAT  ('+1  or  2  only!'/) 

GOTO  3 
40  ZERO-O . 

TYPE  SO 

30  FORMAT  ('+Flisht  eath  (ca.K  '»*> 

ACCEPT  60. PATH 
60  FORMAT  <F9,4) 

TYPE  70 

70  FORMAT  ('+Period  (msec.)!  ',*> 

ACCEPT  60, PERIOD 
TYPE  80 

80  FORMAT  <'+Eduetor  delay  (usee.)!  '.*> 

ACCEPT  60. DELAY 

TO*GTO-<  DELAY+3. ) * . 103/PERIOD 

PER I OD*P£R 10  0/ 1 000 . 

ITO*INT ( TOf  1 • ) 
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IZ£R0»INT<3.*<  T0-ICHNL1 )  ) 

IZ-IZERO 

IF  ( IZERO  .LT.  0)  IZERO-O 
TYPE  33 

S3  FORMAT  < ' +Mo  leojlar  uai3htl  '»♦) 

ACCEPT  60, WEIGHT 
GOTO  <35 >120)  NSPEC 
33  TYPE  90 

90  FORMAT  (  '  rTa—aratura  CcantiSrada )  t  '»*) 

ACCEPT  60,  TEMP 
TEMP*TEMP+273.13 

TSCALE-SORT  <  HEIGHT/ TEMP*6 . 01 639E-9 ) *PATH 
TYPE  113 

113  FORMAT  <'+Mach  number:  ' ,*) 

ACCEPT  60 ,  SR 

U*SR*1 . 12837?  ! Sat*  SR  (Mach  no.)  to  b»  ratio  of  U  to  «auara 

GOTO  130  Toot  of  <  8k  T/C* i 3* ) 

120  TYPE  130 

130  FORMAT  < ' +Tran*latina  Ma:(wallian —  ' 

+/'  Ta»-aratura  ( cantis rad* ) I  '»*> 

ACCEPT  60, TEMP 
TEMP-TEMP+273.13 

TSCAL1-SGRT  <  U£1GHT/TEMP*6 . 01639E-9  )  *PATH 
TYPE  133 

133  FORMAT  < ' +  Mach  mjabar:  ',»> 

ACCEPT  60, SR 
U»SR*1. 128379 
TYPE  140 

140  FORMAT  t'+Raaular  Maxwall lan—  V '  Taa-aratijra :  ',*) 

ACCEPT  60, TEMP 
TEMP-TEMP+273.13 

TSCAL2-SQRT  <  UEXGHT/TEMP«6 . 016S9E-9 ) 4PATH 
150  TYPE  160 

160  FORMAT  <'+Hou  ran*  channel*  art  'zero’?  ',*> 

ACCEPT  20.NZER0 

OQ  170,  I-1T0-NZERQ. ITO+NZERQ-1 
170  ZERO-ZERO+DECONV  < I ) 

ZERO-ZERO/NZERO/2 

IY*INT<  < ZERO-MINI ) *201 ./ (MAXI-MINI ) )+17 
DO  130,  1*0 , IZERQ+3 
ISO  GRAPH<I+1)*IY 

GOTO  (200,230)  NSPEC 
C  SINGLE  MAXHELLI AN  FIT! 

200  S1«0. 

S2-0. 

204  TYPE  203 

203  FORMAT  f'+Uhat  fraction  do  aou  want  to  fit  (atartins 

+fi*oa  riaht)?  '»*) 

ACCEPT  206 ,FRACT 

206  FORMAT  (F6.4) 

IF  (FRACT  .GT.  0.  .ANX).  FRACT  .LE.  1.)  GOTO  208 
TYPE  207 

207  FORMAT  ( '+Huh?  Between  0.  and  1.  alaaaa!'/) 

GOTO  204 

208  START-INT  < 103* < 1 -FRACT ) ) 

TMAX— 100000. 

IF  (START  .LT.  NZERO/2)  START-NZERO/2 
00  210.  I-1T0+START.IT0+102 
V-103 . / ( I-TO ) /PERIODFTSCALE 
TMAX-AMAX1 ( TMAX , TRANS  <  V , U ) ) 

SI -S 1 + (  DECONV < I ) -ZERO ) * TRANS ( V ,  U  > 

210  S2-S2+TRANS  ( V , U ) **2 

AMP-SI /S2 

00  220.  I-IZERO+1.311 

IF  (I  .LE.  0)  GOTO  220  . 

V-312./(I-IZERO)/PERI0D*TSCALE  f 
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y»AMP*TRANS> V.U) +  ZERO 

GRAPH  <  I  >4  )  •=  INT  (  < Y-MIN1  >  *20 1 .  /  C MAXI  -MIN  t )  )  +  t7 
IF  ( GRAPH  < 1  +  4 )  .LT.  0)  GRAPH  ( I  +4 ) *0 
IF  (GRAPH!  1  +  4  )  .GT.  233)  GRAPH  (  1+4 )  *233 
220  CONTINUE 
GOTO  300 

C  DOUBLE  MAXWELLIAN  FIT: 

230  51*0 . 

S2*0. 

S3*0. 

S4«0. 

S3*0. 

TMAX*-100000. 

DO  260 • I* IT0+NZER0/2, ITO+ 102 
Vl*l03./( I-T0)/PERIOD*T3CAL1 
V2*103 . / ( I-TO  > /PERI0D*T3CAL2 
TMAX*AMAX1 < TMAX , TRANS ( VI ,U) > 

S1»S1+TRANS( VI » U ) **2 
S2*S2+TRANS ( VI »U ) *MAXU£L ( V2> 

S3*S3+MAXUEL ( V2 > **2 

S4*S4+  (  DECONV  < I ) -ZERO > #TR ANS ( V 1 , U ) 

260  S3*S3+ ( DECONV  < I ) -ZERO ) *MAXWEL ( V2  > 

DEN»S1 *S3-S2*32 
AAMP* ( S4*S3-S2*S3 > /DEN 
BAMP- <  S 1*S3-S4*S2 ) /DEN 
DO  270,  I*IZER0+1 ,311 
IF  <1  .LE.  0)  GOTO  270 
Vl«312 . / ( I-IZERO) /PER10D*TSCAL1 
V2*512  •  /  (  I-1ZER0 )  /P£RI0D'*TSCAL2 
Y«AAMP*TRANS  <  VI ,U ) +BAMP*MAXUEL  (  V2 ) +ZERO 
GRAPH ( I+4)*INT(( Y-M1N1 ) *201 • / ( MAXI— MINI ) >+17 
IF  < GRAPH ( I +4 )  .LT .  0)  GRAPH< 1+4  >  =0 
IF  ( GRAPH < 1+4 )  .GT.  233 >  GRAPH! 1+4 > *233 
270  CONTINUE 

300  VAR-0 . 

DO  330,  I-IT0+NZER0/2, IT0+102 
GOTO  (310.320)  NSPEC 
310  V-103./(I-T0)/PERI0D*TSCALE 

TERM»AMP*TRANS  <  V , U ) 

GOTO  330 

320  Vl-iO3./(I-TO)/PERI0D*TSCALl 

V2-103./(I-T0)/PERI0D*TSCAL2 
TERM«AAMP*TRANS ( VI , U  >  +BAMP4MAXUEL  <  V2  > 

330  VAR-VAR+  <  DECONV ( I > - ZERO-TERM ) **2 

C  PLOTTING  ROUTINE 

CALL  PL0T33< 13,72.) 

CALL  PLQT33< 10 . , > 

CALL  PLT33 (1,1.) 

CALL  PLT33 (2,167,314) 

CALL  PLT33( 4, 1 , IY ) 

CALL  PLT33(7.0,0) 

CALL  PLT33(3, -312, GRAPH) 

CALL  PL0T33( 1,0.) 

CALL  PLT33(6,IZ,1) 

CALL  PL0T33 (9,44,0) 

CALL  OATE(TDATE) 

TYPE  333.TDATE 
333  FORMAT  (,+',?Ai) 

CALL  PL0T33( 9,44,1) 

TYPE  340 

340  FORMAT  <'+Rur.  »',*> 

CALL  PUTS  TR ( 7 , RUNNO , '  + ' ) 

CALL  PLT33(9,44,3) 

GOTO  (330.370)  NSPEC 

330  CALL  SCOPY (  'Fit  1  i 5  M'.TDATE) 

CALL  REPL(TDATE,AMP*TMAX) 
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CALL  PLTTXT ( TDATE ) 

CALL  PLT33 (9.44.6) 

CALL  SCOPY(  'T«»Mntur«:  M '  .  TDATE ) 

CALL  PLTREI ( TDATE. INT < TEMP ) > 

CALL  PLTTXT (TDATE) 

C  Print  desree  sian  an  terminal  end  Plotter. 

CALL  PL0T33  (13.70.  ) 

TYPE  366 

366  FORMAT  ( ' YfK  ,  '.*> 

CALL  PL0T33<  13.71 » ) 

CALL  PLTSTA(ISTAT) 

GOTO  (363.367)  I3TAT+1 

367  CALL  3C0PY( ' UC1 . 3 . 99 . 1 » 1 .0 , 1 . - 1 . 1 . -1  > - 1 . 0 . - 1 . 

+1.-1.-99JL3K  > ' , TDATE ) 

CALL  PLTSND( TDATE) 

C  Continue  athsr  printine. 

368  CALL  3C0PY ( ' Mach  *:  M'.  TDATE) 

CALL  PLTRERC TDATE. 3R) 

CALL  PLTTXT (TDATE) 

GOTO  400 

370  CALL  SCOPY( 'Coefficient  for  Maxwellian:  #♦' .TDATE) 

CALL  REPL ( TDATE . BAMPS4 ./EXP ( 2. ) ) 

CALL  PLTTXT (TDATE) 

CALL  PLT33( 9 .44 .4 ) 

CALL  SCOPY( 'Coefficient  for  othar  meek:  M'. TDATE) 

CALL  REPL ( TDATE » AAMP4TMAX ) 

CALL  PLTTXT (TDATE) 

400  CALL  PLT33( 9 .44.5) 

VAR«SQRT (VAR/( 102-NZER0/2) ) / (MAXI -MINI > *100 . 

VAR«INT(VAR#10. )*. 1 

CALL  SCOPY ( ' Standard  deviation:  MX. '.TDATE) 

CALL  PLTRER( TDATE. VAR) 

CALL  PLTTXT (TDATE) 

GOTO  (413.423)  NSPEC 
413  CALL  PL0T33( 9.33. 10 ) 

C  Allow  chanae  of  temperature  1/or  MacM  number  without  aivins  all  else. 

TYPE  420 

420  FORMAT  (  '+Another  temperature  l/or  Mach  number’’  '»*> 

CALL  GET3TR  (3 . TDATE . 3 . ERR > 

CALL  SC0MP< 'YES', TDATE, IVAL) 

IF  ( IVAL  .EO.  0)  GOTO  430 
423  CALL  PL0T33 (13.72.) 

CALL  PL*T33(2»1.) 

RETURN 

430  TYPE  90 

ACCEPT  60. TEMP 
TEMPmTEMP+273 . 13 

TSCALE*SQRT (UEIGHT/TEMPS6 . 01659E-9 ) 4PATH 

TYPE  113 

ACCEPT  60. SR 

U-SRal. 128379 

S1«0. 

S2-0. 

GOTO  208 
END 

SUBROUTINE  SUBMAX (ITYPE) 

C  Subtracts  a  translating  or  static  Maxwellian  from  data.  Most  details 

C  are  as  thaw  are  in  MAXFIT  and  TRAFIT. 

BYTE  TDATE (81). RUNNO ( 3 ) • ERR 
COMMON  /RUNUM/  RUNNO 

REAL  MAXWEL . DECONV ( 309 >  > MAX 1 .MINI .DIFFER (309) .MAX2.MIN2 
MAXWEL ( V . U ) »£XP ( - ( ( V-U )  »*2  > >*V**4 

INTEGER  GRAPH (312)  | 

COMMON  /DECO/  DECONV, MAXI. MINI .ICHNL1.GT0  5 

40  ZERO-O.  1 
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TYPE  50 

50  FORMAT  t  '  fFlisht  -atn  '»*> 

ACCEPT  60 • PATH 
60  FORMAT  (F9.4) 

TYPE  70 

70  FORMAT  <'+Period  (»s*e.>:  '•*) 

ACCEPT  60 .PERIOD 
TYPE  80 

80  FORMAT  C'+Eductor  d#l*%  Ojs«c.)5  '»*) 

ACCEPT  60. DELAY 
T0*GT0-<DELAY+3 . >*. 103/PERI0D 
PERIOD-PERIOD/lOOO . 

ITO*INT ( TO+l . ) 

TYPE  90 

90  FORMAT  ( '+T«»*eratijr»  i  c#ntiarad*  >  1  '»*) 

ACCEPT  60. TEMP 
TEMP-TEMP+273. 15 
GOTO  (95.100)  ITYPE 
95  U*0 . 

GOTO  105 

100  TYPE  101 

101  FORMAT  <'+M*ch  nuaber!  '«•; 

ACCEPT  60. SR 

U-SR*1. 12337? 

105  TYPE  110 

110  FORMAT  <  '  +Molactjlar  ueiaht:  '»*> 

ACCEPT  60.UEIGHT 

TSCALE-SQRT ( WEIGHT /TEMP*6 . 0 1639E-9 >  *PATH 
AHPMAX—10000. 

DO  120.1*1.103 
0*103 . /I/PERIODATSCALE 
120  AMPMAX*AMAX 1 <  AMPMAX  » M AXUEL ( V . U ) ) 

185  TYPE  190 

190  FORMAT  ( '+Maxuallian  aa^litudaS  '»*) 

ACCEPT  60.  AMP 
AMP*AMP/ AMPMAX 

C  SUBTRACTION  BEGINS  HERE. 

MAX2*- 10000. 

MIN2*10000 . 

DO  210.  I*IT0+3 » IT0+105 
0*103 . / < I -TO ) /PERI0D*T3CALE 
DIFFER! I >*DECONU<  I )  -AMP*MAXWEL< V.U> 

IF  <1+206  .LE.  309)  DIFFER < 1+206 ) -DIFFER < I ) 

IF  <1+103  .LE.  309)  DIFFER < 1+103 ) -DIFFER < I ) 

IF  <1-103  .GE.  1)  DIFFER < 1-103 ) *BIFFER < I ) 

MAX2-AMAX 1 <  MAX2 . 0 1 PFER  < I > ) 

210  MIN2*AMIN1 (MIN2. DIFFER! I ) ) 

DO  220.  IX-0.505.S 

1Y-INT  t ( DIFFER  < ICHNLl  +  IX/5)-MIN2) *201 ./ <  MAX2~MIN2 ) )  +  l7 
DO  220.  1*1.5 
220  GRAPH! IX+I )*IY 

IY*rNT< < DIFFER <ICHNL1+102)-MIN2) *201. /( MAX2-MIN2 )) +17 
GRAPH<511 )*IY 
GRAPH! 512 )*IY 
C  PLOTTING  ROUTINE 

CALL  PL0T55< 13.72. > 

CALL  PL0T55! 10. . ) 

CALL  PL0T35( 1.0.) 

CALL  PL0T55< 2.512. ) 

CALL  PLT55  <2.11. 500 ) 

CALL  PLT55<7.0.0) 

CALL  PLT53<3»-512.GRAPH) 

CALL  PLT55<9.A4.0) 

CALL  DATE<TDATE> 

TDATE<10)*0 
CALL  PLTTXT<TDATE) 
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call  plT33(9,44. i ) 

CALL  SCQPY  <  '  Run  ♦ '  , T  DA  TE ) 

CALL  CONCAT ( TDA  TE  , RUNNO  ,  TDA  TE ) 

CALL  PLTTXT<TDATE) 

CALL  PLT33  (9,44.2) 

CALL  SCOPY (  'Full  scale  aeel  i  tude  :  *♦ ' » TDATE ) 

CALL  REPL ( TDATE »  MAX2-MIN2 ) 

CALL  PLTTXTf TDATE) 

CALL  PLQT33<  9  » 33 » 10 ) 

C  Allow  another  a»r litude  choice* 

TYPE  400 

400  FORMAT  ('tAnother  try?  ',*> 

CALL  GETSTR (3 »  TDATE , 3 .ERR ) 

CALL  SCOMP( 'YES' .TDATE, IUAL) 

IF  ( XVAL  .£0.  0)  GOTO  133 

C  Replace  appropriate  values  with  their  new  values. 

MAX1-MAX2 
MIN1-MIN2 
DO  410.  1-1.309 
410  DECQNV  <  I ) -DIFFER ( I ) 

RETURN 

END 

SUBROUTINE  PEAKS 

C  Picks  out  peaks  by  looking  for  chanse  in  direction  (averaged  over  a 

C  few  channels).  Also  calculates  average  velocity. 

REAL  DEC0NVC309) , MAXI. MINI 
COMMON  /DECD/  DECONV. MAXI .MINI , ICHNL1 .GTO 
40  ZERO-0 . 

TYPE  30 

30  FORMAT  < '+Flisht  path  <c*.)l  ' .*> 

ACCEPT  60, PATH 
60  FORMAT  <  F9 . 4 ) 

TYPE  70 

70  FORMAT  (' 4-Period  (msec.)l  ’,») 

ACCEPT  60. PERIOD 
TYPE  80 

80  FORMAT  ('4-Eductor  do lav  (usee.)!  '»*) 

ACCEPT  60. DELAY 

TO-GTO-C BELAY 4-3.  ) 4.103/PERIOD 

PERIOD-PERIOD/ 1000 . 

ITO—INT (T04-1 .  ) 

TYPE  160 

160  FORMAT  < '+How  many  channels  are  ‘aero’?  '»*> 

ACCEPT  163.NZER0 
163  FORMAT  (13) 

DO  170,  I-ITO-NZERO » IT0+NZER0-1 
170  ZERO-ZERO+DECONV ( I ) 

ZER0-ZER0/NZER0/2 
CALL  PL0T33< 13,72.) 

C  Erase  text  for  clarity.  Keep  identifying  labels. 


TYPE  173 

173  FORMAT  < ' +  '/ 

♦  '  '/ 

+'  '/ 

+  '  '/ 

4-'  V 

♦  '  '/ 

♦  '  '/ 

♦  '  V 

♦'  '/ 

+'  '/ 

♦  '  '/ 

♦  '  '/ 

'/ 

♦  '  '/ 


+  ' 

CALL  PLQT33  (  9 » 0 « 1 A  ) 

IFLAG-0 

C  IFLAG  identikits  'direction  of  sloe*.  6  »oints  on  either  side  of 

C  current  location  are  sij»»ed . 

SUMLQ-OECONU ( ICHNL 1 1 +OECGNY  < ICHNL 14-1) +DECONU < ICHNL 1+2) 

++DECONU ( ICHNL1 +.3  ) +DEC3NU ( ICHNL1 +4 ) +DECONU ( ICHNL 1+3) 

SUMHI -0EC3NU ( ICHNL 1 + 7 ) +DEC0NU ( ICHNL 1 +3 ) +DEC3NL ( ICHNL1+9 ) 

++DECONU ( ICHNL1 +1 0 ) +OECONU  < ICHNL1  +  11 ) +DECGNU ( ICHNL1  +  12) 

DIFF-SUMHI-5UMLG 

IF  (DIFF  .GT.  0)  IFLAG- 1 

DO  190 »  I * ICHNL 1 +6 » ICHNL 1+70 

3UML0-SUML3-DEC0NU< I-4)+DEC0NU< I) 

SUMHI»SUMHI-OECQN<J( 1+1 >+DEC0NV< 1+7) 

DIFF  -SUMHI -SUMLG 

IF  < IFLAG  .EG.  0  .AND.  DIFF  .LT.  0 
+  .OR.  IFLAG  .EG.  1  .AND.  DIFF  .uE.  0)  30+0  190 
IF  (IFLAG  .EQ.  0  .AND.  DIFF  .GE.  0)  GOTO  133 

TYPE  180. I  +  1-ICHNL1 . ( I  +  1 . 3- TO  )  /l . 03E-44F £R I0D  .  DECONV < I  +  1 ) -ZERO 
130  FORMAT  (  '  +Peak  in  channel '« 14 »  '  ('.F3.0»'  m icrosecori.d» >  with 

+  amplitude!  '.F9.0/) 

XFLAG-0 

CALL  PLQT33 ( 1 >  0  * ) 

CALL  PL0T33 (2. 128  >  ) 

C  Mark,  peaks. 

CALL  PL0T33C 6 » 3*  ( I+1-ICHNL1J -3.1) 

GOTO  190 
183  I FLAG- 1 

190  CONTINUE 

C  Calculate  and  print  average  velocity. 

51- O. 

52- 0. 

DO  210.  I-IT0+NZER0/2. ITO+102 
V-PATH/(I-T0)*103. /PERIOD 
Sl-Sl+U* ( DECQNU  < I ) -ZERO ) 

210  S2-S2+0ECQNV ( I ) -ZERO 

TYPE  220.S1/S2 

220  FORMAT  ('+Avera*e  velocity!  '»F7.0.'  centimeters  per  second.') 

CALL  PLQT33 (13.72.  ) 

RETURN 

END 

SUBROUTINE  SMOOTH 
BYTE  TDATE ( 81 ) . RUNNO (3 ) .ERR 
COMMON  /RUNUM/  RUNNO 
INTEGER  GRAPH (312) 

REAL  IN(309> »0UT(30?) .MAXI .MINI 

COMMON  /DECD/  IN. MAXI .MINI , ICHNL1  /BNS/  9M103) 

TYPE  3 

3  FORMAT  ('  This  routine  smooths  data  by  disital  low-pass 

+  filtering. ' ) 

C  Lambda  is  set  to  30.  dB.  Delta  is  set  to  about  2.3  components. 

C  This  sets  MP  to  103.  aivina  aood  accuracy  and  reasonable  run-time. 

C  See  Kaiser  and  Reed.  R.S.I.  4311447  (1977)  for  details  and  for 

C  modifications. 

5  TYPE  10 

10  FORMAT  ('  How  many  Fourier  components  to  be  kept?  '.*> 

ACCEPT  20.  BETA 

20  FORMAT  (F8.2) 

IF  (BETA  .LT.  32.)  GOTO  22 
TYPE  21 

21  FORMAT  ('-^Maximum  number  is  32!') 

00T0  3 

22  BETA-2. *BETA/103. 

IFLAG-0 

CALL  NER(BETA) 
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00  40 , N* 1 , 103 
OUT ( N  >  «0 • 

00  30.  K*-l03 > 103 

0UT(N>-aUT(N)+6K<  rABS<K)*l )*IN<N+103-K) 

QUT(N+206)*aur<N) 

40  0UT<N+103)»0UT(N) 

C  PLOTTING  ROUTINE 

DO  30.  IX»0 » 305 . 5 

IV*INT<  <OUT< ICHNLl+rx/3) -MINI) *201 ./(MAXI -MINI )  >  +  17 
SO  SO.  1*1.3 
50  GRAPH(IX+I>*IY 

IY»INT < ( OUT ( ICHNL1+102)-MIN1 >*201 ./(MAXI-MINI >  >  +  17 

GRAPH ( 3 1 1  >  *  I Y 

GRAPH < 312 >»IY 

CALL  PL0T33 (13.72.) 

Call  plotssco,,) 

CALL  PLOTSSa.l.  > 

CALL  PLT33 (2.31. 500 ) 

CALL  PLT33 17.0.0) 

CALL  PLT35 (3.-312. GRAPH ) 

CALL  PLT35 (9.44.0) 

CALL  SATE (TOATE) 

TDATE(10)»0 
CALL  PLTTXT(TDATE) 

CALL  PLT53 (9.44.1) 

CALL  SCOPV ( ' Run  t' .TOATE) 

CALL  CONCAT ( TOATE. RUNNO.TD ATE) 

CALL  PLTTXT(TDATE) 

CALL  PLT33  <9,44.2) 

CALL  SCOPY t  'Full  seal*  am«»liUid*:  M' .TOATE) 

CALL  REPL(TDATE.MAXl-HINl) 

CALL  PLTTXT< TOATE) 

CALL  PL0T33  <9,33,10) 

TYPE  90 

90  PORMAT  ('+AnotM*r  tr*?  ',») 

CALL  GETSTR 13, TOATE » 8 » ERR ) . 

call  scomp< -yes' .tdate.ivad 
if  <ival  .ea.  o>  goto  3 

CALL  PL0T33  <  V » 33 . 1 1 ) 

C  Either  r*i»lac*  DECONV  with  smoothed  version,  or  simAl*  abort  ro.jtin«. 

TYPE  100 

100  PORMAT  ('+S#v*  smoothed  version?  '.») 

CALL  GETSTR <3, TOATE, 3, ERR) 

CALL  SCOMPC 'YES' .TOATE » I VAL) 

CALL  PLT33 (2,11, 300 ) 

IF  (IVAL  .NE.  0)  RETURN 
CALL  PL0T35U.0.) 

CALL  PL0T33<3. -312, GRAPH) 

00  110,  1*1,309 
110  IN( I ) *QUT  < I ) 

RETURN 

END 

SUBROUTINE  NCR (BE) 

COMMON  /BKS/  BK ( 103) 

PI*3. 141392634 
FK-10.0344 
ET*7. 83726 
CALL  INO(ET.FIA) 

DO  10,  K-1,103 
8K-PI*FL0AT(K> 

3E«ET*SQRT ( 1 . - ( FLOAT ( X ) f 1 03  > **2 ) 

CALL  INO(GE.E) 

BK(K)*(3IN< BE*GK ) /GK >* (EXFIA) 

10  CONTINUE 

DO  20,  1*2 .104 


109 


K-104-r-t-c 

J-K-l 

20  3K(K)-BK(J> 

BK(1  )»BE 

RETURN 

END 

SUBROUTINE  INO(X.S) 

S-l. 

OS-1 . 

D-O. 

1  D-D+2. 

os-os*x*x/(d*D) 

S-S+DS 

IF  (OS  .ST.  0.2£-3*S)  GOTO  1 

RETURN 

END 

SUBROUTINE  READER 

C  Read  in  new  spectrum  via  a/d  or  old  from  disk. 

INTEGER  DAT  A ( 206 ) » I DATA (  100 ) 

BYTE  ANSWER ( 3 > » RUNNO ( 3 ) .TEMP(U)  .FILNAHC 13) >ERR,NAM£(7) 

COMMON  /RAUD/  DATA . MIN . MAX  .  /RUNUM/  RUNNO 

ICMF-0 

MIN-10000 

MAX— 10000 

C  Set  values  for  missina  channels  to  Oe  interpolated. 

DATA< 101 >=2048 
DATA<  102)=»2048 
DATA< 103) *2048 
3  TYPE  10 

10  FORMAT ('+NEU  or  OLD''  '.») 

CALL  GETSTR  <  3 » ANSUER  *  4 . ERR  ) 

CALL  SC0MP( " NEU ' .ANSWER  » I UAL ) 

IF  (IVAL  .EQ.  0)  GOTO  100 
CALL  SC0MP( 'OLD' .ANSWER. I UAL) 

IF  (IUAL  .EO.  0)  GOTO  200 
TYPE  20 

20  FORMAT  ('+What?'/) 

GOTO  3 

100  TYPE  110 

110  FORMAT  ( ' +Run  no.t  '»*> 

C  Create  temporary  data  file  • DX1 : RUn . DAT* 

CALL  GETSTR ( 3 .RUNNO . 4 .ERR > 

CALL  CONCAT ( ' 0X1 :RU ' .RUNNO . TEMP ) 

CALL  CONCAT ( TEMP  » ' .  DAT ' .FILNAM ) 

CALL  ASSIGN( 2.FILNAM .0 ) 

TYPE  120 

120  FORMAT ( '  Read*  with  data?  '.»> 

CALL  GETSTR ( 5  *  TEMP » 1 0 » ERR ) 

CALL  RTS( I DA TA.l 00 »»»».. 4. XCMF» IBEF )  (read  in  data  at  clock  rate 
CALL  LWAIT(ICMF.Q)  (provided  at  ST1  input. 

DO  123.  1-1.100 
123  DATA< I)-IDATA<I) 

CALL  PUT3TR  (  2 » RUNNO . 0 ) 

WRITE  (2.130)  (DATA(I) .1-1.103) 

130  FORMAT  (2313) 

GOTO  300 
200  TYPE  210 

210  FORMAT ( ' fFi le  name:  '»*> 

ERR-. FALSE. 

CALL  GETSTR (3. NAME. A. ERR) 

IF  (.NOT.  ERR)  GOTO  230 
TYPE  220 

220  F0RMAT('+6  characters  msx imum  for  file  name!'/) 

GOTO  200 
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23C  CALL  CONCAT; ' 3X1 t ' .NAME. TEMP > 

CALL  CQNCA T ( TEMP » ' .DAT'  .FILNAM ) 

CALL  ASSIGN(2. FILNAM. 0) 

232  TYPE  110 

CALL  GETSTR ( 3  * ANSWER  »  4  .ERR  ) 

235  READ  ( 2 » 240 . £ND=260 )  RUNNO(l) 

240  FORMAT  !A1) 

BACKSPACE  2 

CALL  GETSTR  <  2  .RUNNC >  4 . ERR ) 

READ  (2.250) ( DATA ( I ) . 1*1 > 103) 

230  FORMAT ( 2313 ) 

IF  ( ISCCMP(RUNNO. ANSWER > >  233.300.233 
260  TYPE  270 

270  FORMAT ('+No  such  rur.no.'  /) 

REWIND  2 
GOTO  232 

300  CALL  CLOSE ( 2 ) 

DO  310.  1=1.103 
data; i+io3)=data( r  > 

MIN=(1IN0;  MIN*  data;  I  )  ) 

310  MAX=MAXO< MAX. DATA ( I )  ) 

RETURN 

END 


SUBROUTINE  JOIN 

C  Produce  composite  file  far  a  set  of  temporary  files  by  reading  in 

C  each  file  and  writing  them  out  into  one  new  file.  This  is  necessary 

C  to  avoid  to  inefficient  concatenation  achieved  by  the  monitor  COPY 

C  command,  which  does  not  eliminate  blank,  spaces. 

BYTE  ANSWER; 11! .ERR. IVAL. JUAL.FILNAM; 15) .TEMP; 11 ) 

3  TYPE  10 

10  FORMAT  C'+Outmut  filet  '.*> 

CALL  NAME; ERR. FILNAM) 

IF  tERR  .EQ.  .TRUE.)  GOTO  5 
CALL  ASSIGN<2. FILNAM, 0) 

13  TYPE  20 

20  FORMAT  ('+Are  files  consecutively  numbered  starting  with  1?  ',*) 

CALL  GETSTR;3. ANSWER, 10. ERR) 

CALL  SCOMP ; ' YES ' , ANSWER , IVAL ) 

IF  (IVAL  .£0.  0)  GOTO  30 
CALL  SCOMP ( ' NO ' . ANSWER  » J UAL ) 

IF  ( JVAL  .EQ.  0)  GOTO  30 
TYPE  23 

25  FORMAT  ; '+What?'/) 

GOTO  3 

30  TYPE  AO 

40  FORMAT  ; ' -HHou  many  files?  ',*) 

ACCEPT  30,  NFILES 
30  FORMAT  ;I3) 

IF  ( IVAL  .NE.  0)  GOTO  200 
IF  (NFILES  .LT.  100)  GOTO  70 
TYPE  60 

40  FORMAT  ('^Maximum  of  99  consecutive  files  allowed.'/) 

GOTO  13 

70  DO  100.  1=1, NFILES 

IF  <1  .GE.  10)  GOTO  SO 

CALL  CQNCAT ( '0X1 tRU' ,49+1 .TEMP) 

GOTO  90 

80  CALL  CONCAT; 'DX1 tRU' ,48+1/10, FILNAM) 

CALL  CONCAT (FILNAM»40+I-1O*( 1/10) .TEMP) 

90  CALL  CONCAT(TEMP. ' .DAT' .FILNAM) 

100  CA»  TRANSF; FILNAM) 

CA..L  CLOSE; 2) 

RETURN 

200  DO  230.  1*1. NFILES 
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210 

220 


TYPE  220.1 

FORMAT  ('4-Filo  »'.I4.':  '.i> 
CALL  NAME (ERR .FILNAM) 

IF  (ERR  .EQ.  .TRUE.)  GOTO  210 
230  CALL  TRANSF< FILNAM) 

CALL  CLOSE ( 2 ) 

RETURN 

END 


SUBROUTINE  NAME (ERR .FILNAM > 

BYTE  ERR .ANSWER ( 7 ) .TEMPU1 > »FILNAM( 15) 

ERR-. FALSE. 

CALL  OETSTR ( 5 . ANSWER  >  6 » ERR > 

IF  (ERR  .£0.  .FALSE.)  GOTO  20 
TYPE  10 

10  FORMAT  <  ' +A  character*  maximum  for  file  namm,'/) 

RETURN 

20  CALL  CONCAT( '0X1 1' .ANSWER. TEMP) 

CALL  C3NCAT ( TEMP . ' . OAT ' , F ILNAM ) 

RETURN 

END 


SUBROUTINE  TRANSF( FILNAM) 
INTEGER  DATA (103) 

BYTE  MLNAM(  13)  ,ANSUER(6)  .ERR 
CALL  -  3IGNC,. FILNAM. 0> 

CALL  GET3TR ( 3 .ANSWER » 4 .ERR) 
READ  (3.10)  (DATA(I). 1-1.103) 
10  FORMAT  (2513) 

CALL  Pin  STR< 2. ANSWER. 0) 

URITE  (2.10)  (DATA(I) .1-1,103) 
CALL  CLOSE (3) 

RETURN 

END 
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APPENOIX  C.  SOME  THOUGHTS  ON  FITTING  THEORY  TO  TIME -OF- FLIGHT  DATA 

Early  time-of-fl ight  work  (eg.  [45])  was  able  to  extract  only  a  few 
moments  of  a  velocity  distribution,  in  part,  because  of  the  shutter  func¬ 
tion  deconvolution  problem,  and  also  because  the  higher  moments  magnify 
the  effects  of  even  small  amounts  of  noise  in  the  first  part  of  a  time- 
of-flight  spectrum.  Most  workers  [eg. 46, 61 ,62]  now  fit  Maxwellians  where 
appropriate  (as  we  have  done)  with  varying  degrees  of  mathematical  sophis¬ 
tication  (ranging  from  peak  height  matching  to  non-linear  least-squares 
fitting  to  determine  temperature).  For  nozzle  beams,  tranlating  Maxwel¬ 
lians  generally  work  (as  in  our  data;  see  Figure  15);  others  have  used 
the  sum  of  two  translating  Maxwellians  [64]  or  even  a  translating  Max¬ 
wellian  multiplied  by  a  Hermite  polynomial  expansion  [65]. 

Initially,  it  seemed  attractive  to  devise  a  scheme  which  would  fit 
some  sort  of  expansion  in  terms  of  some  set  of  orthogonal  functions  to 
the  data.  This  gives  a  convenient  analytic  form  which  can  be  arbitrarily 
manipulated  for  comparison  with  theory,  and  in  effect,  smooths  the  data 
providing  a  representation  in  terms  of  a  small  set  of  fitting  parameters. 
In  practice,  there  are  problems. 

Basically,  given  sufficient  parameters  to  adjust,  and  a  functional 
form  which  is  approximately  correct,  one  can  readily  make  a  fit  to  most 
any  set  of  experimental  data.  The  problem  is  to  choose  a  function  whose 
parameters  can  readily  be  fitted  to  the  data  (i.e.,  the  function  should 
be  linear  in  the  fitting  parameters),  and  which  approaches  the  data 
with  just  a  few  terms.  For  velocity  distributions,  a  polynomial  times 
an  exponential  is  the  obvious  choice.  Hermite  polynomials  are  attractive, 
because  the  expansion  coefficients  are  closely  related  to  the  moments  of 


the  distribution.  However,  evaluating  these  coefficients  amounts  to 
calculating  moments,  which  is  difficult  to  impossible  beyond  the  first 
three  or  so,  if  there  is  appreciable  noise  in  the  early  steep  high 
velocity  part  of  of  the  time-of-fl ight  curve.  Examples  o *  such  an 
attempt  are  shown  in  Figure  40.  A  better  choice  is  to  select  a  set  of 
polynomials  which  are  explicitly  orthogonal  over  the  data  spacing  with 
respect  to  the  exponential  weighting  function: 


N 


.-2(vu)‘ 


nm 


(C.l) 


i»l 


Such  a  set  of  polynomials  can  be  readily  generated,  although  they 

★ 

must  be  regenerated  for  each  particular  data  spacing  .  The  advantage 
of  this  orthogonality  relation  is  that  the  matrix  inversion  necessary 
to  calculate  the  expansion  coefficients  becomes  trivial  (see  below), 
and  successive  terms  in  the  expansion  can  be  added  without  recalculating 
previous  ones.  If  we  choose  to  set  the  (somewhat  arbitrary)  translational 
velocity,  u,  in  the  exponential  equal  to  zero,  we  have  an  expansion  in 
time  domain  of  the  form: 


£ 


(C.2) 


We  determine  the  by  least  squares  fitting;  i.e., 
for: 

Cf<vj>  •  ^  8,Pi(vj)e-vj2]2 


we  require  a  minimum 


(C.3) 


where  we  take  only  m  terms  of  the  expansion  and  use  v»l/t  for  convenience. 
This  gives  a  set  of  fth-1  equations: 


*  These  polynomials  are  related  to  the  Gram  polynomials  which  have  a 
constant  data  spacing  and  no  weighting  function  [66]. 


114 


For  moderately  noisy  data,  the  expansion  generated  in  this  way  behaves 
well  to  about  5  terms.  For  better  data,  about  10  terms  can  be  obtained 
(see  Figure  41).  After  these  terms,  the  problem  of  magnifying  the  effect 
of  small  amounts  of  noise  again  becomes  important,  and  the  fitted  curve 
starts  to  diverge  from  the  data. 


8eing  a  purely  mathematical  construct,  this  kind  of  curve  fitting  or 


smoothing  tends  to  obscure  rather  than  illuminate  the  physics  of  the 
situation.  The  exponential  contains  important  physical  information  (Mach 
number,  mass,  temperature) ,  which  must  be  chosen  arbitrarily.  In  the 
present  experiments,  many  if  not  most  time-of-fl ight  spectra  seem  to 
consist  of  two  to  four  components  ("fast  peak",  "Maxwellian  peak",  two 
"slow  peaks";  see  Section  3. 4. 2. 2).  What  best  distinguishes  these  com¬ 
ponents  is  the  exponential,  not  some  polynomial.  Especially  in  the  case 
of  incident  beam  data  of  the  present  experiments  where  the  presence  of 
particles  of  widely  varying  mass  produce  a  distinctly  non-Maxwell ian 
spectrum,  there  seems  to  be  nothing  to  be  gained  by  finding  an  arbitrary 
analytic  function  to  fit  the  data. 

A  final  note  of  caution  is  in  order  for  the  interpretation  of  time- 
of-fl  ight  spectra.  Often  the  peaks  are  sufficiently  broad  and  noisy  that 
any  function  of  approximately  the  same  shape  and  enough  (even  3  or  4) 
fitting  parameters  can  be  made  to  fit.  This  does  not,  of  course,  give 
any  indication  of  a  reasonable  theory.  One  must  either  reduce  the  number 
of  fitting  parameters  to  one  or  two  by  independently  measuring  the  others, 
or  at  very  least,  show  that  calculated  fitting  parameters  behave  in 
predictable  ways  for  different  data. 


APPENDIX  0.  SAMPLE  EXPERIMENTAL  PROCEDURE 

Since  the  data  have  been  presented  in  this  document  largely  in  a 
processed  form,  we  present  here  an  example  of  a  typical  experimental 
session.  The  details  obviously  depend  on  what  particular  experiment  is 
being  run,  but  an  example  can  still  give  a  feel  for  the  sorts  of  pro¬ 
cedures  involved.  We  consider  an  experiment  to  survey  a  range  of  Incident 
and  scattering  angles  with  the  time-of-fl ight  spectrometer  for  constant 
beam  and  surface  conditions  (cleaned  silicon  target). 

Starting  from  standby  status  with  the  chamber  evacuated  to  5  x  10"®Pa, 
the  start-up  procedure  is  as  follows:  The  target  is  heated  to  130°  to 
drive  off  adsorbed  material.  The  target  is  lowered  out  of  the  path  of 
the  beam.  The  liquid  nitrogen  cooling  systems  are  started  (valves  are 
opened  on  a  pressurized  supply).  The  source  boiler  is  refilled  with 
distilled  water,  if  necessary,  and  its  heater  is  turned  on.  When  all  the 
liquid-nitrogen-cooled  surfaces  are  near  77°  K.  (by  thermocouple  readings), 
and  the  boiler  is  up  to  the  desired  temperature,  the  source  valve  is 
opened.  The  source  temperature  controller  is  then  activated  and  a  nozzle 
temperature  is  selected.  The  target  temperature  is  adjusted  to  the 
experimental  value.  These  procedures  generally  take  about  an  hour. 

The  time-of-fl ight  equipment  is  turned  on:  ionizer,  chopper,  electron 
multiplier,  amplifiers.  Waveform  Eductor,  oscilloscope,  and  computer.  An 
appropriate  chopper  period  is  selected  (just  long  enough  so  that  there  Is 
no  overlap  between  consecutive  spectra  after  deconvolution  [52])  and  gains 
are  adjusted,  and  a  delay  is  set  on  the  Waveform  Eductor  (typically,  so 
that  the  3  channel  dead  time  between  sweeps  spans  the  trigger  pulse  which 
may  otherwise  be  picked  up  as  signal).  Before  bringing  the  target  into 
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position,  a  spectrum  of  the  incident  beam  is  recorded.  (The  chopper  and 
ionizer  are  positioned  to  look  at  the  beam.)  When  a  satisfactory  spectrum 
can  be  seen  on  the  oscilloscope,  it  is  read  into  the  computer,  deconvoluted , 
stored,  and  copied  onto  paper. 

The  target  is  now  brought  into  position.  For  cold  surface  work,  it 
is  exposed  to  a  normally  incident  beam  for  15  minutes  before  proceeding. 

Then  the  target  is  set  for  the  first  incident  angle  (say  75°),  and  the 
time-of-fl ight  apparatus  is  positioned.  For  a  broad  survey  a  spectrum 
might  be  recorded  at  10°  increments  in  both  scattering  and  incident  angles. 
For  each  spectrum,  recorded  information  might  include  time,  run  number, 
angles,  period,  temperatures  of  target,  source,  and  boiler,  chamber 
pressures,  and  any  changes  in  other  settings.  The  spectrum  is  read  in 
and  processed  as  for  the  incident  beam.  Spectra  can  typically  be  recorded 
at  3-4  minute  intervals.  This  allows  the  Waveform  Educator  to  average 
to  its  noise  limit.  A  survey  of  the  sort  we  describe  might  take  about 
5  hours  total. 
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APPENDIX  E.  EFFECTS  OF  MASS  SPECTROMETER  AMPLITUDE  CORRECTIONS  ON  DETER¬ 
MINATION  OF  SURFACE  AND  SULK  ENERGIES  OF  CLUSTERS 


Let  N  be  the  actual  number  of  clusters  of  size  n  and  A  be  experi 
n  n 

mental  mass  spectrum  peak  heights.  Neglecting  the  possibility  of  frag¬ 
mentation  on  ionization,  we  relate  An  to  Nn  by  applying  the  corrections 
discussed  in  Section  3.1.1: 


(E.l) 


for  18n+l  >m^.  The  quadrupole  transmission  is  unity  for  18n+l  <mth  (the 
last  factor  in  equation  E.l  is  set  equal  to  1  for  18n+l  <mth).  md  is  the 
decay  constant  for  the  transmission  for  18n+l  >mth-  Divide  by  ^Njo, 
where  a  is  a  number  obtained  by  setting  n*l  in  equation  E.l,  and  take 
the  logarithm: 


A  N  f3/7  n  +  0.4’ 

lnA7  *  1niT  ■  1na  +  M/TOTH  +}sln28  -  2.3  (18n+l-mth)/md 

(E.2) 

Nn  * 

Now  set  ln~-  equal  to  -A  n'3  +  Bn  and  rearrange  terms: 

’"st  *  " +  *  t*1"28 

-  lno  +  2.3  (mth  -  1) /md  ]  (E.3) 


The  expression  in  brackets  is  constant  and  does  not  affect  the  determina¬ 
tion  of  other  fitting  coefficients.  The  In  term  is  seen  to  be  negligible 
if  A  and  8  are  of  at  least  order  unity.  It  could  be  included  for  a  more 
precise  fitting.  What  remains  is  a  small  offset  in  8  (since  md  is  typi¬ 
cally  >  150).  There  is  no  error  in  A  as  long  as  the  In  term  can  be  neg¬ 
lected. 


We  can  thus  fit  to  the  A  's  without  correction  and  make  only  small 

n 


errors  in  the  determination  of  A  and  8,  and  fairly  crude  estimates  of  the 
correction  factors  allow  quite  precise  determinations  of  A  and  8.  In 
retrospect  this  makes  sense.  Since  we  are  interested  in  the  logarithm 
of  the  amplitudes,  only  corrections  which  go  exponentially  with  n  are 
important,  and  if  that  exponential  can  be  made  to  decay  much  slower  than 
the  Nn  themselves,  almost  no  correction  is  necessary. 


120 


FIGURE  1.  SCHEMATIC  OIAGRAM  OF  APPARATUS 

A  -  Temperature  controlled  boiler;  B  -  Shielded  heat  pipe;  C  -  Nozzle; 
D  -  Oven  around  nozzle;  E  -  Liquid  nitrogen-cooled  inner  wall; 

F  -  Collimating  chamber;  G  -  Collimating  plate;  H  -  Beam  flag; 

I  -  Target;  J  -  Chopper;  K  -  Ionizer  and  electron  multiplier; 

L  -  Quadrupole  mass  spectrometer  (shown  in  the  two  positions  used); 

M  -  Main  vacuum  chamber;  N  -  Liquid  nitrogen -coo led  flow  through 
enclosure. 


1 

i 
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FIGURE  2.  EARLY  NOZZLE  SOURCE  ASSEMBLY 

A  -  Vacuum  pump;  B  -  Gas  supply  manifold;  C  -  Needle  valve; 

0  -  Vacuum  wall;  E  -  Water  reservoir  with  external  heater; 

F  -  Motor-driven  needle  valve;  G  -  Nozzle;  H  -  Collimator; 

J  -  Liquid  nitrogen  reservoir;  K  -  Capacitance  manometer  (pressure  gauge). 
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CONTINUE 
DO  20 t  I»2,104 


109 


FIGURE  3.  FLOW  THROUGH  CYLINDRICAL  NOZZLE 

The  growing  boundary  layer  through  the  cylindrical  channel  constricts 
the  flow  to  give  a  converging-diverging  flow  pattern  as  in  "supersonic 
nozzles. 


FIGURE  4.  FINAL  SOURCE  ASSEMBLY 

A  -  Stainless  steel  boiler;  B  -  1000W  cartridge  heater;  C  -  Thermocouple; 
D  -  Fill  tube;  E  -  Distilled  water  supply  bottle;  F  -  Overpressure 
relief  valve;  G  -  Collimating  chamber;  H  -  Shielded  heat  pipe; 

J  -  100W  source  oven;  K  -  Nozzle  (with  another  thermocouple,  not  shown); 

L  -  Collimator;  M  -  Beam  flag  in  main  vacuum  chamber;  N  -  Liquid 
nitrogen-cooled  inner  wall. 
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FIGURE  5.  TARGET  ASSEMBLY 


Not  to  scale.  Target  -2x  actual  size,  reservoir  -1/5  actual  size. 

A  -  Vacuum  wall;  B  -  Liquid  nitrogen  reservoir;  C  -  Flexible  (stainless 
steel  bellows)  tubing;  D  -  Liquid  nitrogen -coo led  copper  block; 

E  -  Alumina  spacers;  F  -  Copper  plates;  G  -  Mica  sheets;  H  -  Stainless 
steel  screws;  J  -  Target;  K  -  Thermocouple;  L  -  Nichrome  wire; 

M  -  Thermocouple. 


FIGURE  6.  IMPORTANT  DIMENSIONS  OF  BEAM  SYSTEM,  TARGET,  AND  TIME -OF -FLIGHT 
APPARATUS 

All  dimensions  in  cm.  Dimensions  next  to  labels  qive  diameter  of  opening. 
A  -  Mass  spectrometer  ionizer  positioned  for  scattering  measurements. 

B  -  Mass  spectrometer  ionizer  positioned  for  beam  measurements. 

Ionization  region  is  tapered  from  .47  cm  to  .64  cm. 
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FIGURE  7.  TIME-OF-FLIGHT  SIGNAL  PROCESSING 
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FIGURE  8.  EXAMPLES  OF  MASS-FILTERED  TIME-OF-FLIGHT  SPECTRA  FOR  BACKGROUND 
SPECIES 

(a)  H,0  beam,  P0*l . 1  MPa,  T0»458°  K.  Mass  filter  set  at  282  AMU,  flight 
path  /0.6  cm*,  period  1.75  ms*.  Note  that,  although  examde  (a)  does  not 
show  a  dip  below  zero,  the  shape  of  the  curve  is  not  typical  of  a  velocity 
distribution;  the  steep  portion  is  on  the  wrong  side  of  the  peak. 

(b-d)  H20  beam,  Pq*830  kPa,  Tc*466°  K,  scattering  from  CuS  at  241°  K  (e^e 
*72°).  Period  1.92  ms*,  flight  path  41.9  cm*,  mass  filter  set  as  marked.  s 


*In  these  and  all  tlme-of-f! ight  spectra,  time  increases  to  the  right, 
velocity  to  the  left.  They  represent  particle  density  in  the  ionizer 
versus  time.  To  obtain  velocity,  divide  the  flight  path  by  the  time. 

There  are  103  channels;  each  channel  represents  a  time  of  1/103  x  period. 
The  time  origin  is  at  or  near  the  beginning  of  the  trace  and  Is  marked  on 
traces  where  a  theoretical  curve  has  been  fitted.  All  soectra  are  plotted 
at  the  same  amplitude  (i.e.  relative  amplitude  Information  is  not  shown). 
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FIGURE  9.  MASS-FILTERED  TIME -OF -FLIGHT  SPECTRA  FOR  MIXED  Ar/He  BEAM 
WITH  N2/O2/H2O/CO2  BACKGROUND 

952  He  +  52  A r  beam,  P0»7.1  kPa,  T<j*284°  K.  (a)  spectrum  without  mass¬ 
filtering,  period  1.12  ms*,  flight  path  20.6  cm*,  (b-g)  mass-filtered 
spectra,  period  2.92  ms*,  flight  path  70.6  cm*,  mass  filter  set  as  indi¬ 
cated.  Examples  (e)  and  (f)  might  be  mistaken  for  normal  tlme-of-fl ight 
spectra,  but  note  that  the  rising  slope  Is  not  as  steep  as  In  (b)  and  (c).t 

*See  footnote  for  Figure  8. 

+The  spectra  could  also  be  Interpreted  as  an  indication  of  air  contamina¬ 
tion  In  the  source. 
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FIGURE  10.  TIME -OF -FLIGHT  SPECTRA  OF  ARGON  BEAM 
Tq«308°  K,  P3  as  marked,  period  1.62  ms*,  flight  path  20.6 
*See  footnote  for  Figure  8. 


a.  3.2  fePa 


d.  270  kPa 


b.  16  kPa 


e.  740  kPa 


c.  91  kPa 


f.  1150  kPa 


FIGURE  11.  TIME  OF  FLIGHT  SPECTRA  OF  WATER  BEAM  FOR  VARYING  SOURCE 
PRESSURE 

T0*456®  K,  PQ  as  marked,  no.  mass  filtering,  periods  (a,b)  1.07  ms, 
Cc-f)  .81  ms*,  flight  path  20.6  cm*.  Note  esoecially  the  appearance 
of  an  increasingly  prominent  low  velocity  tail  at  high  pressure. 


*See  footnote  for  Figure  8. 


1.  6.  11.16,21.26.33.36.41.46.51.56. 
Cluat«r  Sixm  (n) 


FIGURE  12.  MASS  SPECTRUM  OF  A  HIGHLY-CLUSTERED  WATER  BEAM 

Pq*M  MPa,  T0»458°  K.  See  Section  4.1  for  an  explanation  of  fitted  curve. 
Chart  recorder  output  illustrates  peaks  for  half-integer  cluster  sizes. 
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FIGURE  13.  STAGNATION  GAUGE  INTENSITY  VERSUS  SOURCE  PRESSURE 

Range  of  points  shown  cover  a  temperature  range  of  400-500°K.  The 
variation  of  Intensity  with  source  temperature  is  relatively  small. 
Higher  temperature  results  in  lower  intensity.  Thus  the  lowest  points 
at  any  fixed  source  pressure  In  the  graph  are  for  500°K,  and  the  highest 
are  for  400°K. 


133 


FIGURE  14.  TIME-OF-FLIGHT  PEAK  AMPLITUDE  VERSUS  BEAM  INTENSITY 

Intensity  at  right  corresponds  to  source  conditions  where  beam  appears 
highly  clustered  as  indicated  by  the  low  velocity  tail  in  the  time-of- 
flight  spectrum  (Figure  11)  (Poal.l  MPa). 


FIGURE  15.  MASS-FILTERED  TIME-OF-FLIGHT  SPECTRA  OF  A  CLUSTERED  BEAM 

H20  beam,  P0*1.1  MPa,  T0*458°  K,  periods:  (a)  .98  ms,  (b-t)  1.75  ms, 
(u-ff)  1.46  ms*,  flight  paths:  (a)  20.6  cm,  (b-ff)  70,6  cm*,  (a)  no 
mass  filtering,  (b-ff)  mass  filter  set  as  Indicated  (cluster  order  also 
Indicated),  (h-ff)  translating  Maxwellian  shown  fitted  to  data  using 
Indicated  time  origin  (obtained  with  an  assumed  delay  of  157  ys). 

(a-t)  data  of  25-Jan-80,  (u-ff)  data  of  31-Uan-80.  (ff)  uses  a  half- 
Integer  mass  peak  (see  Figure  12). 


*See  footnote  for  Figure  8. 
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FIGURE 


FIGURE  15  (.concluded). 
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FIGURE  16.  PEAK  AMPLITUDE  VERSUS  CLUSTER  SIZE  FROM  MASS-FILTERED 
TIME-OF-FLIGHT  SPECTRA  OF  BEAM 

Taken  from  data  of  Figure  15.  Two  short  ranges  shown  here.  Squares 

are  peak  amplitudes,  triangles  are  amplitudes  from  mass  spectra  (Figure  12). 

Amplitudes  are  scaled  so  that  the  two  amplitudes  match  at  n=4  and  n=»l  1 . 
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FIGURE  17.  TRANSLATIONAL  VELOCITY  VERSUS  CLUSTER  SIZE  FROM  MASS- 
FILTERED  TIME-OF-FLIGHT  SPECTRA  OF  BEAM 

Line  after  n«16  divides  data  from  separate  days.  Data  obtained  from 

fitted  curves  In  Figure  15.  Plotted  is  u  from  v2e"m^v"u^  /^kT 

Note  that  there  are  no  points  for  n»l-*-3  since  it  was  not  possible  to 
fit  a  translating  Maxwellian  to  the  time-of- flight  spectrum  for  these 
cluster  sizes  (see  Figure  15). 


*  Actually,  since  the  data  are  time-of-fl ight,  the  fitted  curve  is  a 
time  distribution  function  determined  as  follows: 

f(v)dv«f(^)£jdt  -  (k)4e-*([L/t]-u)2/2kTdt 
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FIGURE  18.  TRANSLATIONAL  TEMPERATURE  VERSUS  CLUSTER  SIZE  FROM  MASS- 
FILTERED  TIME-OF-FLIGHT  SPECTRA  OF  BEAM 

Line  after  n*16  divides  data  from  separate  days.  Data  obtained  from 

fitted  curves  in  Figure  15.  Plotted  is  T  from  v2e"m^v'u^2^^*  where 
m  was  taken  as  the  monomer  mass  throughout. 

Note  that  there  are  no  points  for  n»l-*3  since  it  was  not  possible  to 
fit  a  translating  Maxwellian  to  the  time-of-fl ight  spectrum  for  these 
cluster  sizes  (see  Figure  15). 

*  See  footnote  for  Figure  17. 
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FIGURE  19.  MASS  18  TIME -OP-FLIGHT  SPECTRA  OF  BEAM  VERSUS  SOURCE 
TEMPERATURE  (T0) 

H?0  beam,  (a-g)  Pn*680  kPa,  (h)  P0*330  kPa,  T0  as  marked,  flight  Dath 
70.6  cm*,  period  1.45  ms*.  Note  gradual  dlsapearance  of  slower  peak. 
The  mass  spectrum  for  example  (h)  still  showed  clusters  to  n*20. 


*See  footnote  for  Figure  8. 
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Emiaaion  Angla  (dagraaa) 


FIGURE  20.  PEAK  AMPLITUDE  VERSUS  EMISSION  ANGLE  FOR  EVAPORATION 
FROM  COVELLITE 

No  beam,  CuS  (1000)  target.  Frost-like  ice  layer  grown  from  beam  at 
T  *180°  K.  Oata  shown  for  T  *204°  K. 


FIGURE  21.  TIME-OF-FLIGHT  SPECTRA  VERSUS  EMISSION  ANGLE  FOR  EVAPORATION 
FROM  COVELLITE 

Same  conditions  as  for  Figure  20.  Flight  path  20.6  cm*,  period  1.61  ms*. 
Fitted  curves  are  static  Maxwell ians  at  204°  K.  Emission  angles  are 
indicated. 


*$ee  footnote  for  Figure  8. 
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FIGURE  22.  TIME-OF-FLIGHT  SPECTRA  FOR  EVAPORATION  FROM  SINGLE  CRYSTAL  ICE 

T$  and  e$  as  marked.  Mass-filtered  spectra,  mass  filter  set  at  13  AMU. 

Flight  path  41.9  cm*,  period  2.96  ms*.  Fitted  curves  are  static  Maxwel- 
Hans  at  T$  using  a  delay  of  90  us. 


*See  footnote  for  Figure  8. 
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FIGURE  24.  (concluded) 

b.  Amplitude  of  Maxwellian  component  of  scattered  stream.  H20  beam, 
To^SS®  K,  P0«560  kPa,  scattering  from  SI  at  223°  K.  For  the  colllmatlon 
geometry  of  the  system,  a  cosine  distribution  gives  a  constant  amplitude. 
The  apparent  variation  with  ei  Is  a  consequence  of  a  drift  in  the  gain  of 
the  electron  multiplier  for  which  corrections  have  not  been  made. 
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a.  chopping  scattered  stream 


FIGURE  25.  COMPARISON  OF  TIME-OF-FUGHT  SPECTRA  CHOPPING  INCIDENT 
BEAM  VERSUS  CHOPPING  SCATTERED  STREAM  WHEN  MOST  OF  BEAM 
STICKS  TO  TARGET 

HjO  beam  Pq»880,  T0*556°  K  scattering  from  SI  at  212°  K,  e,»65®,  9  *75°, 
flight  path  20.6  cm*  (for  chopping  scattered  stream),  period  1.76  msec*. 
Chopping  position  as  marked. 


♦See  footnote  for  Figure  8. 
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INCIDENT 


FIGURE  26.  FAST  PEAK  AMPLITUDES 

Source  conditions:  To*455°  K,  Pqs560  kPa ;  Target:  Si(m), 

T  *223°  K.  Dashed  lines  connect  points  of  constant  "turning  angle" 
(i,  +  0  *  constant).  The  amplitudes  are  normalized  by  dividing  by 
th6  associated  Maxwellian  peak  amplitude.  This  eliminates  any 
drift  effects  and  automatically  corrects  for  the  flux  variation 
with  incident  angle. 
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FIGURE  27.  VELOCITIES  AT  MAXIMUM  AMPLITUDES  FOR  THE  FAST  PEAK 

Same  conditions  as  for  Figure  26.  Data  for  only  a  few  incident  angles 
are  shown  for  clarity.  Note  that  the  peak  velocity  is  also  constant 
for  constant  "turning  angle"  (see  Figure  26). 


FIGURE  28.  EXAMPLES  OF  TIME-OF-FLIGHT  SPECTRA  FOR  THE  SPECULARLY 
DIRECTED  SLOW  PEAK 

H20  beam  P0»57O  kPa,  Tg«456°  K  scattering  from  Si  at  241®  K.  Flight 
path  20.6  cm*,  period  1.58  ms*.  There  are  three  peaks  present  in  each 
spectrum:  a  fast  peak  and  a  Maxwellian  peak  (unresolved)  and  a  slow 
peak.  Note  the  large  Intensity  of  the  slow  peak  and  its  variations 
with  angle.  All  four  spectra  are  scaled  to  the  same  height.  The 
actual  changes  in  the  amplitudes  of  the  slow  peak  with  e,  can  be  realized 
by  noting  that  the  peaks  on  the  left  (Maxwellian  and  fast  peaks)  have 
nearly  consistent  amplitudes  with  changing  8. 


♦See  Figure  8. 
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FIGURE  29.  AMPLITUDES  OF  THE  SPECULARLY  DIRECTED  SLOW  PEAK 

Same  conditions  as  for  data  in  Figure  28.  Amplitudes  are  normalized 
by  dividing  by  the  Maxwellian  peak  amplitude  as  in  Figure  26.  Note 
specular  reemission  angles. 
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FIGURE  30.  STAGNATION  GAUGE  AMPLITUDES  FOR  SPECULARLY  OIRECTED 
SLOW  PEAK 

at  241°  K.  Stagnation 
ition  with  the  same 
colllmation.  Since  the  stagnation  gauge  measures  flux,  amplitudes  are 
divided  by  the  average  velocity  obtained  from  time-of-fl ight  data  to 
give  a  density  which  can  be  compared  to  the  time-of-fl Ight  peak 
amplitude  data  of  Figure  29. 


H20  beam  P0a5OO  kPa,  T0*455°  K  scattering  from  S. 
gauge  located  at  usual  time-of-fl ight  ionizer  pos 
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FIGURE  31  FITTING  OF  TRANSLATING  MAXWELLIAN  TO  SPECULARLY 
OIRECTED  SLOW  PEAK 


1,0  beam  Pn»560  kPa,  T0»455°  K  scattering  from  Si  at  241°  K,  flight 
>ath  20?6  cm*,  period  ?.58  ms*.  Angles  as  marked.  Maxwellian 

:omponent  subtracted  out.  Fitted  Note 

nonomer  mass,  Mach  number  2.2,  temperatures  (a)6  K,  (»)J5  K-  Note 
that  a  Mach  number  of  2.2  also  provides  an  approximate  - fit  to  the 
time-of-fl ight  spectrum  (without  mass  filtering)  of  the  incident 

beam  (Figure  11). 


♦See  Figure  8. 
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FIGURE  32.  PEAK  POSITIONS  IN  TIME-OF-FLIGHT  SPACE  FOR  THE  SPECULARLY 
DIRECTED  SLOW  PEAK 

Same  conditions  as  for  data  in  Figures  28  and  29.  The  ordinate 
represents  the  shifts  noted  In  Figure  28.  The  fast  peak  and 
Maxwellian  peak  positions  are  marked  for  comparison. 
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FIGURE  33.  PEAK  AMPLITUDES  VERSUS  BEAM  INTENSITY 

H20  beam  P0*25(M120  kPa,  T0*471°  K  scattering  from  Si  at  239°  K. 
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FIGURE  34.  AMPLITUDES  OF  THE  NON-SPECULAR  SLOW  PEAK 

H20  beam  P0*960  kPa,  T0»452°  K  scattering  from  Pt  at  293°  K.  Data 
for  only  a  few  incident  angles  are  shown  for  clarity.  Amplitudes 
are  normalized  by  dividing  by  the  Maxwellian  peak  amplitude  as  in 
Figures  26  and  29. 
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FIGURE  35.  FITTING  OF  TRANSLATING  MAXWELLIAN  TO  NON-SPECULAR  SLOW  PEAK 

H«0  beam,  P0«950  kPa,  T0*452°  K,  scattering  from  Pt  at  291°  K,  flight  path 
20.6  an*,  period  1.58  ms*,  angles  as  marked.  Maxwellian  component  subtracted 
out.  Fitted  with  translating  Maxwell  Ians,  monomer  mass,  Mach  number  1.6, 
temperatures  (a)  80°  K,  (b)  35°  K.  Compare  Figure  31. 


♦See  footnote  for  Figure  8. 
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FIGURE  36.  PEAK  POSITIONS  FOR  THE  NON-SPECULAR  SLOW  PEAK 

Same  conditions  as  for  Figure  34.  Plot  Is  as  a  function  of  incident 
angle,  since  data  can  be  obtained  for  only  a  few  scattering  angles. 

The  fast  peak  and  Maxwellian  peak  positions  are  marked  for  comparison. 
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a.  18  AMU 


b.19  AMU 


c.  37  AMU 


FIGURE  37.  VELOCITY  SPECTRA  OF  THE  ONLY  THREE  WATER  SPECIES  FOUNO  IN 
SCATTERING  FROM  SI 

H20  beam,  To-4520  K,  Pn«1.Q  MPa,  scattering  from  SI  at  29V*  K.  flight 
path  41.9  cm*,  period  2.99  ms*,  e^»67°,  e$*79*.  No  spectra  could  be 
generated  for  mass  filter  set  higher  than  37  AMU . 

*See  footnote  for  Figure  8. 
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FIGURE  39.  COMPONENTS  OF  FAST  PEAK  VELOCITY 

Taken  from  data  shown  in  Figures  26  and  27.  Data  are  normalized  by 
dividing  by  the  Incident  velocity.  For  reference,  the  components  of 
the  Incident  velocity  are  shown  as  a  smooth  curve.  Note  that  all 
points  above  this  curve  represent  conditions  where  there  has  been 
a  gain  In  that  particular  component.  This  page  shows  normal 
components;  tangential  components  are  on  the  next  page. 


FIGURE  40.  EXAMPLES  OF  HERMITE  POLYNOMIAL  FITTING  TO  DATA 

thV“f!r  tern»  OF  the  expansion  used.  For  the 
good  data  of  example  (a),  six  terms  gives  a  good  fit;  more 

Introduees  extraneous  peaks.  No  satisfactory  fit  Is  achieved 
In  the  noisier  data  of  example  (b).  acmeveo 


FIGURE  41.  EXAMPLES  OF  MODIFIED  GRAM  POLYNOMIAL  FITTING  TO  DATA 


Numbers  Indicate  the  number  of  terms  of  the  expansion  used.  For  the 
"good"  data  of  example  (a),  10  terms  give  an  excellent  fit  to  the 
data.  For  noisier  data  (the  same  spectrum  used  In  Figure  40b),  4  or 
5  terms  give  a  good  fit.  For  6  or  more  terms,  there  are  again 
extraneous  peaks. 
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